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Interactions 


Between Wholly Laminar or 


Wholly Turbulent Boundary Layers and 
Shock Waves Strong Enough to Cause 
Separation’ 


National Physical Laboratory, 


SUMMARY 


discussion is presented covering both wholly 


transition occurs 


A qualitative 
laminar and wholly turbulent cases (i.e., 
upstream of the region of interaction). 
for the pressures at with 
boundary layers. Semiempirical relations permitting the rapid 
calculation of the overall characteristics of turbulent interactions 


where 
Relations are 
and 


obtained 


separation laminar turbulent 


are also obtained 


List OF SYMBOLS 


distance from leading edge measured parallel to wall 
distance from wall measured normal to wail 

mean x-component velocity in boundary layer 
x-component velocity at edge of boundary layer 
fluctuating x-component velocity 

mean y-component velocity 

fluctuating y-component velocity 

mean pressure 

mean density 

density at edge of boundary layer 

fluctuating density 

viscosity 

viscosity at edge of boundary layer 

heat at constant to specific 


ratio of specific pressure 


at constant issumed to equal 1.4 


heat volume :: 
viscosity temperature relation index (ua7@, where T is 


absolute temperature ) 
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above was carried out in the Aerody- 
National Physical Laboratory, 


Aeronautical Research 


work described 
and it is 


GADDt 


Teddington, England 


= Prandtl] Number (uc,/k, where cp is specific heat at 
constant pressure and k& is thermal conductivity) 
= Mach Number at edge of boundary layer 
Reynolds Number, pio l’ox/ m0 
boundary-layer total thickness 
6 for zero pressure gradient laminar boundary layer 


; F pu 
displacement thickness 1 — — | dy 
0 pil 


6* for zero pressure gradient laminar boundary layer 


. pu u 
momentum thickness -tl- -} dy 
0 pil l 


parameter in assumed velocity for laminar 
boundary layer [see Eqs. (4)] 
as X 
t = friction stress 
y angle between edge of boundary layer 


profile 


and local mean 
flow stream lines 
The suffices F, 0, R, s refer to the positions F, 0, R, 
and 3 [see Section (2)]. 


s of Figs. 2 


(1) INTRODUCTION 


r I HIS PAPER IS CONCERNED WITH the two-dimensional 
interaction between the boundary layer on a flat 


plate and an incident oblique shock wave strong enough 
The configuration envisaged is as 
1, and it is assumed that the boundary layer is 
i.e., that 

interac- 


to cause separation. 
in Fig. 
either wholly laminar or wholly turbulent 
transition does not occur in the region of 
tion. 

Wholly laminar cases are physically similar to wholly 
turbulent ones, and in Section (2) a qualitative interpre- 
tation for both conditions is presented. 

In Section (3) a relation is obtained for the pressure 


at the separation point with laminar boundary layers. 


729 
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Shapes. of pressure distribution at the wall typical 
with (a) laminar (b) turbulent boundary layers. 














By a similar but much cruder method, an expression for 
the pressure at separation in wholly turbulent cases is 
derived in Section (4). Semiempirical relations are also 
obtained permitting the rapid calculation of the overall 
effects produced in a turbulent boundary layer when 
the shock is strong enough to cause a large extent of 
separated flow. 
tained in Section (3) for wholly laminar cases with 


Similar relations could have been ob- 


separation, but this was not considered worth while, since 
the obtaining of such cases experimentally is fairly dif- 
ficult so that they are in a sense more academic. More- 
over, semiempirical methods can only make a limited 
contribution to the theoretical understanding of any 
problem, so that their application to laminar cases, for 
which an adequate theoretical treatment may well be 
feasible, seems less justifiable than applying them to 
turbulent cases, where some admixture of empiricism is 
probably almost inevitable. 
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The method of Section (4) for turbulent boundary 
layers do not explicitly deal with all the relevant fae. 
tors. Thus, while it is shown in the qualitative dis. 
cussion of Section (2) that the turning of the flow toward 
the wall which occurs in the immediate vicinity of the 
shock is an important factor, the methods of Section (4 
take no explicit account of it. They account for it im. 
plicitly in an empirical constant. Such a procedure js 
possible because interactions between turbulent bound. 
ary layers on flat plates and incident oblique shock 
waves represent a fairly restricted class of conditions 
even with varying shock strength and Mach Numbers, 


(2) QUALITATIVE DISCUSSION 


The following discussion is in part based on the 
analysis of Section (3) below and also partly on reference 
1. It is consistent with the experimental results of 
reference 2 but not entirely so with those of reference 3.} 


When a sufficiently strong oblique shock wave is in- 
cident uport the boundary layer on a flat plate, as in Fig. 
the 
boundary layer to thicken and separate upstream of it. 


2, it causes both pressure to increase and_ the 
The thickening of the boundary layer generates a band 
of compression waves that determine the pressure dis- 
tribution acting on the boundary layer upstream of the 
shock, and this pressure distribution in turn governs the 
rate of thickening of the boundary layer. The two proc- 
esses must adjust themselves to be in equilibrium, so 
that the pressure distribution upstream of a point J; 
where a shock impinges on a boundary layer would pre- 
sumably remain unaltered if the shock were increased 
in strength and simultaneously moved downstream to 
some point /, such that the separation point did not 
move. Hence, the ratio of the pressure at separation to 
the undisturbed free stream pressure should be a func- 
tion of Mach Number and Reynolds Number only, 
independent of shock strength. 


The shock is reflected locally where it strikes the 
boundary layer as an expansion just as it would be if it 
impinged on a constant pressure boundary. This is be- 
cause without the expansion there would be a discon- 
tinuous jump in pressure which the separated boundary 


layer could not withstand. The flow at the edge of the 








boundary is turned in toward the wall by the shock- | 


expansion system, downstream of which accordingly 


—— 


there is a further pressure rise as the flow turns back | 


approximately parallel to the wall. 

+ When comparing these two sets of experimental results, one 
should bear in mind the fact that the boundary layer thickness as 
the same as the thickness of the 
The latter 
thickness is probably roughly the same as the thickness of the 


defined in reference 3 is not 
boundary layer as it appears in a schlieren photograph. 


zone in which the friction tangential stresses on surfaces normal 
to the local isobars are appreciable, and it is this zone which is de 
fined as the boundary layer in the present paper. 
boundary layer so defined is almost coincident with the local mean 
flow stream lines 
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INTERACTIONS BETWE 

The pressure distribution at the wall is as in Fig. 3. 
There is a relatively steep (for turbulent boundary 
lavers very steep) initial rise of pressure up to the 
separation position s. Downstream of s for some dis- 
tance there is a “dead air’’ region next to the wall.7 
Where this dead air region is thin, it can sustain an ap- 
preciable adverse pressure gradient largely because of 
the friction forces acting on it. Where the dead air 
region is thick, however, the rate of change of friction 
stress normal to the wall must become small in the dead 
air region, and the boundary layer will be unable to 
Now the 
dead air region increases rapidly in thickness down- 
stream of s because the edge of the boundary layer is 


withstand a large adverse pressure gradient. 


deflected away from the wall at a considerable angle. 
Accordingly, the pressure gradient soon falls off down- 
stream of s. At /, the point where the incident shock 
strikes the boundary layer and is reflected locally as an 
expansion, the external flow is turned in toward the 
wall, and the dead air region begins to reduce in thick- 
ness and to become capable of sustaining a larger ad- 
verse pressure gradient. Hence, the pressure gradient 
again begins to rise until the reattachment point R is 
reached. 

In the case of turbulent boundary layers, mass ac- 
celeration effects may play a significant part in the 
equilibrium of the so-called dead air region. Where the 
dead air region is reducing in thickness, as near reattach- 
ment, the mass acceleration forces will balance part of 
the pressure gradient forces, all of which accordingly 
will not have to be met by the friction forces. On the 
other hand, where the dead air region is increasing in 
thickness, the mass acceleration effects will be such that 
friction will have to meet than the pressure 
gradient forces. This may explain why, with turbulent 


more 


boundary layers, the pressure increase between the 
shock and reattachment is much greater than that be- 
tween separation and the shock, whereas there is much 
less difference with laminar layers where mass accelera- 
tion effects in the dead air region are probably less im- 
portant. 

Downstream the R the 
boundary-layer profile begins to approach the zero pres- 
sure gradient form, and the pressure gradient once more 
falls off. The pressure asymptotes to the peak pressure 
br, which is reached shortly downstream of R at F. 
The ratio of this peak pressure to the original free- 
stream pressure is approximately the same as if there 
were no boundary layer present, when the incident 
shock would simply be reflected as another single shock 
as in Fig. 4, a pattern known as “regular reflection.” 


of reattachment point 


t The dead air region is defined as being of thickness t(x), 


ba 
f pu dy = 0 
0 : 


The velocities in the dead air region, though not large, are not 
zero, since there is reversed flow near to the wall and forward flow 


where 


further out. 
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With the help of shock tables it is easy to calculate the 
regular reflection pressure rise. This is approximately 
equal to the actual pressure rise because the waves re- 
flected from the boundary layer must have the overall 
effect of turning the external flow deflected through the 
incident shock back approximately parallel to the wall, 
just as the single reflected shock does when there is no 
boundary layer present, and though the change of 
entropy in the external flow is slightly different in the 
two systems, this makes little difference to the pres- 
sures reached. 

The fact that an increase of shock strength moves the 
separation point upstream may be interpreted physi- 
cally as follows: At J, the position where the shock 
strikes the edge of the boundary layer (see Figs. 2 and 
3), the pressure is little affected by the distance of the 
separation point s upstream of J, provided s/ is not too 
small, because the pressure gradient becomes small be- 
tween s and J some distance downstream of s. Hence, 
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just downstream of /, where the incident shock is re- 
flected locally as an expansion as at a constant pressure 
boundary, the external flow will be more steeply inclined 
toward the wall if the shock strength is increased. With 
increase of shock strength, it is unlikely that the bound- 
ary layer will be able to withstand any increase of pres- 
sure gradient between J and the reattachment point R, 
so that the curvature of the external flow stream lines 
between J and R cannot increase. The boundary-layer 
thickness at R will probably increase with increase of 
shock strength, while the external flow direction at R 
will remain nearly parallel to the wall because the pres- 
sure at Ris near to the maximum. In order, therefore, 
for the external flow to accommodate itself to the 
change of direction between / and R without violating 
the conditions imposed by the boundary layer, the 
thickness of the boundary layer at J must be greater 
when the shock strength is greater, as in Fig. 5. This 
requires the separation point to move farther up- 
stream of the shock. 


The possibility that mass acceleration effects may 
play a more important role in the equilibrium of the 
so-called dead air region with turbulent boundary 
layers than with laminar has been pointed out. Never- 
theless, in broad outline the laminar and turbulent 
cases are physically similar, and the principal difference 
lies in the fact that to produce comparable effects shock 
strengths some five to ten times greater are required for 
turbulent cases than for laminar ones.t Also, of course, 
for attaining wholly laminar flow is 
Number should extremely 


a condition 
that the Reynolds 
low, while it must be relatively high for turbulent 
flow. 


be 


(3) BouNDARY LAYERS LAMINAR OVER THE WHOLE 
REGION OF INTERACTION 


An approximate method of predicting the behavior of 
the laminar boundary layer is developed. It is a variant 
of an unpublished method for incompressible flow due to 
B.S. Stratford. Besides being easier to apply than the 
Pohlhausen method, it is probably more accurate for 
cases where small pressure increases with sharp pressure 
gradients occur after a considerable length with zero 
pressure gradient. The method gives grounds for the 
physical interpretation presented in Section (2) above, 
and by its use a fairly simple expression for the pressure 
ratio at separation is obtained. 


(3.1) The Method of Analysis 


The following approximate method of analysis is 
applicable to cases where a boundary layer, after grow- 
ing for a distance x) with constant pressure fp», ex- 
periences pressure increases small compared with f» in a 
distance not greater than, say, xo/2. 


t The shock strength is defined as the ratio of the pressure rise 
through the shock to the pressure upstream of the shock. 


AERONAUTICAL 


SCIENCES NOVEMBER, 1953 

In the region where the pressure (which is assumed 
not to vary through the thickness of the boundary layer 
in the direction normal to the wallf{) is equal to py, the 
velocity profiles are assumed to be of the form 


= sin (17/2)(y/6.) (1) 
1 and there is no heat transfer to the wall, 


tia=«¢ = 


so that 


a pw0P 


a u-\y> 1 
rf) . ( 7 A 2 us | 
a\y~— I : 
H=poll+ {1 —- Uy? 2 Mo" 


the displacement thickness with zero pressure gradient 


is 
_ 1 1/o 
2 tan 1(° " ) Ms] 


Tv [(v “a 1) eT . M) 
But from reference 4 


6.* = 1.721[1 + 0.693(y — 1) Mi2](x/VR,) (2 


where R, is the Reynolds Number, pil/ox/u. Hence, 


1.721 [1 + 0.693(y — 1)Mg?]x 


er | 1 . 
0.636 tan~! (7 - ) Mo] 


1/ j R, 
[(y — 1)/2]'7 Mo * 
(3) 


Where the pressure differs from po, the approximate 
form of the velocity profiles can be deduced as follows: 
In Fig. 6a the velocity profiles at two points 0 and P are 
sketched. 0 is the point downstream of which the pres- 
sure begins to increase above fp. The points By and B 
are on the same stream line, and B is in the region of the 
point of inflexion of the profile at P. If the pressure 
gradient between 0 and P were very steep, the pressure 
gradient forces acting on the fluid in the outer part of 
the boundary ,layer}would be much greater than the 
viscous forces, which could accordingly be ignored. 
Hence, the outer profile at P could be derived from that 
at 0 by the application of Bernouilli’s equation and 
continuity considerations. Only the total change of 
pressure between 0 and P would be relevant; the way 
in which the pressure varied with x would not matter so 
long as the pressure gradient were everywhere suf- 
ficiently steep. If, however, the pressure gradient be- 
tween 0 and P is not steep, it is not legitimate to ignore 
the viscosity effects, but with a little modification the 
method can take account of these approximately. Fig. 
6b shows the velocity profiles as they would be if the 

t This assumption is probably not seriously in error for laminar 
layers except in the immediate neighborhood of the shock. 
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INTERACTIONS 


pressure remained equal to fy; B’ ison the same stream 
jine as By. (The profile at 0 is, of course, the same as 
in Fig. 6a.) Outside of B in Fig. 6a the profile does not 
differ greatly from that outside of B’ in Fig. 6b, so the 
viscous forces acting between 0 and P on the stream 
lines passing outside of By at 0 are little different in the 
two cases. It seems plausible to treat the effects of 
viscosity and pressure gradient as separable and to as- 
sume that the profile outside of B at the point P, where 
the pressure is p, is approximately the same as if the pro- 
file outside of B’ in Fig. 6b were suddenly subjected to 
the pressure rise Pp — fp. On these assumptions inviscid 
fluid theory indicates that outside of B the profile is of 


‘ (y+1)/27 /y 
“<4 = sin ~! (?) (; - s) 
l 2 \pbo 6. 


where 6, is given by Eq. (3), with x the distance of P 
from the leading edge of the boundary layer, and ) is an 


the form 


arbitrary constant dependent on the distance that the 
“outer” profile at P is displaced from the wall. For the 
complete profile at P, to satisfy the conditions that the 
velocity should be zero at the wall and that the velocity 
and its first two derivatives with respect to y should be 
continuous through the thickness of the boundary layer, 
it is assumed that 


u J Vv Po Gt) /27 ) 
- = %, formX <—- <A 
l é- p 
=7,+ 972, for 0<— < md 
where } (4) 


r, = sin 
2 \po 6. 
yy 3 (y+1)/2y 
r= E —_— | sin : (”) nN 
nn; 2 \po 


(see Fig. 7). The outer edge of the boundary layer is at 


¥y é- =X + (Po prt " 
and the position y/6, = md is the point B. Between this 
point B and the wall there must be the same quantity of 


fluid flow as between Bp and the wall at 0 and, hence, the 
same as between B’ and the wall in Fig. 6b. If 


es (: T : : Mt )(p = Po) 
i 2 
NW(n — 2) = " 5 
16 “ } 7M." po ( ») 


this continuity requirement is approximately met pro- 
vided that (p — o)/p) and X are small so that 1 < 


b/pPo < 1.2, say, and cos (7/2)X ~ 1 — (md?/4). 
Finally, near the wall the pressure gradient must be 
balanced by the viscous stresses—i.e., 


dp/dx = p(0°u/dy’) 


Hence, 


BETWEEN 
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dp 37u10l 0 - a l 
= -{)+ - M,? (6) 
dx n*6," 2 


Eqs. (4) to (6) determine the velocity profile in terms 
of the pressure and the pressure gradient. The pressure 
is related to the displacement thickness of the boundary 
layer. If (m — 2)A < 0.25, say, the displacement thick- 
ness can be shown from Eqs. (4) to be approximately 
equal to 


56* = 6,* + 6. (7) 


where 6,* and 6, are given by Eqs. (2) and (3). Since 
the pressure changes are assumed to be small, the angle 
of deflection of a of the external flow from the free- 
stream direction parallel to the wall is related to the 


pressure by the equation 
a = VM? — 1(p — po)/yMopo 


upstream of the shock. Where there are sharp longitu- 
dinal pressure gradients a is, in fact, somewhat in- 
definite, since the divergence of the stream lines in the 
external flow is of the same order as that within the 
boundary layer; but it is assumed that with a flat wall 


a = (dé*/dx) — (dé,*/dx) 
Hence, upstream of the shock 


d (* — &*) = VM — l(p — po) (8) 
dx y Mo" po : 

In the above analysis the division of the boundary 
layer into an inner viscous flow and outer inviscid flow 
corresponds closely to the exact mathematical analysis 
for small perturbations in reference 5. Moreover Eqs. 
(5)-(8) above can be solved for the logarithmic decre- 
ment with x of a small disturbance, and results nearly 
the same as those of reference 5 are obtained. 


(3.2) The Interrelation Between the Pressure, Pressure 
Gradient, and Velocity Profile 
The physical interpretation presented in Section (2) 
above is partly based on the following facts deduced 


from Eqs. (5) and (6): 


Eqs. (5) and (6) show that for a given m the pressure 
gradient is inversely proportional to the square root of 
(p — po)/po. This suggests that the pressure gradient 
in the region of the reattachment point is at any rate un- 
likely to increase greatly with increase of shock strength, 
although Eqs. (5) and (6) are not really applicable for 
pressure ratios so great as occur in practice at reattach- 
ment. 


Eq. (5) shows that m must be greater than 2. It also 
implies that for a given pressure ”°\ is a minimum when 
n = 8/3. When nis3, so that the velocity gradient with 
respect to y is zero at the wall and the profile is of the 
type which occurs at a separation or reattachment 
point, °X is little different from the minimum. If, 


however, 7 is close to 2, corresponding to a_ well- 
’ Z 
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separated profile, or is large, corresponding to a profile 
not much different from the zero pressure gradient case, 
n°) can be considerably larger than the minimum value. 
(However, the theory would break down if m were as- 
sumed to be extremely large or close to 2 with (p — py) + 
po» not very small.) From Eq. (6) it follows that the 
pressure gradient is near to the maximum for a given 
pressure at the separation or reattachment condition 
but is considerably smaller when is near 2 or is large. 
In physical terms, when the profile is of the separation 
type, as in Fig. Sb, 0°”/Oy? at the wall is bigger than 
when there is a thick dead air region, as in Fig. Sc, or 
when the profile is near to the zero pressure gradient 
form, asin Fig. Sa. It should be noted that the pressure 
gradient can be large when the boundary layer is 
separated provided the dead air region is not too thick; 
in fact, the maximum pressure gradient for a given 
pressure occurs when 2 = 8/3—1.e., when there is a 
certain amount of backflow. 


(3.3) The Pressure at Separation 
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at s, and 
(d/dx)(6* — 6:*) = 0 
at 0. Accordingly, if (6* — 6.*) can be assumed ty 


vary roughly parabolically with x in the interval 0s, the 
length Os, 


(6* i 3.7), =" (6* = a." 6 


lj d 
(6* — 6,*) 
2 Ldx , 


2625As Y Mo" po 
Vv Me? a 


Xe Xe = 


l(ps — po) 


by Eq. (7). Also, since the pressure gradient is small at 
0 and the pressure varies roughly parabolically with x in 


the interval Os, 


1 (dp 
rms (“?) ~- 


Xp) 


The pressure does not begin to change appreciably rur0Uo (1 + = Met) 1M 
. . » 
from fp) until downstream of the point 0. Hence, by Eq. = = 
(8), 3625V My? — 1(ps — po) 
d (3* — 6,*) = V Mo? — 1(s — Po) by Eq. (6) and the fact that » = 3 at the separation 
dx j yMo*po point s. Hence, 
= " ce 
, 0.636 tan 1(? = ) | : 
_ : (: += ue) we , “ 
Ps — Po _ 0.7807 Mo? 2 [(y — 1)/2] 7*Mo 
bo ss [(My? — 1)R,]'“4 1 + 0.693(y — 1) Mi? 
by Eq. (3) and the fact that piol/o/yAo?po = x. Rr. overall characteristics of interactions for shocks so 


The values for (p, — po)/ po calculated by means of this 
relation for the cases y = 1.4, R,, = 0.25 X 108, and AM, 
= 1.5, 2, 3, and 4 are 0.065, 0.095, 0.175, and 0.272, 
respectively. The first three results are not far from 
the experimental values,’ but the figure for M1) = 
probably too high. This may be due to the fact that the 


4 is 


zero pressure gradient velocity profile departs consider- 
ably from the assumed form, Eq. (1), at the higher 
Mach Numbers. It would have been practicable to 
have used the actual zero pressure gradient profile in- 
stead of Eq. (1), but the analysis would have been con- 
siderably more complicated. 


(4) BOUNDARY LAYERS TURBULENT OVER THE WHOLE 
REGION OF INTERACTION 


A method, similar to but much cruder than the 
analysis given in Section (3) for laminar cases, is first 
presented for the pressure at separation in Section (4.1). 
Semiempirical relations are then obtained covering the 


strong that the distance between /’, the peak pressure 
position, and 0 is several times greater than the dis- 
The dis- 


tance x» — X» between F and 0 and the thickness 6, of 


tance between s and 0 (see Figs. 2 and 3b). 


the boundary layer at F are regarded as the primary un- 
known parameters of the interaction to be determined in 
terms of the boundary-layer thickness 4) at 0, the peak 
pressure pr, and the free-stream Mach Number JM. 
The position of 0 relative to the shock does not enter 
into the equations, but, when comparing with ex- 
perimental data on the distance between 0 and J or be- 
tween 0 and /’, it is assumed that 0/ is about (1/2) X 


(xe — X) and OJ’ about (2/3)(x- — X»). 


(4.1) The Separation Pressure 


In incompressible flow and also at a Mach Number of 
2.5, a good fit to the turbulent zero pressure gradient 


velocity profile is given by the relation 











$$ —__ — 


which 


up to, 


This 

u/ Us 
presst 
steep! 
pressi 
fluid 

excep 
presst 
part : 
mate’ 
the a 
acros 
Ap ji 
fluid 
() the 
0.6 p 
than 
thick 
for p 
high 
be bi 
regio 
roug! 
tion 

poin' 


At \ 
the \ 
are C 
4 ca 
This 
chie! 
satis 
lines 
ary | 
at t 
tem) 
relat 
shay 


(4.2) 


It 
Fig. 
y/ 0; 
bou: 
laye 
pon 





ned to 


Os, the 


all at 
h x in 


ition 








SO 





ure 
lis- 
lis- 


of 





in 
ak 
f). 


er 





INTERACTIONS 


u/U = (y/8)” 


which is here assumed to obtain for all Mach Numbers 


up to, say, 4. Thus, at the position 0, 


u U = (y by) 7 


This profile, sketched in Fig. 9, has a “‘shoulder’’ at 
u/Uy ~ 0.6. Experiments show that the bulk of the 
rise to the separation point takes place 


pressure 
so it is reasonable to suppose that in this 


steeply, 
pressure rise region the pressure gradient forces on the 
fluid are considerably greater than the friction forces 
except near the wall. If this is so, the changes that the 
pressure increase causes in the velocities in the outer 
part of the boundary layer can be calculated approxi- 
mately by the application of Bernouilli’s equation, on 
the assumption that the total temperature is constant 
across the boundary layer. Consider a pressure increase 
Ap just sufficient to bring to rest without friction the 
fluid at the shoulder u/ U) = 0.6 of the profile at 0. At 
( there is only a small flow of fluid between the u/l = 
0.6 position and the wall, so for pressure increases less 
than Ap it is perhaps unlikely that there can be any 
thick dead air region near the wall. On the other hand, 
for pressure increases greater than Ap all the relatively 
high velocity fluid at the shoulder of the profile at 0 will 
be brought to rest, and there will be a thick dead air 
region. Hence, it is plausible to assume that Ap is 
roughly equal to the increase of pressure at the separa- 
tion point-—i.e., that the pressure p, at the separation 


point s is given approximately by 


l 
1+ ~— Mt? 
S : (9) 


¥-— 1 

1 + 0.64 —~— M,? 

At Mach Numbers of 1.5, 2, 3, and 4 this ratio assumes 
the values 1.52, 1.84, 2.51, and 3.08, respectively, which 
are close to the experimental values except for the Wo = 
4 case where the experimental value is rather lower.* 
This approximate experimental verification is really the 
chief justification of Eq. (9), which is not based on 
satisfactory theoretical grounds. An analysis on the 
lines of that given in Section (3) above on laminar bound- 
ary layers, taking a more adequate account of the fluid 
at the wall, would be preferable. This was not at- 
tempted, however, because of uncertainties as to the 
relation between the turbulent friction stress and the 


shape of the velocity profile. 


(4.2) The Reattaching Boundary Layer 


It is assumed that at the reattachment point R (see 
Fig. 2), u/U, and u’v’/ Up? are invariant functions of 
¥/drx, where L’p is the velocity at the edge of the 
boundary layer, 5 is the thickness of the boundary 
layer, and wv’ and v’ are the fluctuating velocity com- 
ponents parallel and perpendicular to the wall. 
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This assumption seems plausible for the following 
reasons: 

(1) The zero pressure gradient profile does not vary 
appreciably from the form 


u/U = (9/6) 


at least up to Mach Numbers of 2.5,° so it is reasonable 
to suppose that the profile of uw at reattachment is not 
greatly affected by Mach Number. 

(2) It is shown in the Appendix that the distribution 
of u'v'/U* with y/6 does not vary greatly with Mach 
Number in the zero pressure gradient case except where 
sharp gradients of the mean density occur. This result 
is, moreover, to be expected because where the mean 
density gradient is small something like von Karman’s 
incompressible flow ‘‘similarity hypothesis,” relating 
the turbulent fluctuations to the lower derivatives of 
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the mean u-velocity, should be applicable. Since the re- 
attachment w-profile will probably be roughly linear 
with y, there will not be any sharp gradients of mean 
density at the reattachment point such as there are 
close to the wall in the zero pressure gradient boundary 
layer. Hence, it is reasonable to suppose that at the 
reattachment point, “‘v'///,? will be a function of y/é, 
only independent of Mach Number. 

If there is no heat transfer to the wall and the 
Prandtl Number does not differ greatly from 1, the 
total temperature will not vary much through the 
boundary layer, and the density near the wall at the re- 
attachment point will be given approximately by 


p = pwPr/t1 + [(y — 1)/2]Me?} bo 
where pio is the free-stream density and p, is the pres- 
sure at the wall at R. Now the friction stress 7 is equal 


to —pu'v’, so that near the wall at the reattachment 
point (1/p )(O7r/Oy) will, according to the foregoing as- 
sumption, be proportional to 


Pio RU,” a YM? pRUR* 
(1 + : 9 Me) btby (: + : 9 ue) poU oer 


since pjo/Po = yMo?/ U0’. 

From the boundary-layer equations, it follows that, 
at the reattachment point, (1/9)(O7/Ov) near the wall 
is equal to (1/p))(dp/dx). At R, (1/po)(dp/dx) is 
roughly proportional to 


(be — ps)/Po(Xp — Xo) 


Also pr, Up, and 6, are not very different from p,, U,, 
and 6,, respectively. Hence, 


ky Mop; l Te? 
| 
Me) pol 0°Op 


where k is a constant independent of Mach Number. It 
is found that this relation fits the experimental facts 
0.01. Thus, on semiempirical 


py = 2, 
Po\Xr — X) (: 4 Y = 


approximately if k = 
grounds, 
op ; 0.01 yMo* pel lp? 


ace — 1 10) 
sii (: + 7 Mt (br meee * 


In this equation the right-hand side can be evaluated, 
since p, is given by Eq. (9) and velocities at the edge of 
the boundary layer such as U;,, can be calculated 
roughly in terms of the pressures at the wall? from the 
isentropic relation 

8 1 (p/ po)» y¥— |] 11) 

— = 3 = : ( 

l 07 [(y — by; 2|Mo? 

t At the peak pressure position F, the pressure at the edge of 
the boundary layer will probably not differ greatly from that at 
the wall, although nearer the shock there is probably a consider- 
able variation of pressure across the boundary layer. 
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(4.3) Continuity Consideratians 


In Eq. (10) there are essentially two unknowns, 3, 
and (x» — x), so another relation between them is re- 
quired. This is obtained from the continuity condition, 


flow across FF’ — flow across 00’ = flow across ()' F’ 


(see Fig. 2). 
The flow across FF’ is calculated in terms of 5, on the 
assumption that at F 
u/Up = (y/5,)'/" 
This is a profile somewhat “‘fuller’’ than the probable 
profile at the reattachment point R, where the profile 
will probably be roughly of the form 


u/Upz = (y/b,)'" 


with ” between | and 2 and nearer to 1. For a profile 
of this sort is found to occur at a separation point in in- 
compressible flow The 


profile at / must be somewhat fuller than that at R be- 


see, for example, reference 7. 


cause ultimately, of course, the profile returns to the 
(y, 6) 7 form. However, this ultimate form is probably 
not reached until a distance of the order of 106, down- 
stream of R. 
ceptionally high rate of entrainment of the external 
flow at the edge of the boundary layer, continuity re- 


For otherwise, unless there were an ex- 


quirements would cause the boundary layer to diminish 
considerably in thickness downstream of R, and ex- 
periments show that no such reduction in thickness 


occurs. Hence, it is reasonable to assume that 


u/Up = (y/b,)'”” 


at F, which is only a short distance downstream of R. 
From this assumption and the assumption that at F the 
pressure and total temperature are constant across the 
boundary layer, it follows that 


“bp 
flow across FF’ - f pu dy = Opp 5p 
0 


at F 
where 
2p,U; J A l 
QO; = mae sa -1| A log a _ 2| 
pol PLY 7 1) M,? = i A = 1 
71/2 
1+ 7— ln? (12) 
A= — 
fy l 
=: Me 
o 2 7 





At 0 the velocity profile is assumed to be of the form 


r s \1/7 
Uu U) = (y, 00) ef 


* 


flow across 00’ - f pu dy = QopyoU 60 
0 


at U 
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it is possible to calculate an expression for Q» similar to 
Eq. (12) for Q,. However, the expression involves the 
difference of nearly equal quantities, and it is easier to 
use the relation 


Qo = 1 — (6)*/6o) (13) 


where 6)* is the displacement thickness and values of 
§)*/59 as a function of Jp are given in references 8 and 9 
and reproduced in Table 1. 

Finally. an expression for the rate of entrainment of 
flow across the edge of the boundary layer between 0 
and F is needed. Ina zero pressure gradient boundary 
layer, the edge of the boundary layer is inclined to the 
wall at an angle d6/dx, while the mean-flow stream lines 
at the edge of the boundary layer are inclined at an 
angle dé*/dx to the wall. Hence, there is an angle of 
entrainment; y = angle between the edge of the 
boundary layer and the local mean-flow stream lines = 
(dé/dx) — (dé*/dx) with zero pressure gradient. Now 
references 6 and 10 show that, for the zero pressure 
gradient boundary layer, 6* is approximately given by 


0.0475 (1 + 0.35.M;?)x, 


iat | 0.44 j - 
*ue) R,' 4 
9 


6* = 


(1 + 0.88 ~ 


where x, is the x-distance from the (fictitious) leading 
edge of the ‘“‘equivalent fully turbulent flat plate” and 
R, is the Reynolds Number based on x, and the free- 
stream density and velocity. Hence y can be calculated 
from the data of Table 1, and in Table 2 values of y are 
presented for various Mach Numbers and Reynolds 
Numbers. It can be seen that over a wide range of 
these variables y for the zero pressure gradient boundary 
layer is the region of 0.010. Between the positions 0 
and F, y should be considerably higher than this, be- 
cause the presence of a dedd air region near the wall will 
presumably increase the general scale of the turbulence. 
The case should be intermediate between that of the 
zero pressure gradient boundary layer and the boundary 
of a two-dimensional jet emerging into still air. For this 
latter case in incompressible flow y = 0.083 on the side 
of the jet.'' It is assumed here that the mean value of 
between 0 and F is 0.025. This is not very different 
from the value y = 0.03 assumed in reference 12 for a 
separated turbulent boundary layer. (In reference 12 
the concept of the mixing at the edge of the boundary 
layer is used in a more fundamental way than it is used 
here and is related to the changes in the boundary layer 
though of greater 


profile. However, the method, 


theoretical interest than the present semiempirical 
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TABLE 1 


5o* . . 
for u/ Uo = (y/5y)!/7 


bo 
Mo 0 1 2 3 4 
5o* /5o 0.13 0.16 0.24 0.33 0.42 
TABLE 2 i 


¥ for Zero Pressure Gradient Boundary Layer at Various 
Reynolds Numbers R; 





R:\ 0 1 2 3 4 
106 0.0167 0.0158 0.0148 0.0133 0.0123 
5.108 0.0121 0.0114 0.0107 0.0096 0.0090 
107 0.0106 0.0100 0.0094 0.0084 0.0078 


method is more complicated analytically.) It is further 
assumed here that the mean value of p,U between 0 and 
F is 1/2(p1,U, + pieUy), where the densities p; at the 
edge of the boundary layer are given roughly in terms of 
the pressures at the wall by the isentropic flow relation 


pi/ pw = (p ‘po)'/ (14) 


and velocities U/ are given roughly by Eq. (11). 
Hence, the flow across 0’F’ is approximately equal to 


0.025 [(m,U; 4- pip Up) 2] (Xx, — Yes 


and continuity therefore requires that 


,U, + prU; ? 
Q,5~ — Qodo = 0.025 (° " fi ") (xp — %0) (15) 
=~Pi0o 


(7) RESULTS 


From Eqs. (9)-(15) it is easy to calculate 6, and 
Xp — Xo in terms of dy) (and, hence, by Table 1, 69*), 
Pp/ po, and Mo. The results of several such calculations 
are presented in Table 3. 

It is assumed in comparing with experiment that the 
upstream distance measured from 0 to the point /’ where 
the plane of the shock meets the plate is (2/3)(x- — Xo). 
(This assumption is found to fit closely the experimental 
facts.) Fig. 10 shows the calculated points plotted for 
comparison with the experimental curves of upstream 
distance OJ’ divided by 6)* against p,/po for the values 
1.5, 2, 3, and 4 of Mo.2. It can be seen that the agree- 
ment is good apart from the My) = 4 case. Comparison 
with schlieren photographs of the flow shows that the 
theory predicts 6, for My = 4 to be smaller than it 
really is, x» — Xo is also underestimated, but the ratio 
6,/(xp — Xo) is not far wrong. Hence, the errors seem 
to arise from Eq. (15) and not from Eq. (10). 


TABLE 3 


Calculated Results 


M 1.5 1.5 2.0 2.0 
pr Pp 1.60 l 70 2.00 2.25 
br /d 1.4 Se 4 1.2 1.5 
(Xz Xo) / bo 5 13 4 13 
(XF xo) /6 25 70 18 56 
(2/3)(xr — x0)/6 17 47 12 37 


2.0 3.0 3.0 3.0 3.0 1.0 4.0 4.0 

2.50 3.0 4.5 6.0 7.0 7.0 9.5 12.0 

; = 0.9 ee 2.0 3.4 0.9 1.0 1.0 
25 } 19 42 84 13 17 21 
106 13 58 27 254 30 10) 50 
72 9 39 85 169 0) 27 3 
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APPENDIX.—-THE DISTRIBUTION OF y’v’ IN A ZERO TABLE 4 
PRESSURE GRADIENT TURBULENT BOUNDARY LAYER M 0 l ) 3 } 
n 2 23 33 1.3 5.4 


The equations of continuity, motion, and energy for 
the zero pressure gradient boundary layer are assumed 


to be TABLE 5 
—u'v'/U%(db/dx) at Various u/U in Zero Pressure Gradient 
O( pu) O(pv) (p'v’) Boundary layer 
+ 4 SS a @ (16) = 
Ox Oy Oy Se P ‘i 
_ r u/l 0 l Z » + 
1.0 0.0 0 0 0 0 
pu ou + pv ou -_ 0 (pu'v’) —_ p’v’ ou (17) 0.9 0.0594. 0.0588 0.0571 0.0550 0.0550 
or Ov Ov Ov 0.8 0.0840 0.0849 0.0875 0.0915 0.0970 
. : 0.7 0.0933 0.0966 0.1047 0.1150 0.1269 
u2\ y — 1 0 0.0972 0.1092 0.13880 0.1725 0.207 
p=p/jlt+{1-s= M? (18) 
Tig 2 
— : ‘ TABLE 6 
where quantities with dashes are fluctuating compo- 
. oe 2 2 o - 
nents. It is assumed on empirical grounds® that ~! ~~ . 7 ; 
é 7 (5 )/( ) 8 7.98 8.07 8.30 8.80 
u/l’ = (y/6)‘ (19) ox dx ic 
From Eqs. (16), (18), and (19), it follows that 
v = U(db/dx)(u/U)SE — (p’v'/p) 
where 
a . 
: l + Cl + s°) Vf? 
E=(“)i, <_< ee a = : 
ik + a U? 9 ° =a as’ d 
u 2 0 re ° 
1+ (1 — s?) = VW 


FE varies monotonically between 1/8 at the wall and 5* ’6 at the edge of the boundary layer. Eq. (17) becomes 


(0/Ov)(pu'v’) = (pu?/75)(db/dx)(1 — FE) 











so that 


u\y—-1.. dd ff) s8(1 — E) ds 
uv’ = —-J}14+(1-— M? | U? 
U2} 2 dx Ju/t ae] 


1+ (1 — s*) MW 


The factor 

(1 — E)/f1 + (1 — s?)[(y — 1)/2]M?! 
occurring in the integral assumes the values 7/8{1 + [(y — 1)/2]\/°} at the wall and 1 — (6* 
the boundary layer. It can therefore be replaced by the approximation 


— |] ¥ ; 
i s(1 + 7 wr) | 7 ! ~ ~ - }7/8(1 +7 . ur) [fs 


6) at the edge of 
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where the constant 7 is chosen to give the correct value of 4’v’ at the wall. Here the friction stress 


Got piu’v’ 


where 8 is the momentum thickness. Hence, 


d6 db 6 


in ws 2 
ha oo eee 


1 — (6*/6) — 9(6/5) 


(0/5) — (7/72 


11 + [(y — 1)/2]M?}) 


The ratios 6*/6 and 6/6 are given as functions of .\/ in references 8 and 9, so that m can be calculated. The 


resulting values for Mach Numbers up to 4 are given in Table 4. 


Hence, the quantity r = —u’v’/ L*(db/dx) can 


be calculated for various values of the velocity ratio u/ l’ and Mach Number /, and results are presented in Table 


od. 


It can be seen that in the outer part of the boundary layer r does not vary much with .\/, though it varies a great 
deal near the wall. The gradient d5/dx of boundary-layer total thickness also varies little with Mach Number; 
the gradient dé* ‘dx of displacement thickness increases, but (dé*/dx)/(d6/dx) also increases. According to ref- 


erences 6 and 10, 


0.0475 (1 + 0.35.M*)x, 


5* = 


{1 + O.88[(y — 1)/2]M2}°-44R 1/9 


where x, is the .-distance from the (fictitious) leading edge of the ‘‘equivalent fully turbulent flat plate’ and R, is 
the Reynolds Number based on x, and the free-stream density and velocity. Hence, 


d6* d6,* 7 
dx dx 


where 6,* is the displacement thickness in incompressible flow at the same Reynolds Number &,. 


i 


1 + 0.35M? 


{1 + O.SS[(y — 1)/2]M?} o44 


By the use of 


Table 1, (dé dx) /(dé,;*/dx) can be calculated, and the values are given in Table 6. It can be seen that there is 


little change up to Mach Number of 4. 


Hence, the rough constancy of r in the outer part of the boundary layer 


means that the distribution of 4’y’, LU? with y 5 for zero pressure gradient boundary layers with the same Reynolds 
Number R, does not vary greatly with Mach Number except near the wall, where there are large gradients in the 


mean density at the higher Mach Numbers. 
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Some Basic Concepts for Analyzing 
Dynamic Flight-Test Data 
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ABSTRACT 

The problem considered here is that of finding, through the 
use of observed data, the coefficients of the differential equations 
that describe a dynamic system. This problem occurs in the de- 
termination of the aerodynamic stability derivatives of an air- 
plane or missile from flight-test data. There is a maximum num- 
ber of coefficients which may be determined. A method is given 
for finding this number, and conditions are given under which 
the maximum number is obtained. These results are summa- 
rized in a working rule and are applied to two examples. It is 
shown that the free oscillation type of flight testing can never 
yield all the coefficients. Suitable forcing functions are de 
scribed which do permit all coefficients to be evaluated. 


INTRODUCTION 


_— USUAL PROBLEM IN DYNAMICS is the following: 
Given the forces acting on a system, find the re- 
sulting motion of the system. The solution to this 
problem is obtained by setting up and solving a set of 
differential equations subject to given initial condi- 
tions. This will be called the direct problem. The in- 
verse problem would then consist of finding the in- 
ternal forces acting within a system given the motion 
of the system. The solution to the inverse problem 
would require, in effect, constructing the equations of 
motion, given the solution of these equations. The 
type of problem considered in this paper may be called 
a semi-inverse problem, since the form of the equations 
of motion is assumed given. This means that in every 
equation each term involving the generalized coordi- 
nates, their derivatives, and the forcing functions is 
given except for an unknown multiplicative constant. 
The solution to the semi-inverse problem consists of 
finding these unknown constants, the constant co- 
efficients in the equations of motion. 

The importance of the semi-inverse problem in aero- 
nautics lies in its application to dynamic flight testing. 
These tests are normally accomplished by subjecting 
the aircraft to a disturbance (usually a control deflec- 
tion) and measuring the resulting motion. The object 
of such testing is the determination of the aerodynamic 
forces and moments, or their effects, which are con- 
ventionally represented by stability derivatives, transfer 
functions, or indicial responses. The stability deriva- 
tives are the coefficients of the differential equations 

Presented at the Aircraft Design Session, Twenty-First Annual 
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describing the motion of the aircraft. The transfer 
functions are the responses to a sinusoidal input. The 
indicial responses are the responses to a step function 
input. The convenience of each of these descriptions 
of a dynamic system depends upon what use is to be 
made of the results. Table 1 gives a comparison of 
certain features of these descriptions from the point of 
view of the problem of flight testing for aerodynamic 
data. 

Because of the desirability to an aerodynamicist of 
knowing the coefficients of the differential equations, 
the work described in this paper was undertaken. 

This paper considers the theoretical aspects of the 
flight-testing problem exclusively from the point of 
view of obtaining the coefficients of the differential 
equations. These coefficients are found by solving 
sets of linear algebraic equations. It is shown that it 
may be impossible to find all coefficients and that there 
exists an upper limit to the number that can be found. 
A method for finding this number is given, and three 
conditions are stated which must be satisfied in order 
to find this maximum number of coefficients. The re- 
sults are summarized in a working rule whose applica- 
tion is demonstrated by two examples. A method is 
pointed out for increasing the number of coefficients 
which may be determined. 

For the theoretical considerations in this paper, it is 
assumed throughout that the equations of motion are 


TABLE | 
Comparison of Descriptions of Dynamic Systems 


' Type of Description 

Consideration Transfer functions or Coefficients of 
indicial responses differential equations 
Test Time history must be Time history may be 
conditions taken in a period in taken in any period of 
which steady-state time. Steady-state 
conditions exist at conditions need not 
beginning and end exist at begimning or 

of period end of period 

Can be applied only to Can be used in linear 


Application 
linear systems and nonlinear systems 


to nonlinear 

problems 
Usefulness of Results are useful for Results are useful for 
synthesis of control checking aerodynamic 
systems. Responses design calculations 
can only be com- Responses for any 
puted for some in kind of input (gust, 
put variables—e.g., control, etc.) may be 
gustresponsescannot computed by well- 
be computed from known methods 
measured responses 
to controls 


results 
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linear with constant coeflicients. Equations of this 
type result from the linearization of the exact equa- 
However, in certain cases lineariza- 
In those cases the method 


tions of motion. 
tion may not be justified. 
of analysis discussed here can be applied to the follow- 
ing more general types of equation of motion. If the 
equation consists wholly of terms that are composed of 
an unknown coeflicient times a known function of time, 
these coefficients may be obtained as the solution to the 
semi-inverse problem. This known function of time 
may appear naturally in the problem, as for example, 
the dynamic pressure; or it may be a nonlinear func- 
tion of the generalized coordinates that would occur if 
the linearization process is not valid. To determine 
the maximum number of coefficients that may be found 
in the solution to the semi-inverse problem for cases in 
which the more general types of equations occur, the 
same method of analysis given here is applicable. How- 
ever, two of the three conditions that are found to en- 
sure a unique solution cannot be applied to a system 
that is described by nonlinear equations. Throughout 
the discussion it is also assumed that the data giving the 
time history of the generalized coordinates, their deriva- 
tives, and the forcing functions are free from experi- 
For any practical experiment this will 

The problem of developing a suitable 
other questions related to 


mental error. 
not be the case. 
averaging process and 
handling inaccurate data are being studied and will 
be considered in a later paper. 

The idea of determining coefficients in the equations 
of motion from experimental time histories is not new. 
A limitation on the number of coefficients which may be 
determined was pointed out by Greenberg! for a two 
degree of freedom system. The present paper considers 
systems with an arbitrary number of degrees of free- 
A summary of the status of this problem has 
been given by Milliken.? Previous work on this prob- 
lem is separated into two categories. The first cate- 
gory, ‘““Equations of Motion,” is essentially what is 
called here the semi-inverse problem. The various 
methods of approach listed differ only in the details of 
the necessary calculations. The category, 
“Response Curve Fitting,’ is a collection of methods 
using various devices to extract from the data certain 
combinations of the differential equation coefficients. 
These combinations are called transfer coefficients in 
reference 1. These methods are usually applied to the 
case of a system with two degrees of freedom. 

Attempts have been made to generalize the response 
curve fitting techniques to a system with more degrees 


dom. 


second 


of freedom by using the concept of transfer functions. * 
There are inherent difficulties in obtaining the transfer 
coefficients for a system with more than two degrees of 
freedom. First, there is a practical difficulty since a set 
of nonlinear algebraic equations must be solved. 
Second, there is a theoretical difficulty, since the trans- 
fer coefficients may satisfy one or more functional rela- 
uons that can apparently be discovered only by trial- 
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Obtaining the differential equa- 
irom 


and-error methods. 
tion coeflicients—1.e., the stability derivatives 
the transfer coefficients is also difficult, because in the 
general case another set of nonlinear algebraic equations 
must be solved. Furthermore, the transfer function 
concept can only be used for a system described by 
linear, constant coefficient differential equations. Evi- 
dently the transfer function concept, which is useful for 
other purposes, is not well adapted to the solution of the 
semi-inverse problem with several degrees of freedom. 
With regard to the notation employed, it was not 
practical to use conventional aerodynamic notation 
when discussing the general problem. The coefficients, 
a;j;, can represent stability derivatives and the inertia 
properties of an airplane or missile.. The correspond- 
ence will be apparent by referring to the section on ex- 
amples. The arrangement of the paper is such that, if 
one is interested only in applying the results to specific 
examples, the section on “Practical Application and 


Examples” can be read at once. 
THEORY 


One Degree of Freedom Case 

The semi-inverse problem for a simple example will 
now be discussed in detail. This will illustrate the prin- 
ciples involved in the general problem and the methods 
to be used, while avoiding extensive calculations. 
Problems of greater generality will be considered in the 
next section. 

A linear spring-mass-damper system executing free 
vibrations is the example to be discussed. Data giving 
the displacement, velocity, and acceleration of the mass 
as functions of time are measured directly or computed 
from measured data by integration or differentiation. 
This system is assumed to be described by a second- 
order, linear, constant coefficient differential equation. 
It is assumed that all three terms are present in this 
equation—i.e., the coefficients of the mass, damping, 
and spring terms are not zero. If all three terms are not 
present in the equation that actually describes the sys- 
tem, this will be discovered in the calculations. 

If x denotes the displacement of the mass, the equa- 


tion of motion assumed is 


AX + dox + a3x = 0 


and the problem is to determine a), d2, and a3. Not all 
three of these constants can be found because of the 
homogeneity of the equation. The trivial solution, 
Q, = a2 = a3 = 0 


isexcluded. But it would seem possib!e to solve for two 
of the coefficients in terms of the third or equivalently 
solve for the ratios of two of the coefficients to the third 
For example, suppose it is deciced to solve for 
The above equation would 


one. 
the ratios a;/d. and da3/ds. 


be written 
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(a;/a2)* + (a3 Qa)x = —x 


Select two values of time, /; and f,, and substitute the 
corresponding known values of x, x, ¥ in this equation, 
giving 

(a,/d2)¥(t)) + (a3 d2)X(ty) = —%x(t) 


(dy, 2) (to) + (a3 a2) x(te) = —X(te) 


If the determinant 


is not zero, then it is possible to solve for the ratios 
uniquely. 

The condition of a nonzero determinant is necessary 
and sufficient for a unique solution, but there is no as- 
surance that it is possible to choose ¢; and f, such that 
X| #0. There are three situations that will cause |X 
to vanish. First of all, an unfortunate choice of the two 
values of time will give a zero determinant; for example, 
if ¢; and f) are such that x(t;) = x(f2) = 0. But this can 
be circumvented by choosing other values of time. 

Second, as mentioned above, the assumed form of the 
equation of motion may be wrong. Suppose, in that 
case, that a; and az are solved for in terms of a2, but the 
physical system is actually one without damping 
i.e., dg = 0. Then [x = 0 for any 4; and f. Again 
this difficulty is not serious, since it is known then that 
ad, = 0 and the ratio a,/a3; can be found. The data thus 
still give two of the coefficients of the system. On the 
other hand, suppose that the assumption of nonzero 
mass, @, is false. Then the equation of motion would 
be of first order, and it would be possible to solve for no 
This is so because the solution would 





two coefficients. 
be 
NM 
x = ce 
where ¢ is an arbitrary constant and \ is the ratio, 
—(d3/d2), and the three possible second-order deter- 
minants from the matrix 


x(t,) 
x(to) 


¥(t) 
¥ (te) 


X(te) 
would all be zero for any ¢; and f,. In this case it is 
possible to solve for only one coefficient in terms of two 
others, one of which is zero. The difference between 
the cases a; = O and a. = O is that the first reduces the 
order of the equation while the second does not. The 
case a; = 0 is more crucial, since all determinants of the 
above matrix are zero. This second situation shows 
that the assumed form of the equation of motion deter- 
mines what coefficients may be solved for and whether 
the expected maximum number can be solved for. 
An equivalent way of expressing this latter conclusion, 
which is more convenient for multiple degree of freedom 
systems, is: If the order of the true characteristic equa- 
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tion describing the system is greater than, or equal to, 
the order of the characteristic equation describing the 
assumed system, then the maximum possible number of 
ratios can be obtained from the data provided certain 
other conditions are met. 

There is a third situation that will make it impossible 
to solve uniquely for two ratios. It is known that the 
general solution of the second-order differential equa- 
tion is 


nat ~ prel 
x = (3é + C2 


and that the initial conditions determine c; and ¢. 
The roots of the characteristic equation, \, and \2, may 
be real, imaginary, or complex. Consider first that 
Further, suppose the initial condi- 
Then the solution is simply 


A; and Xz» are real. 
tions are such that c. = 0. 


x = ce 

Again, all of the possible second-order determinants of 
the coefficients arising in solving for two ratios will be 
zero. Note that the above solution for x satisfies any 
number of distinct second-order differential equations, 
all with the same initial conditions so that the semi- 
inverse problem would not have a unique solution in 
this case. If \, and \, are pure imaginary, this implies 
that a2 = 0, which has already been considered. If 
A; and 2 are complex with nonzero real part, then all 
of the three possible second-order determinants are 
The conclusions concerning this third situa- 
To be able to ob- 


nonzero. 
tion can be summarized as follows: 
tain two (the maximum possible number) of the co- 
efficients of the system in terms of a third, the initial 
conditions must be such that all modes are excited. 

In summary, the results for this example show that 
to be able to solve for a maximum of two coefficients of 
the given system from the data alone, the following three 
conditions are necessary and sufficient. 

(1) A fortunate choice of values of ¢ is made. 

(2) The order of the true characteristic equation is 
greater than or equal to the order of the assumed char- 
acteristic equation.* 

(3) The initial conditions are such that all modes are 
excited in the motion. 

These three conditions are of importance to any at- 
tempted solution of the semi-inverse problem. How- 
ever, it should be noted that-it cannot be determined 
whether or not they are satisfied except by evaluating 
the determinants. In practice the first condition will 
not be troublesome, since for actual data some averaging 
This, in effect will nullify an 
It is possible 


process should be used. 
accidental choice of unfortunate times. 
that the second condition might be violated in the 
study of airplane or missile dynamics. However, there 
will usually be sufficient knowledge of the physical 

* The “‘correct’’ values of the coefficients can be found only if 


the orders of the assumed and true characteristic equations are 
equal 
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system so that it would be known how to change the 
The third 
determine 


assumed form of the equations of motion. 
condition presents a dilemma, since, to 
whether all the modes are excited or even what the 
modes are, the coefficients of the system must be known. 
However, from all possible types of initial conditions, 
only a relatively small number will force one or more 
of the arbitrary constants to vanish. Thus the prob- 
ability that condition three is violated is extremely 
small. 

For this simple example, two was the maximum num- 
ber of coefficients of the differential equation that 
could be obtained from the data. To be able to discuss 
more complicated examples, a systematic way of finding 
the maximum number of coefficients is needed. The 
maximum number of coefficients that can be solved for 
is determined by the rank of a matrix. The kth row 
of this matrix consists of the quantities x, x, and * evalu- 
ated atatime?/,. The rank of this matrix is the number 
of coefficients that may be solved for in terms of other 
coefficients. Evidently the rank must be less than 
three. The maximum possible rank is two, though 
from the previous discussion it may be less than two. 

A systematic way of deducing the maximum possible 
rank of a matrix is offered by the geometrical concept 
This concept will be introduced by apply- 
Let the three coefficients 


of vectors. 
ing it to the above example. 
of the differential equation, a), a), and a3, be the com- 
ponents of a vector A, and let *, x, x be the components 
of a vector X = X(t). Let 


Ay = X(t,) 


In this notation the equation of motion is then 


Multiple Degrees of Freedom Case: 


Consider first the free vibrations of an m degree of freedom physical system. 
For an “‘ordinary”’ system the equations of motion will be 


be x), ...Xn- 


Q1X) + ae + GinXe + Qa, n+1 xy + oe + QA, 2n =. + Q1, 2n+1X1 +... 


Guts Po ss TP Geka 


The term ‘‘ordinary”’ is used to emphasize the fact that 
it is possible to have higher order derivatives in the 
equations of motion of certain physical systems. This 
is true, for example, if downwash lag effects are com- 
pletely considered. However, the presence of such 
terms introduces no essential changes in the discussion. 
Define the following sets of numbers as vectors in a 
vector space of 3” dimensions: 


A;: (dit, Qj2,.- 


_ 
bo 


+» diy 3n) 
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A-X = 0 


Thus the fixed vector A is perpendicular to the X vector 
at all times so that all vectors X, lie in the plane per- 
pendicular to A. The end point of X traces out a 
curve in this fixed plane, as shown in Fig. 1. 

The concept of the linear independence of vectors 
will be used in the following discussion. The m vectors 
X; for j = , n are linearly independent if it is 
impossible to find m scalars, C1, C2, . . . C,, Where at least 
one is not equal to zero, such that 


) 
ie: 


1X, + OX. +... + CnX, = 0 


It can be shown that, for vectors lying in a plane, at 
most two are linearly independent. For vectors in 
three-dimensional space, at most three are linearly 
independent, so that any vector may be expressed as a 
linear combination of three linearly independent vec- 
tors. 

The rank of the matrix defined above is equal to the 
number of linearly independent rows—4.e., the maxi- 
mum number of lineariy independent X, 
Since all X, lie in a plane, there are at most two linearly 
independent X, vectors. 

The same result can be obtained analytically as fol 


vectors.‘ 


lows: In a three-dimensional space there are at most 
three linearly independent vectors. Since A must be 
one of these (any two vectors that are perpendicular 
are linearly independent) there can be at most two 
linearly independent X, vectors. Thus the maximum 
possible rank of the matrix is two. 

The above methods will now be applied to more gen- 


eral physical systems. 


Let the generalized coordinates 


+ Q, 3axXn = O 
(1) 


a Qn, 3nXn - 0 ) 


a Wins: I (3) 


oe: ee 


Then the equations of motion, Eqs. 


and Z, = X(t;). 


(1), can be written as 


A,X =0 (2a eee 
There are » fixed (but unknown) vectors A; The de- 
termination of the components of A; is the semi-inverse 


problem. These » vectors are linearly independent 
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because the equations of motion are linearly inde- 
pendent. Therefore, in the 3n-space, there are at 
most 2m linearly independent vectors X,. Call the 
matrix, which has for its kth row the components of X;, 
the X matrix. The rank of this matrix is at most 2n. 
This means that from each equation of motion at most 
2n of the a;; can be solved for in terms of the remaining 
n. This is a generalization of the conclusion reached 
in studying the one degree of freedom case. 

In order to be able to solve for this maximum num- 
ber, the previously stated conditions 1, 2, and 3 must 
be satisfied. That these conditions are necessary and 
sufficient in order to obtain the maximum number of aj; 
from the data will be demonstrated by examples. It 
is obvious that condition 1 must be satisfied, but, since 
some averaging process should be used with real data, in 
practice it is of no great concern. To illustrate why 
conditions 2 and 3 must be satisfied for multiple degrees 
of freedom, assume that the physical system is repre- 


sented by Eq. (1) form = 2. 
AnX, + Ay. + Aish + AuX. + A15X1 + Aix, = O 
AX, + Aro. + Ao3X, + Aske + Ao5Xi + AoeX2 = O 


According to the rule for determining the maximum 
number of coefficients that can be found, it is possible 
to solve for at most four of these in terms of the re- 
maining two in each equation. But suppose that for 
the actual physical system 


lau Aye 
| |= 0 (4) 
|Qer 22 


and that the characteristic equation is of degree three. 
Then x; and x2 are linear combinations of three expo- 
nentials. It can be shown that the X matrix for this 
example cannot have rank four. In fact it is at most 
of rank three, and thus only three of the a,; can be 
solved for in terms of the remaining three in each equa- 
tion. Similar considerations hold for the n degree of 
freedom case. 

To show why condition 3 must be satisfied, some 
further calculations are necessary. It is no longer as- 
sumed that Eq. (4) is satisfied. Assuming that the 
roots of the characteristic equation, d,, for p = 1, 2, 
3, 4 are distinct, it can be shown that 


~ -Apl 
x ce 


pt 


. 
Xe CpK ye 


1 
p=1 
4 
CF 
where the c, are arbitrary constants and the A, are 
constants involving the a,;;. Suppose it is desired to 
solve for a1, and a»; for 7 = 3, 4,5, 6. Pick four values 
of time ¢, and denote the matrix, whose kth row is 


X(t; \xto(t,) x4 (ty) Xo(ty ) 


by Y. Define the matrix 
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A 


a 








Fie. 1. 


e = (e*) (p, k = 1, 2, 3, 4) 


where k denotes the row and p the column, and the 
matrix 


A Kw 1 Ky 
as sae Az Korg l Ky 
= wes Ks 1 K; 
NM Kina 1 Ky 
Then it can be shown that 
a 0 O O 
0 ag 0 0 
7". 6 a o1* 
0 0 O & 


To solve for the desired a;;, Y must be nonsingular 
If condition 1 is satisfied, e is nonsingu- 
Then 


i.e., of rank 4. 
lar, and it can be shown that A is nonsingular. 
the rank of Y is equal to the number of nonzero ¢,. 
The rank of Y is four if all modes are present in the 
solution—i.e., if allc, # 0. These results can be gen- 
eralized to any number of degrees of freedom. 


Returning to the multiple degree of freedom problem, 
it must be emphasized that, when it is decided that the 
rank of X is 2n, this means that there is at least one de- 
terminant of order 2” in X which is not zero. There- 
fore it cannot be said that any set of 2” of the a,; can be 
solved for in each equation. It may happen that only 
certain sets of 2” coefficients may be solved for in terms 
of the remaining 7. 


In determining the maximum number of a, that can 
be obtained from the data, the number of degrees of 
freedom and the dimension of the suitable vector space 

The former would be known in any 
For the general » degree of freedom 
However 


must be known. 
given problem. 
system the latter has been taken to be 3n. 
this requires some discussion. Suppose that a; 3, = 0 
for all 7 so that x, does not appear in any equation. 
Should the appropriate vector space be one of dimen- 


sion 3 or 3n-1? Either of these can be chosen if a 
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certain precaution is taken. If a 3n-1 space is chosen, 
then by following the same reasoning as above it would 
be deduced that, at most, 2”-1 of the a;; can be solved 
for in terms of the remaining ” from each equation of 
motion. If a 3n-space is decided upon (by keeping in 
the equations a literal term for a;, 3,%,), it is found that, 
at most, 2” of the a;; can te solved for in terms of the 
However, 
Other- 

Thus 


remaining ” from each equation of motion. 
included in the set of 2” of the a;; must be aj. 3p. 
wise the coefficient matrix would be singular. 
there is no contradiction in the arbitrary choice of di- 
The same amount of informa- 
2n-1 of the nonzero 
It is not 


mension of the space. 
tion is obtained in both cases 
coefficients and the result that a; 3, = 0. 
difficult to prove the statement that a; 3, must be in- 
cluded in the unknowns solved for when a 3n-space is 
chosen. However, only an illustrative example will 
be given here. 

Consider a two degree of freedom system with no 


damping. 


AX + Q12X2 + 143X) 0 


1X) + QoX2 + 3X, = 0 
where x2 does not appear in either equation. A solution 
to the equations is 


x1 = c3 sin wt + c, cos wt 


xo = b, + dot + b3 sin wt + 4 cos wt 


7 


If the appropriate vector space is taken to be a 3-space, 
only one of the a;; can be found in terms of the other 
two from each equation. But suppose a 4-space is 
decided upon (by carrying along a term @4X2 in each 
equation), with the result that two of the a;; can be 
solved for in terms of the remaining two from each 
equation of motion. This is possible if at least one sec- 
ond-order determinant from the matrix 


x1(t,) Xo(t1) 
x1 (te) Xo(te) 


lay (t1) 
¥1 (ty) 
isnonzero. It can be shown that the following relation 
holds between the constants appearing in the expres- 


¥o(t1) 
Xo(te) 


sions for x; and x». 
bz t= bs C4 


Using this fact and the expressions for x; and x2 given 
above, it can be shown that the only nonzero deter- 
the matrix are those that include the 
This means that, to 


minants from 
last column from the matrix. 
solve for two a;; in each equation of motion, a; must 
be one of these. It will be found that aj = 0. 

In discussing practical applications of the foregoing 
theory, which will be done in the next section, it hap- 
pens that some of the a;; are zero or accurately known 
in advance. 
considered, but all of these have been reduced to a 
simple working rule that is stated in the next section. 
In essence this rule says the following: If conditions 


There is a large variety of cases to be 


FLIGHT-TEST DATA 745 
1, 2, and 3 are satisfied so that the maximum number 
of a;; may be obtained, contained in this maximum 
number may be some a;; that are already known (zero 
or nonzero). This occurs regardless of the chosen di- 
mension of the space. 

Next consider the case of an n» degree of freedom 
system undergoing forced vibrations. The system’ of 
Eqs. (1) will be modified only by the presence of terms 
on the right-hand sides. Suppose the forcing function 
for the 7th equation is 6,6, where 6 is given experiment- 
ally and the 0; are either known or unknown. Un- 
known 3b; arise in the case of a control deflection in air- 
plane dynamics where the control effectiveness may be 
unknown. To reduce this case to the type considered 
previously first redefine A; and X as 


» Diy 3ny b;) 


A;: (da, oes 
a: (X1, ese 


Xn, —6) 


Now the appropriate space is of dimension 3n + 1. 
There are nA vectors and thus at most 3n + 1 — ” = 
2n + 1 linearly independent X vectors can be found. 
The number of independent X vectors was shown to 
be the number of coefficients that may be solved for 
in each equation. Therefore from each equation at 
most 2” + 1 of the a;; and 5; can be solved for in terms 
of the remaining ». If all the ); were known, a total 
of n more of the unknown a;; could be found than in the 
case of free vibrations. 

The next question to be considered is: If all the aj, 
are unknown, can an experiment be devised so that it is 
possible to determine all the a,;? The answer to this 
question is an affirmative one, at least in theory. 

By having forcing functions of the form 0,6, the di- 
mension of the appropriate space was increased by one, 
and it was possible to solve for a total of m more of the 
Therefore, further ways of increasing the dimen- 
‘his must be 


Qij. 
sion of the vector space will be sought. 
done with some care. For example, it is not possible to 
write the forcing functions 6/6 as (1/2)b6 + (1/2)d6 
and conclude that the appropriate vector space is of 
dimension 3 + 2. 

To make certain that the dimension of the vector 
space is increased, the concept of the direct sum of vec- 
The definition and a discussion 
However, it is 


tor spaces will be used. 
of this concept is given in reference 5. 
possible to gain an intuitive feeling for this concept 
which is sufficient for the purpose of this paper. Let 
Y represent a vector space of dimension ry. Then Y is 
completely specified if r linearly independent vectors 
from Y are given, since any other vector in Y is a linear 
These r vectors are 
, Y, where 


combination of the given vectors. 
Let the vectors be Yj, ... 


: f = (Vu, oe 


said to span Y. 


» Vir) 


d 


in component form. It is convenient to represent this 


set of vectors as a matrix 
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Suppose another vector space, Z, of dimension s is 


given—i.e., s linearly independent vectors Z,, .. . Z, are 
given. Let the vectors in component form be repre- 
sented by the matrix 
2 TT. 
201 eee a 








i Zs1 


Now define the matrix 





[ ys Vu Se t3298 
2) 2 Oo iS 0 
ane: @ Ba ..+ Be 

ma sox @ Ba ss Be 





It can be proved that these 7 + s vectors with r + s 
components are linearly independent. Therefore they 
span a vector space of dimension r + s. This vector 
space is called the direct sum of Y and Z, which is rep- 


resented by the symbol 
Y@Z 

The result that will be used in the following discussion 
is. The direct sum of two vector spaces is a vector space 
whose dimension is the sum of the dimensions of the two 
spaces. 

The following result will be proved. 
independent functions of f, say 6, .. . 
independent vectors B; 


B;: (by, ee 


Given v linearly 
, 6,, and n linearly 


» Din) 


where the },;; are known. If the equations of motion of 


a physical system are 


where 
F; = bd) + eae +r b indy 


with A; and X defined as in Eqs. (2) and (3), then the 
maximum number of a;; which can be determined is 
3n from each equation of motion. Note that the com- 
ponents of X as defined in Eq. (3) are now obtained from 
This maximum 


the solutions to the system Eq. (5). 


1953 


SCIENCES NOVEMBER, 
number can be determined if conditions 1, 2, and 3 are 
satisfied. 

Denote the original 3n-space by A and the n-space 
spanned by B; by B. The vector space defined by the 
direct sum of A and B 


S=AQ@B 


is a 4n dimensional vector space. The space S can be 


spanned by the » B; vectors, which are now written as 


we GOO... ode) 


with the first 32 components zero, the m A; vectors, and 
2n X,. vectors 


A;: (it, Ae ay 3ns 0, ea ® 0) 


Be: (Wille)... «+ Sele) @... 0] 


with the last 7 components zero. 
However, a different set of 4m linearly independent 
vectors that span the space is desired. Define in S the 


vectors 


A;: [0, eee y 0, 6i(t,), e+ 6% 5p (tx) | 


The A, lie wholly 
Therefore, each 
It can be 


with the first 32 components zero. 

within B, which is spanned by the B;. 

A, must be a linear combination of the B,. 

shown that this is possible if the set of functions 6, 
. , 6, form a linearly independent set. Let 


(2 eee 


, 
A, _ A; + B, 
These vectors are linearly independent. Then it is 
possible to choose a maximum of 3 other linearly inde- 


pendent vectors to span S. For these, take 


X,’ = X, — Ay, Wwe i,... 


, on) 
The A,’ are linearly independent of the X,’, since 
A,’-X,’ =0 


Since 3” is clearly the maximum possible number of 
linearly independent vectors X,’, the maximum pos- 
sible number of a;; which can be determined is 3 from 
each equation. 

Therefore it has been shown that to completely de- 
termine the physical system, » essentially distinct fore- 
ing functions, /; are needed. Since this is in effect 
what was proved, it might be argued that there was no 
point in introducing the function 6;. This was done be- 
cause in the applications to airplane dynamics it is often 
impossible to separate the forcing functions into specific 
degrees of freedom. For example, an elevator deflection 
or gust will produce both a lift and a moment. 

A possible way of introducing additional forcing 
functions is to attach rockets to the airplane or missile 
Then the 6; may be in the nature of 
If the idealized 


at various points. 
short pulses and idealized to impulses. 
impulses are defined to be linearly independent when 
they occur at different times, then the above results 
still hold. 
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APPLICATIONS AND EXAMPLES 


Working Rule 


The result of the foregoing theory may be summa- 
rzed for convenience into a working rule. This rule, 
which presupposes certain stated requirements for the 
data used, allows one to predict which coefficients may 
be solved for, in principle. In practice, inaccuracies 
in the data may cause large inaccuracies in the results, 
but this problem is not considered here. The following 
previously stated conditions are assumed satisfied : 

(1) A fortunate choice of values of ¢ is made. 

(2) The order of the true characteristic equation is 
greater than, or equal to, the order of the assumed 
characteristic equation. 

(3) All modes are excited in the motion. 

The rule states that in each equation the number of 
unknown coefficients which can be determined plus the 
number of known coefficients which may be checked 
are equal to the total number of measured quan- 
tities less the number of equations. This is restated 
below in terms of symbols. The following quantities 
are defined : 


number of equations describing system 
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c¢; = number of coefficients (zero or nonzero) in 
ith equation known in advance and which 
may be checked 

u; = number of coefficients in ith equation that 


remain unknown 


The rule for finding d,, c;,, and u, may be stated as 
If p; Z (m — n) and m; > p,, then d; = p,. 
If 


follows: 
If p; Z (m — n) and m,; = p,, then d; = p; — 1. 


pi > (m — n), thend; = (m — n).  Inany event, 
Cc; = (m — n) — d; 
u; = (pi — dj) 


Solving for Coefficients 

Devising a method of solving for the coefficients, after 
having decided (by the foregoing working rule) which 
coefficients are to be solved for, is a separate problem. 
In addition, the preceding discussion has not considered 
the effect of inaccuracy in the data. In order to mini- 
mize the effect of such inaccuracies upon the calculated 
coefficients, it is evidently desirable to adopt some pro- 
cedure that averages over a great number of times. 
The selection of the “‘best’’ averaging procedure is a 


n = 
m = number of nonzero column vectors in coefli- formidable task, which is not examined here. 
cient matrix 
m, = number of nonzero coefficients in ith equation : aed , 
: . pos ‘ Examples, Airplane Longitudinal Motion 
pf; = number of coefficients not known in advance 
in 7th equation In order to illustrate the application of the working 
d; = number of coefficients in 7th equation that _ rule, two cases of the airplane longitudinal equations of 
may be determined motion are considered. The matrix of coefficients of 
TABLE 2 
(m n) = (4 2)=2 Non- Not Deter- 
zero Given mined Checked Unknown 
a 4 a 4 a 6 m Pp d, C; uy; 
C 
0 0 I — 0 3 1 1 l 0 
. — 9 9 
0 —h Coz Coie C 0 $ 4 2 0 2 
TABLE 3 
Non- Not Deter- 
(m —n) = (5 —2) =3 zero Given mined Checked Unknown 
a 6 a 6 a 6 be mM; pi d; C; uy 
cz ‘ 
0 0 1 1 0 0 3 ] 1 Ss 0 
. : 5 : 2 ; a 
0 —h Cd Cap Cme 0 Cg 5 5 3 0 2 
TABLE 4 
(m — n) = (10 — 4) = 6 Non- Not Deter 
zero Given mined Checked Unknown 
6 b- 7 a 0 5 u a 6 be m; pi d, C; uy 
: CDa—CL Cl . . ' 
0 0 ] 0 0 0 CD zs - 0 } 3 3 3 0 
. CI 
0 0 0 1 = 0 cz = 0 0 } 2 2 j 0 
—h 0 0 Coe Cnae Cmas. Cm, Cme O Crys. 7 7 6 0 I 
l (1) + he) O —h, (Crago + hy) Chas. Ch, Ci, 0 Chg, Ss Ss 6 0 2 


All symbols are take 


n from reference 6 








748 JOURNAL OF THE 
the equations of motion are taken from reference 6 in 
both cases. 

(1) Short-Period Motion Only. 
motion for this case are given nondimensionally in ref- 
The degrees of freedom are 


The equations of 


erence 6 on page 403. 
plunging, a, and pitching, @. 


a — 6+ (C,,/2)a = 0 
—h6 + Cra + Crna + Cra = 0 


where all symbols are defined in reference 6. Rewriting 
the coefficient matrix and applying the working rule 
gives Table 2. 

For the example in Table 2, it may be noted that in 
each equation d; coefficients may be determined in terms 
of u; other unknown coefficients. For example, Cz, 
may be found outright, whereas two of h, Crga, Cae, and 
Cnq may be found in terms of the remaining two. No- 
tice also that adding a moment due to elevator deflec- 
tion would alter the situation, as shown in Table 3. 
Thus, the addition of a forcing function permits check- 
ing one more coefficient from the first equation and 
finding one more in the second. 

(2) Longitudinal Motion Stick-Free. 
for this fairly general case of longitudinal motion with 


The equations 


four degrees of freedom are given in reference 6 on page 
2Q9 
OTe. 


city, wu; plunging, a; pitching, #; and elevator position, 6,. 


The degrees of freedom considered are axial velo- 


The coefficient matrix from these equations is given in 
Table 4. 
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From Table 4 it is clear that in the first row, the three 
C,..)/2, and Cz/2 
checked. 


coefficients Cp, (Cp may be 
Also three In the 
second row, both C; and C;,/2 may be found and foyr 


found. zeros may be 


known coefficients may be checked. The Cy, values 
from the first and second row will agree exactly only if 
the data are perfect. In the third row, six of the seven 
coefficients may be found in terms of the seventh. In 
the fourth row six of the eight coefficients may be found 
in terms of the remaining two. If two linearly inde- 


pendent known forcing functions were introduced 
(e.g., rockets at two points fired at separate times), 
then u; would be zero for all equations and all coefii- 


cients could be determined. 
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A Simplified Method to Obtain the Load 


| Distribution Corresponding to the Ackeret 


Region for Wings Having Arbitrary Source 
Distribution at Supersonic Speeds 


JULIAN H. KAINER* 


Ryan Aeronautical Company 


SUMMARY 


A simplified approach to the loading on swept and triangular 
wings of planar and nonplanar form at supersonic speeds was 
presented!s? which required the potential in the region between 
the supersonic leading edge and the leading-edge disturbance 
wave. It is the purpose herein to indicate a simple means of ob 
taining these potentials and to justify the method with mathemat 
ical and physical reasoning. The process involved is inspec 


tion. 


INTRODUCTION 


rem THEORETICAL STUDIES’? have been com- 
pleted as a result of the current interest in the 
supersonic loading distributions due to aeroelastic 
and/or twist and camber effects on swept and triangu- 
lar wings. The means of applying the quasi-conical 
solutions’: * to the construction of the span loading on 
swept wings, in particular, requires a knowledge of the 
loading in the region bounded by the supersonic leading 
edge and the disturbance waves emanating from the 
root and tip chords. Hereinafter this region will be 
referred to as the Ackeret region, since it corresponds 
toa swept wing of infinite span. Experience has shown 
that the desired solutions are obtained by inspection of 
the existing solutions for the flow within the root and/or 


tip Mach cones. 


SYMBOLS 


V = free-stream velocity 

M = free-stream Mach Number 

8 = VM?—-1 

(x, y, 3) = Cartesian coordinates 

r = a point inside a Mach cone 

P, = a point outside a Mach cone 

u,v, w perturbations and/or components of veloc 
ity in the (x, y, s) directions, respectively 

or ® = velocity potential 

(, n) or (¢;, ») = Cartesian coordinates of a source disturb 
ance , 

Tp or Ip = source strength of a point inside or outside 
the Mach cone, respectively 

R, R; = area within the fore-cone of a point /? or 
P\, respectively 

6, = are cos [(1 — mt)/(m — t)! 
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are cos [(1 + mt)/(m +.t)] 


0» = 
0 = arc cosh (1/t) 
t = By/x 
€ = wing apex angle 
m = @tane« 
ANALYSIS 
Theorem 


Given the potential for a point P(x, y) within the 
after-cone of the leading-edge disturbance at 0, the 
potential for any point P,; on the wing exterior to the 
leading-edge disturbance is obtained by retaining only 
the pure real expressions when P is replaced by 


rr 


Proof 


Let u, v, and w be the perturbation velocities in the 


x, y, and g directions, respectively. Furthermore, de- 


fine the pure real potential ¢ (or ®), 


(1) 


¢ = f (u dx + vdy + w dz) 


Puckett* has shown that in the plane of the wing the 
potential at a point P(x, y) is represented by a summa- 
tion of all sources, (£;, »;) located within the fore-cone 


R of the point P; hence, 


~— / J 
JSR V (xp — 


where (x, y) represents the Cartesian coordinates of a 


(2) 


Tp dé dn 
£,)* — Bye — 2)? 


point within the leading Mach cone, (£;, ,) are the 
Cartesian coordinates of source disturbances in R, and 
Tp is the source strength at (x, y). 

Then, by definition, the potential at a point Pi(xp,, 
yp,) exterior to the leading-edge Mach cone is defined 


similar to Eq. (3) as 


[ [ Tp, dé dn : 
(3) 
J Ri V (xp, — &:)? — B*(ve, — 1:)° 


Eq. (2) is now expanded as follows: 


Pp, = 
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LEADING EDGE 








AREA OBPD=R 


AREA CAR=R, = 
* ai 
YN 
Vv - Xs 





i a a a a a a aa ae ae a a a 


Fig. 1. Typical configuration. 


2 F Ip dé dn 
| ] a = = (2) 
g V (ep — &)° — Pye — Hi)" 


Lf Tp dé dn 
ms Ri aN “(xp 
I, J V (xp — 


where all expressions under the integrals are pure real 
since (vp — &;) > B(yp — ni). 

The potential at P;(xp,, yp,) is obtained by changing 
proper subscripts in Eq. (4), thus 


Pp 








gr = 


— §)? — B*(yp — n,)? 
Tp dé dn 


E,)? — B(yp — 1)” 


(4) 


Tp, dé dn 


oP, sad J J / 9 
Ru \ (Xp, = £,)° 
> . Ip, dé dn , 
; \° ‘ "N\6 (9) 
E R= By Witte, = ° — POer — Hi )* 


where (£;, 7;) are sources located in R; (area CA P,) and 


— B*(yp, — 13)? 


(¢,’, n;’) are sources located in R-R,, where (vp, 
E,’) < Bye, — 1’). 

By examination it is noted that the expression im- 
mediately to the right of the equality signs in Eq. (5) 
is identically the expression ®p, defined in Eq. (3); 
hence, 

Tp, dé dn 


bp, + | [ ; - = 
R—Ry - (te = 8) — FOR ~- eT 


(6) 


op, = 


Since (xp, — &:’) < Blype, — 1’), Eq. (6) may be 


written in the form of the complex potential, 


F ff i Ip, dé dn 
op, = PP a fre 9 4\9 
R-—R V B2(yp, — ni’)? — (xp, — &:’)? 


(7) 
By physical identity, 


Pp, 


op, = 
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so the only mathematical means of expressing the po. 
tential at P, from the expression for ¢@p, derived from 
Eq. (7) is that 

real part ®p (9) 


op, = 


In order to secure physical consistency for the expres- 
sion just derived, it is necessary to prove only that the 
imaginary potential has no physical contribution to the 
potential at P;. This is immediately apparent upon 
consideration of the manner in which the disturbances 
in a uniform supersonic stream are propagated (see, 
Puckett’). A disturbance 
(x, y) in a uniform supersonic flow creates a stationary 


for example, located at 
conical wave front in which the effects of the disturb 
ances are felt. Hence, a point in a source field is dis 
turbed only by those sources that lie within the fore- 
cone of the point in question. It may be shown geo- 
metrically that the imaginary function corresponds 
to the source distribution (£,’, n;’) in the field (R 
R,) that does not lie within the fore-cone of P,(xp 


yp,). Therefore, it may be implied that the imaginary 
functions appearing in ¢@p, = function of ®p do not 


contribute to the velocity potential at P;. This con- 
stitutes a physical proof for the mathematical deriva- 
tion that the potential in the Ackeret region may be ob- 
tained by retaining only the real part of the three- 
dimensional potential. 

The form of the loading functions! ? within the 


Mach cone is quasi-conical—4.e., 


d = x'd(t) (10) 
where the conical functions ¢(¢) consist of various com- 
binations of polynomials in m and ¢ multiplied by the 
various antitrigonometric functions (i.e., 0; and 6s), 
various inverse hyperbolic functions (i.e., 63), and vari- 
In the Ackeret region 


9 


ous irrational functions in ¢.! ? 
(1 < ¢< m) all the above functions become complex, 
the real part of which is zero with the exception of 4, 
which becomes equal to 7. Hence, ¢(¢) reduces simply 
to that polynomial in (m, ¢) which is the coefficient of 


6, with 6; replaced by 7. 
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Airflow About Cone-Cylinders with 
Curved Shock Waves 


V. E. 


BERGDOLT 


Ballistics Research Laboratories, Aberdeen Proving Ground 


ABSTRACT 


Flows about 35° cone-cylinders in free flight at Mach Numbers 
ranging from the régime of detached shocks to the region where 
completely supersonic Taylor-Maccoll flow has been established 
have been investigated with interferometric techniques. Meas- 
ured fringe shifts from the interferograms, density distributions, 
and the corresponding pressures have been compared with theo- 
data. Good agreement 


retical predictions and wind-tunnel 


has been found with theory and other experiments. 


(1) INTRODUCTION 


bps FLOW ABOUT an infinite solid cone traveling 
through air at supersonic velocity has been deter- 
mined mathematically, with certain assumptions, by 
Taylor and Maccoll! and others.?, Depending on the 
free-stream Mach Number, there are three possible 
types of solutions: The flow behind the shock wave 
can be completely subsonic, subsonic on the cone and 


supersonic behind the shock, or wholly supersonic. The 
limits for these flow régimes are tabulated by 
Kopal? as functions of Mach Number and cone 
angle. 


For the cone-cylinder or any practical model with a 
finite cone, only the completely supersonic flow can be 
realized. Partially subsonic Taylor-Maccoll flow can- 
not occur on a finite cone because the influence of the 
change in model contour would be propagated upstream 
through the subsonic flow along the cone, thereby dis- 
rupting the conical nature of the flow. However, the 
shock wave remains attached down to the minimum 
Mach Number for Taylor-Maccoll flow but now is 
curved because of the influence of the corner. Below 
this Mach Number, the shock is detached. It was the 
purpose of this program to investigate by interfero- 
metric techniques flows over cone-cylinders in free flight 
from the region of detached shocks to the establish- 
ment of completely supersonic Taylor-Maccoll 
flow. 

In a supersonic wind tunnel at the 
tute of Technology, Mac Kinnon* determined pressure 
distributions over the surface of 35 
cylinder models. His test Mach Numbers ranged from 
1.49 for a detached shock wave to 1.997, where the 
flow is theoretically completely supersonic. In the pres- 
ent work, geometrically similar projectiles were fired 
at approximately the same Mach Numbers in order to 
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make comparisons between the wind-tunnel and inter- 
ferometric data. 

For the 35° cone, 
flow régimes are at Mach Numbers 1.681, 
1.911. The wind-tunnel tests were at Mach Numbers 
1.49 for a detached shock wave, 1.636 for a detached 
shock just before attachment, 1.694 for completely sub- 
sonic flow behind the shock wave, 1.86 for mixed sub- 
sonic and supersonic flow, and 1.997 for completely 
supersonic flow. The projectile models were fired 
at Mach Numbers of 1.49, 1.64, 1.69, 1.87, and 2.03. 


the theoretical limits of the three 
1.769, and 


(II) EXPERIMENTAL DETAILS 


The models were brass cone-cylinders incorporating 
a rotating band for spin stabilization and having a 35° 
semicone angle, Fig. 1. They were launched from a 
caliber 0.30 rifle in the enclosed Model Firing Range.*!” 
A velocity-measuring system consisting of photocell 
triggering beams and a 1.6 me. electronic counter gave 
velocity for each round fired. Interferograms, of which 
Fig. 2 at W/ = 1.49 is typical, were obtained with a 2-in. 
Mach-Zehnder interferometer and a rotating mirror 
camera system. The light source was a B-H6 mercury 
tube flashed by the discharge of a 0.07-microfarad con- 
denser charged to 13,000 volts and used in conjunction 
with a carbon-bisulfide liquid prism monochromator. 
This source has an effective duration of approximately 
5 microsec. during which the projectiles move about 
3 mm., but the projectile image is stopped on the film 
with the rotating mirror. The monochromator was 
adjusted to have a pass band of approximately 55A. in 
the region of the 4,358A. mercury line. This gave 150 
good fringes across the interferometer field. 
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(III) Data REDUCTION 


From fringe shifts measured on the interferogram, 
the air density can be computed throughout an axi- 
symmetric disturbance such as that generated by the 
cone-cylinder model.* * When this flow is observed 
by a Mach-Zehnder interferometer, the fringe shift at 
the point (y, s) on the film is related to the ratio of 
disturbed to free-stream density p/p; by 
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2K "tw ((per, z )— ] 
Oy, 2) = S / | st e dr (1) 
r ; Vr— y 


where (y, 2) are cartesian coordinates with the z axis 
on the axis of symmetry, A is the Gladstone- Dale con- 
stant (1.e., the ratio of increments of index of refraction 
to increments of density), A is the vacuum wave length 
of the monochromatic light source, and y = ry at the 
shock on the trace z = constant. The solution of this 
Abel integral equation is 


PG. 2) Kp "ty O67, 2)/OFr 
=]|]—- / dr (2) 
Pi wr J ¥ Vv r” = ¥y* 


For numerical work this can be manipulated into the 


form 
‘ N-1 
Pcr. 2) 2r bj41 — 5; 
= | — — 
pr tKa §=5 Vir? — re + Vee — 12 

(3) 

where r = 7; (j = 1,...,-N) isa monotone sequence of 

values on the relevant part of the trace, y = r;, and 

6; = 4(r;, 2). 


Fringe locations on the interferograms, Fig. 2, were 
measured with a Carl Zeiss toolmaker’s microscope. 
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Fringe shift values 6; were obtained along traces z = 
constant perpendicular to the flow axis at multiples 
of 1 mm. from the projectile tip, Fig. 3. The calcula- 
tion of density from fringe shifts is a lengthy operation, 
involving the determination of as many as 3,500 values 
for an accurate reduction of an interferogram. The 
actual computation has been performed on the ENIAC, 
a machine well adapted for such work. A detailed 
description of the interferometric data reduction proc- 
ess is given in reference 9. 

From a knowledge of the density ratio p/p, the 
pressure, the temperature, and the velocity can be de- 
termined by well-known relations from compressible- 
flow theory.° In the present investigation, ratios of 
surface pressure along the model to the free-stream 


stagnation pressure were required. Now 


p Pio = (p P20) (Pr Pro) (4) 


where p denotes static pressure and where subscript 0 
refers to stagnation values; 1, to upstream conditions; 
and 2, to the downstream side of the shock. 


p/ poo = (p/p) (5) 
where Poy and py» are constant on segments of stream 
lines between successive shocks. Also, 

P/ p21 = (p/ pi) (p1/ pro) (Pro, Pro) (6) 


By isentropic flow relations, 
pi/po = [1 + (1/2) (y — 1) A? ] 


where 1/7; is the Mach Number of the undisturbed 


1/(y-—1) a 
p (4) 


stream. Finally, by shock conditions, 
px P20 (y + 1) M,? sin? o What) 
= = er ee x 
Pro Pio (y — 1) A)? sin? o + 2 
1/(y¥-1) 


: se. 
27M)? sin? ¢ — (y — 1) 


(S) 


where o is the shock-wave angle. Values for Eq. (7) 
are tabulated as functions of J/, and for Eq. (8) as 
functions of 1/, sin o by Emmons’ and in reference 11. 
It will be assumed that the modei profile is a stream line 
along which these formulas are applicable. 


(IV) EXPERIMENTAL RESULTS 


(A) Conicity 

For the Taylor-Maccoll or conical-flow model in the 
region bounded by the cone and the shock wave, ve- 
locity components, density, and pressure do not vary 
along straight lines emanating from the vertex of the 
It follows that, if the fringe shift 6 is measured 
along traces per- 


cone. 
as a function of radial distance Y 
pendicular to the cone axis at various values of z from 
the tip, then the points (V/z, 6/z) should fall on a 
single curve.’ The degree to which these points pro- 
duce a single curve is a measure of the conicity or ap- 
proach to conical flow. 


CONE-CYLINDERS 753 





CONICITY 
30 35° CAL. .30 CONE CYLINDER 038 
~™ | M*2.03 
wo | 
; | 
a © Troce No | 
+ Troce No 2 
* Tr N 3 
ici 1 Trace No 4 
£ t=, * Trace No 5 
u ss =e * Trace No 6 
a 
b 4 eo «6 We, ° 
o 8 » . ‘ = 
> ce » G0 
> ee Se 
~ 
ba 1O- ' - 3 
ie 3: 
=? Theory 
w 
' 
¥o 
t 
| . 
o — in 4 1 ee 
6 8 10 J2 14 
Y 
Reduced Radial Distance, /2 
Fic. 4. 


Figs. 3 and 4 show the development of conicity, 
ranging from the widely scattered data for definitely 
nonconical flow at JJ = 1.49 to approximately conical 
flow at \J = 2.03. The radial distance } is measured 
from the cone axis, and the traces terminate at the 
shock wave where the fringe shift becomes zero. Ap- 
preciable convergence has occurred when the shock wave 
becomes attached, but there is no pronounced tendency 
toward conicity until the limit for completely super- 
sonic Taylor-Maccoll flow is approached. The plot for 
M = 2.03 shows good convergence and agreement 
with theory, except for the last trace which appears to 
reflect some influence of the corner. 


(B) Isopycnals 


For the pressure distribution along the model sur- 
face, it would suffice on each trace to find the density 
only at the body. It is desirable to have the density 
over a substantial part of the trace in order to detect 
errors in the calculation. If these are isolated points, 
which is the usual case, the density values are graphic- 
ally smoothed to eliminate the errors. The calculation 
time is not increased appreciably if the entire trace is 
computed, so all reductions are carried through to 
produce the complete density field, Figs. 5 and 6. 

The conical-flow region for the cone-cylinder model is 
bounded by the cone surface, the shock wave, and the 
characteristic from the shoulder. Beyond this char- 
acteristic is a region of expansion produced by the flow 
turning the corner. If Taylor-Maccoll flow exists 
within the conical region, the isopycnals should theo- 
retically be straight lines emanating from the cone 
vertex. Fig. 6 shows this general radial pattern in 
the conical region, to within the experimental error, 
and the expansion that occurs around the corner. At 
the Mach Numbers below the supersonic régime, the 
radial pattern changes to one where the isopycnals 
gradually loop back to the cone surface with no appear- 


ance of conicity. For these cases with the curved shock 
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wave, it is expected that the flow will attain sonic 
velocity at the corner and that the expansion region 
will develop as in the supersonic case. This expansion 
is shown in all the isopyenals, as well as an early ex- 
This effect 
The high-density 


pansion that occurs ahead of the corner. 
is discussed in the following section. 
gradient associated with the corner expansion results 
in a high fringe shift gradient, which can be seen on the 
interferograms. 

For conical flow, the density along the cone should 
be constant, since the surface defines a straight line 
from the cone vertex. At Mach Number 2.03, Fig. 6, 
the density is not actually constant. Because of prac- 
tical limitations in the construction of the models, the 


Distance in 


Calibers 


5. 


tip is not perfectly sharp so that the flow encounters 
a truncated or rounded cone which means a detached 
The detachment is so slight that it is 
There are 


shock wave. 
obscured by the relatively wide fringes. 
reduction errors’ and probably yaw effects contributing 
to the departure from constant density. These effects 
are manifest in the slightly increased density near the 
tip, which has also been observed independently by 
Ladenburg.° 

Of immediate interest is the location of the sonic 
line in the subsonic flow cases. Along the stream lines, 
the velocity and pressure are specifically related to the 
density by suitable forms of Eqs. (5) and (7). Since 
the stagnation values of density vary from one stream 
line to another, it is necessary to have the system of 
stream lines in the flow. These stream lines can be 
found from the density by a semigraphical method" 
based on the equation of continuity for a stream tube 
bounded by two stream lines. In Fig. 5 the stream 
lines and the sonic line have been superimposed upon 
the isopyenals. The sonic line intersects the cone just 
ahead of the corner at a slant height of 0.98 from the 
vertex. On the pressure distribution plot of Fig. 7, the 
sonic value intersects the 031 curve at an X/S of 0.97. 

The reduction of the interferometric data to densities 
is based on the assumption of zero yaw. This condi- 
tion is usually not realized, although projectiles can be 
fired with yaws of 20 to 40 min. Interferograms are 
considered suitable for reduction if the yaw is approxt- 
The two halves of the 


mately half a degree or less. 
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flow then are not identical and are considered as two 
separate flows, each representing an axisymmetric 
model. If the interferogram of Fig. 2 is studied, it is 
evident that the accuracy of film measurement differs 
for the two halves, since the fringe pattern is more 
favorable on the top side. The difference in fringe 
shift values is such as to prevent averaging the two 
sides. In Figs. 7 to 9 the actual divergence of the two 
halves can be seen for the pressure distributions. 


(C) Pressures 


The general transition of the flow from the detached 
shock wave to the supersonic condition is shown in the 
pressure distributions of Figs. 7 to9. For Mach Num- 
ber 1.49, the shock wave is detached, and there is no 
resemblance to the constant pressure distribution for 
Taylor-Maccoll flow. At the cone vertex the value is 
that for the stagnation pressure behind a normal shock. 
The rapid drop in the pressure curves indicates the 
rapid increase of velocity necessary for the flow to at- 
tain the predicted sonic velocity at the corner. Since 
the flow is theoretically sonic at the corner, there 
should be a Prandtl-Meyer expansion here. The ex- 
perimental values appear to be in good agreement with 
the predicted expansion, but the expansion begins 
ahead of the corner and the sonic point likewise moves 








ahead. This might be expected, since it is conjec- 
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tured that the boundary layer builds up along the cone 
and rounds off the corner, so that the corner seems to 
move forward and the flow does not actually follow the 
cone-cylinder profile. The early expansion of the flow 
at the corner was observed in all the cases. 

The character of the pressure distribution changes 
The 
pressures encountered near the tip are reduced, since 
the tip is no longer a stagnation point, and the curves 


when the shock wave becomes attached, Fig. 8. 


begin to assume the relatively flat appearance that is 
Taylor-Maccoll flow model. As 


r 


associated with the 
the Mach Number increases toward the minimum for 
completely supersonic flow, the theoretical value for 
the cone surface pressure is approached. At Mach 
Number 2.03, Fig. 9, the flow should be completely 
supersonic along the cone, and the experimental results 
show that the pressure is virtually constant. There ap- 
pears to be good agreement at all Mach Numbers with 
Mac Kinnon’s results, as well as with the predictions of 
theory. 


(D) Head-Drag Coefficient 


The head drag on the model is a function of the pres- 
sure distribution over the cone surface. A form of head- 


drag coefficient is 


D " , X x 
Kn = — = se : (*) a(5) we 
pid7v? 2p? Jo po \S S 
MAVYCtY #7! x xX 
2 (1+ ) | P (3) «(5) (9) 
27M, 5 0 po \S mo 


The five values shown in Fig. 10 were computed by this 
definition. For Taylor-Maccoll flow, Kopal’ uses the 
usual definition based on excess pressure, p; — /, to 
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, Ds a," b. l 
K pxopet 7 E _ 1) - - (/ a .) ® (10) 
Pi a" Pi yM* 


where p, is the constant pressure along the cone surface, 
With the constant cone surface pressure, the relation 
between the two definitions becomes 


Kp = (1, 4) [AK pxopat + (1 yM,?) | (11) 


The Kopal (modified) curve was plotted from Eq. 
Gane 

Below the completely supersonic régime, Fig. 10, the 
experimental values do not agree with theory, of course. 
However, the curve fairs smoothly into the theoretical 
plot in the supersonic régime. The failure of the two 
curves to be tangent at the supersonic flow limit is 
probably due to the pronounced rounding off of the 
pressure distribution curves near the model corner, 
which causes a reduction in the drag coefficient. In 
this particular case it should be remarked that extra- 
polation along the tangent to the theoretical curve 
at the completely supersonic limit will lead to values 
of Kp which are too high. 
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Two-Dimensional Unsvmmetric Flow 
Patterns at Mach Number 1 


GOTTFRIED GUDERLEY* ann HIDEO YOSHIHARAt 
Wright Ltr Development Center 


ABSTRACT 


Two-dimensional unsymmetric flow patterns at Mach Number 
1 are discussed by means of the hodograph representation. As 
an example, the flow over a wedge profile is investigated and the 
pressure distribution and the lift curve slope are computed. The 
results indicate that the value of the lift curve slope at Mach 
Number 1 is much lower than what would be expected by an 


extrapolation from the purely supersonic results 


INTRODUCTION 


PATTERNS with the free-stream Mach 


previously are all sym- 


"T Vu FLOW 
Number 1 


metric with respect to an axis parallel to the free-stream 


considered 
direction. In the present paper we shall extend these 
considerations to unsymmetric flows; in particular, we 
shall investigate the dissymmetry introduced by placing 
a symmetric wedge at an angle of attack. 

The investigation will be carried out in the hodograph 
Here the point that represents the free-stream 


plane. 
This singularity has been 


velocity will be singular. 
discussed for symmetric flow patterns, but in the pres- 
ent case it must be re-examined. 

Knowing the properties of this singular point, one 
can formulate the desired boundary value problem. 
In order to obtain the general character of unsymmet- 
ric flows, we shall first discuss the case of a body with 
a rounded nose, since a profile with a pointed leading 
edge (for example, a wedge) will have additional com- 
plications at an angle of attack because of a local super- 
sonic region arising at the nose. The presence of this 
region is not a typical characteristic for a general pro- 
file with an angle of attack, and it adds the difficulty 
that the boundary value problem cannot be formulated 
because of the complexity of the flow over the sharp 
nose. On the other hand, a profile with a rounded lead- 
ing edge gives rise to a difficult boundary condition in 
the hodograph plane when the profile is prescribed in 
the physical plane. To overcome this difficulty, we 
shall consider the flow over the body to be known for 
a given configuration—for example, at zero angle of 
attack. (This flow may be computed indirectly start- 
ing from the hodograph plane.) We then ask for the 
change of this flow due to a small angle of attack. In 
this case, the boundary conditions can be linearized 
to a convenient form. By these investigations, the 
Received March 12, 1953 
* Aeronautical Research Laboratory. 
t+ Wind Tunnel Branch, Aircraft Laboratory 


general characteristics of an unsymmetric flow can be 
found. It will be seen in particular that the flow pat- 
tern is entirely determined by its boundary conditions 
that is, no arbitrariness arises as in the case of incom- 
pressible flows where one required the Kutta condition 
to insure uniqueness. 

The actual computations are carried out for a wedge 
in this regard the present paper may be con- 
sidered as a continuation of The effect 
of the local supersonic region which must be expected 
at the nose will be shown to give a negligible contribu- 
tion to the lift for sufficiently small angles of attack. 
Under these circumstances the wedge has the advantage 
that the boundary value problem in the hodograph 
plane can be formulated and solved with comparative 
ease. One will finally obtain the change of the pres- 
sure distribution over the wedge due to the angle of 
attack. In this way the slope of the lift curve and the 
location of the center of pressure may be found. 


profile ; 
reference 1. 


SECTION (I)-——Basic RELATIONS 


Let us introduce the following symbols: x and y are 
the cartesian coordinates in the physical plane; w and 
6, the magnitude and inclination of the velocity vector; 
p and p, the pressure and density of the gas; a, the ve- 
locity of sound; W, the stream function; &, Legendre’s 
transformed potential; and x, the ratio of specific heats. 
The asterisk (*) will denote the conditions of the sonic 
state. 

The flow is assumed to be such that the usual tran- 
sonic approximations are valid (see reference 2). Under 
these conditions, the hodograph equations for the stream 
function and Legendre’s potential may be simplified to 


Vin — 1 Von = O (1) 
Ou — n Pe = 0 (2) 

where 
n= (x + 1)” [(w — w*)/w*] (3) 


Eqs. (1) and (2) are of the form investigated by Tri- 


comi.* In the subsequent sections, we shall refer to 
either of them as the Tricomi equation. 

The transformation equations from the hodograph 
plane to the physical plane assume an _ especially 
simple form in terms of the Legendre potential; one has 
(4) 


; l ‘ 
x = (x + 1)  P(0, ) 
aw 


4o¢ 
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Fic. 1. Physical plane for the general body. 


y = (1/w*)-, (n, 8) 


Here the subscripts denote partial differentiations. 
The corresponding expressions in terms of the stream 
function may be obtained by using the relation 


Vv = p*-9, (5) 


For the supersonic region, the characteristics and the 
compatibility relations are given by 


6+ (2 3)n”? = const. (6a) 
d(,) — -(dn/d0)-d(%,) = 0 (6b) 


In the physical’plane, the direction of the Mach waves 
will be 


—I/, = 


dy/dx = + (kx + 1)~'7?-n (7) 


Here the signs must be used consistently with those in 
Eq. (6a). 

Consistent with our approximation, the pressure 
coefficient, C,, is expressed in the form 


C, = (p — p*)/(ptw*?/2) = —2(K + .1)7 +n 
(8) 


SecTION (II)—TRE FLOW OVER A SYMMETRIC Bopy 
WITH A ROUNDED NOSE 


(A) The Hodograph 


The general flow configuration to be expected for a 
body with a rounded nose is shown in Fig. 1 with the 
corresponding hodograph plane shown in Fig. 2. Since 
all stream lines must originate far upstream of the pro- 
file with the free-stream conditions, they must all have 
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their origin at the single point A in the hodograph 
plane representing the free-stream velocity. This point 
will, therefore, be a singular point for the stream func- 
tion; we shall call it the free-stream singularity. 

Let us trace the course of the zero stream line in 
order to orient ourselves in the hodograph plane. It 
starts at point A and proceeds toward the profile with 
decreasing velocity. It reaches the body at the stag- 
nation point at B where the stream line branches. The 
course of these branches is then shown by the shaded 
lines. Along the upper branch, the velocity will now 
increase until at some point C the sonic velocity is 
reached. The velocity will then become supersonic 
along the segment of the zero stream line given by CD. 
One has a similar behavior along the lower branch. 

Of the Mach waves starting from the body, one will 
have two classes; the first class is composed of those 
that reach the sonic line at a finite distance from the 
body, while those of the second never reach the sonic 
line (see Fig. 1). The distinction between these two 
classes is made in order to determine the portion of the 
supersonic region which influences the subsonic flow. 
It is clear that the subsonic flow cannot be completely 
determined by specifying only the contour of the body 
up to the sonic line. The segment of the profile from 
which the Mach waves of the first class originate must 
also be given, since a perturbation of this portion of the 
contour will propagate along the Mach waves of the 
first class to the sonic line and thus influence the sub- 
sonic flow. On the other hand, perturbations travel- 
ing along the Mach waves of the second class will never 
reach the sonic line and will not affect the subsonic 
flow. We must accordingly prescribe the contour of 
the profile downstream to the point from which the 
last Mach wave of the first class starts. This wave 
will have a special significance in our discussions, and 
we shall label it as the limiting Mach wave. Thus, in 
Fig. 2, we have shown the zero stream line only to 
points D and F. The limiting Mach waves are the 
characteristics A D and AF. 

The continuation of the flow downstream of the limit- 
ing Mach wave may be computed independently by a 
purely supersonic method—for example, by the method 
of characteristics. Of course, the required initial con- 
ditions will be furnished by the solution for the flow at 
the limiting Mach wave, while the necessary boundary 
conditions are furnished by prescribing the remaining 
contour of the body. Incidentally, for a closed body, 
this portion of the flow will occur on several Riemann 


sheets in the hodograph plane. 


(B) The Free-Stream Singularity 


There exist particular solutions of Tricomi’s equation 
which are singular at one point of the sonic line.? In 
the present section, we shall determine which of these 
particular solutions is suitable to represent the free- 
stream singularity mentioned in the previous section. 
For the stream function we use the hypothesis 
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n 


¥,. = pn SC, #,) (9) 


ns 


where 
¢ = 962 /4n? (10) 


and for the Legendre potential 
d,, = in S(F, ty) (11) 


Here is a constant with the subscripts s and p de- 
noting whether the » refers to the stream function or 
to the potential. The function f(¢, 7) must fulfill the 
hypergeometric differential equation 

2 l 


i = 1) -f” + i- io = “£5 = zs f? 4+ 
oO 


I 
nin — lyf = 0 (8) 

QQ 
Because of Eq. (5), ¥ = p*-%,, particular solutions in 
terms of VW and ® are closely connected. We have es 


pecially the relations 


3(4n, — 5) ,. sen : r 

fb, ".) = 5 77+ f'(E, my) sign (n-@) (13) 
2hgitt, — 1) 

where the symbol, sign 7-0, means a factor one with the 

sign of 7-6. Because of these relations, we may con- 

fine our subsequent considerations to only the stream 


function. 


Before we proceed to determine the parameter 1, let 
us first discuss the solutions of the hypergeometric 
equation, Eq. (12). The singular points of the equa- 
tion are located at the values of ¢ of 0, 1, and o™. 
About each of these singular points there are two 
linearly independent power series solutions. For each 
value of ¢, we find a line in the yn, 6-plane since ¢ = 
(96?) ‘(4n*). Therefore, about each of these lines 
there will be two linearly independent series solutions 
of the hypergeometric equation. If we choose the 
representation about one of the lines, then we must, 
from physical reasoning, determine the coefficients of 
the two linearly independent particular solutions about 
each of the other lines such that the stream function 
and its first derivatives are continuous in the entire 
region under consideration. 

The choice of the solution f must be dictated from 

It is clear that the singular points 
» are due solely to the form of the 


physical reasoning. 
at¢ = Oand ¢ = 
independent variable ¢; the singularity at ¢ = 1, on 
the other hand, is essential. In Fig. 2, if we take the 
origin of the hodograph plane at point A, the lines ¢ = 
1 will represent the limiting Mach waves mentioned 
in a previous section. . It is therefore logical that the 
selection of the correct representation of the function f 
be made by examining the series solution about this 
line. In the vicinity of this value of ¢, one of the line- 
arly independent solutions of Eq. (12), denoted by 
u;({, m), is given by a power series in (¢ — 1) while the 


on 
> 
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other, denoted by 73(¢, ”), is given by a power series 
multiplied by a factor (¢ — 1)°"*"°. Accordingly, 
the expression for Y formed with the portion u; is 
regular along the characteristics ¢ = 1, while that with 
v3 will tend to infinity for negative values of m, the 
negative values of m being necessary to produce a singu- 
larity at the free stream point. The presence of a 
portion of the solution with v3 will cause the character- 
istics ¢ = 1 to map to infinity in the physical plane. 
[See Eqs. (4) and (5) giving the transformation. | 
Since the characteristics ¢ = 1 correspond to the limit- 
ing Mach waves, we are naturally interested only in 
such cases for which the limiting Mach waves are not 
mapped completely to infinity in the physical plane. 

With these remarks the permissible values of m may 
now be found. According to the previous discussion, 
we stipulate that the function f({, ~) must be given by 
u3(¢, mn) along the limiting Mach waves AD and AF 
of Fig. 2. It is found that a solution f(¢, ) for an 
arbitrary value of » which is represented along AD by 
u3({, a) is given along the other characteristic A F by a 
linear combination of “3; and v3, where the coefficients 
are functions of » (see reference 2). The condition 
that the coefficient of v; vanishes along AF will give a 
set of permissible values of 1; we have n, —-|- 
(3/2)k, where k = 0,1, 2,.... 

If we now investigate the physical plane resulting 
from the particular solutions with these values of n,, we 
will find that a single-valued physical plane will result 
only for values of n, equal to —5/2 and —1. We may 
therefore eliminate all but these values of »,. Of these, 
the singularity produced by the particular solution for 
n, = —95/2is familiar; it is this that was used to repre- 
sent the behavior of the flow at a far distance from the 


airfoil for the case of zero angle of attack. The par- 
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Fic. 2. Hodograph plane for the general body. 
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ticular solution with , = —1 taken by itself does not 
give a flow that is of interest at present. It may, how- 
ever, still be used superimposed to the particular solu- 
tion with m, = —5/2, since, near the free-stream point, 
the latter particular solution will predominate. There- 
fore, the free-stream singularity can be represented by 
the expression 


a: (5/2) + Co°V_, (14) 


where c; and ¢, are ‘wo arbitrary constants and W,,, is 
defined in Eq. (9). 

The corresponding representation for the free-stream 
singularity in terms of the Legendre potential may be 
found by using the relationship between , and n, as 
given by Eq. (13). We have, therefore, nm, = —1 and 
1/,, The particular solution with , = '/»2, however, 
is. no longer singular at the free-stream point in the 
sense of Tricomi’s existence and uniqueness proofs.’ 
Consequently, in terms of the Legendre potential, the 
free-stream singularity is represented by only one 
singular particular solution, ®_;. Thus, there is but 
one arbitrary constant in this representation. 


(C) Boundary Value Problem for the General Profile 


As mentioned previously, the problem for a profile 
with a rounded nose will give a better insight into the 
mathematical and physical questions of unsymmetric 
flows, but it will assume a rather difficult form in the 
hodograph plane. For this reason, we must restrict 
our considerations to perturbations of a known flow. 
In this way, a linear boundary value problem will be 
obtained. We shall assume that the flow over a cer- 
tain symmetric profile is known, and we want to find 
the influence of a small angle of attack. For this pur- 
pose, it is convenient to keep the position of the body 
fixed in space and to rotate the free-stream direction. 
In the hodograph, the free-stream singularity will then 
be shifted from the point (y = 0, 6 = 0) to the point 
(n = 0, 6 = a), where a is the angle of attack. In 
order to obtain separately the change of the flow pattern 
due to the angle of attack, we ask which hodograph 
solution must be superimposed to the hodograph repre- 
sentation of the original flow with a = 0 such that the 
flow with an angle of attark is obtained. 

We must therefore subtract from the original solution, 
an expression ®_;(n, #) [see Eq. (11) ], which contains the 
free-stream singularity for zero angle of attack, and 
then add another expression, &_,(7, 6 — a), which has 
the same singularity at the point (7 = 0,0 = a). The 
coefficient with which this latter particular solution is 
multiplied, however, must be considered as a function 
of a, the angle of attack. We have, therefore, for 
the particular solution that performs the shift of the 
free-stream singularity, the expression 


¢g = _,(n, 0) — [1 + k(a)] b_i(n, 0 — a) (15) 


where k(a@) is some function of a and ¢ is the Legendre 
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potential which represents the change of the free- 
stream singularity due to the angle of attack. For 
sufficiently small angles of attack, the second term in 
Eq. (15) may be expanded with respect to a, with the 
result 


g = —k(a)-_,(n, 0) + 
[1 + k(a)]-P_16-a@ + O(a?) 


This expression describes the flow at a far distance 
from the body due to the angle of attack; it must 
therefore be antisymmetric in a. The term _, is an 
even function of 6; consequently, in order for the above 
expression to have the correct symmetry, we must 
take k(a@) proportional to a®. We have, therefore, as 
the terms of lowest order in a, 


g¢=a:P_i4 
This, according to Eq. (13), is equivalent to 
¢g¢ = —(4/9a-P_ (5/2 (16) 


The boundary condition for the problem due to the 
angle of attack may be found as follows: The profile in 
the physical plane is first described in the parametric 
form 


x = p,(0)-w*-! (x + 1)’ 
y = p(0)-w*-! 


where /; and ~2 depend upon the profile shape. Using 
Eq. (4), one then obtains 

®, (IT, 6) = pi(9)\ a7 

i) 


II 


(17, 0) = p.(6)h 


Here ® is Legendre’s potential, and 7 = F7(@) is the 
hodograph of the profile, both for the case of zero angle 
of attack. [The fact that the function /7(@) is an un- 
known function will prevent us from considering the 
general case; we have assumed here that it is a given 
function.] By the change of the free-stream direction, 
the Legendre potential will be perturbed from # to 
@ + gy, where ¢ is due to the angle of attack. The 
hodograph of the body will also be changed from /7 to 
IT + 4. 


physical plane, we must have, therefore, 


In order for the profile to be fixed in the 


,(/7 + 4,0) + 9,77 + 4, 0) = pil(O) 
bo(H7 + 4, 0) + 9, (1 + 4, 0) = p2(8) 


For sufficiently small angles of attack, we may expand 
these expressions with respect to 4, and using Eq. (17), 
we obtain, as the terms of lowest order, 


4 ®,,(77, 0) + ¢,(77, 0) = 0 
7 P,9(17, 0) + ¢,(T7, 0) = 0 


If we eliminate 7 between these two equations, we have 


7 ,,, (7, 8) ‘ 
¢,(7, 0) — = -go(l1, 02) = O (18) 
Px 9 (17, 6) 


Moreover, the contour 7 = //(@) is a stream line; we 
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have, therefore, according to Eq. (5), 
$,(/7, 0) = 0 

or differentiating 

oo(H, 0) + (dH /d0)-,9(H, 0) = 0 
By using this relation and the Tricomi equation 

®,,(17, 0) — H-9(TT, 0) = 0 

one may rewrite Eq. (18) to 

¢,(H1, 0) + H-(di7/d0)-¢ (17, 0) = 0 (19) 


This is now the boundary condition that must be ful- 
filled by the potential ¢ due to the angle of attack. It 
is to be fulfilled along the original hodograph of the 
body 7 = /7(@) for zero angle of attack. 

The particular solution, ®_,5/2) [see Eq. (16) ], which 
gives the required free-stream singularity will not gen- 
erally fulfill the boundary condition (19) at the body. 
Therefore, one must superimpose additional regular 
particular solutions so as to correct the boundary condi- 
tions. The boundary value problem for determining 
the required regular portion of the solution is of a type 
that has not been investigated thoroughly. A few re- 
marks, however, regarding the plausibility of existence 
and uniqueness of this type of problem have been given 
in reference +. (Here, this type of problem has been 
called the boundary value problem of the second type 
because of its similarity to the problem of the second 
type in the potential theory.) According to these re- 
marks, it appears that such problems do not possess a 
solution (at least no physically useful solution) unless 
the boundary values fulfill an additional condition. 
This condition will in general determine the expression 
k(a) in Eq. (15) which, in the present case of a sym- 
metric body, was determined by consideration of sym- 
metry. 


SECTION (III)—-THE WEDGE PROFILE 


(A) The Boundary Value Problem in the Hodograph Plane 


For the actual computations, a wedge profile will 
be used. Naturally, at a far distance from the body, 
the considerations of the previous case will hold; but 
near the body, especially at the nose of the wedge, one 
will have essential differences of the flow. A qualita- 
tive picture of such a flow is shown in Fig. 3 with the 
corresponding hodograph shown in Fig. 4. The main 
difference of this problem in comparison to that for a 
profile with a rounded leading edge is seen along the 
upper branch of the zero stream line. The zero stream 
line proceeds from the stagnation point at B upstream 
along the contour of the wedge to the nose, where a 
Meyer expansion produces either zero absolute pres 
sure in the case of cavitation or some other pressure, 
probably supersonic, in the case of a separation of the 
flow. If the angle of attack is not too excessive, the 
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flow will again reattach, with the zero stream line re- 
turning to the wedge at some point along the side 
DE. A real satisfactory formulation of the boundary 
value problem in the nose region is probably impossible. 
However, it is shown in the Appendix that, if one re 
quires simply that the boundary condition (19) is ful- 
filled along the entire line BDE and disregards the nose 
region entirely, the contribution of the flow at the nose 
to the total lift is negligible for sufficiently small angles 
of attack. 

Disregarding the flow over the nose we have, there- 
fore, the problem of finding a solution of the Tricomi 
equation which has the required singularity at the free- 


stream point, as given by ®_,5/2), and which fulfills 
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the condition (19) along the sides of the wedge and 
along the characteristics representing the Meyer ex- 
pansion at the shoulder of the profile. 


In the present case, the boundary condition (19) 


can be simplified. Along the subsonic boundaries, @ is 
a constant, so that, here, the boundary condition be- 


comes 
go = O (20) 


The supersonic boundary is formed by characteristics 
that are at the same time stream lines. Therefore, the 
image of this boundary will map to single points in the 
physical plane. The boundary conditions along the 


characteristics will thus simplify to 


¢o = 0, ¢, = 0 (21) 


In the case of the wedge, the boundary value problem 
can also be formulated in terms of the stream function. 
Let us denote by y the stream function that must be 
superimposed to the hodograph of the original flow in 
order to obtain the flow with an angle of attack. The 
necessary boundary condition for y can now be ob- 
tained by using the relation Y = p*g, in Eqs. (20) and 
(21). We have, therefore, the condition y = O along 
the given boundary. Because of the symmetry of y, 
we may restrict our considerations only to the upper 
half of the hodograph plane. The required condition 
along the line of symmetry is then y, = 0. 
boundary conditions are illustrated in the upper half of 


These 


Fig. 5. 


The free-stream singularity for the complete hodo- 
graph (original flow plus the superimposed flow) was 
given by Eq. (14). The stream function, y, repre- 
senting the change of the free-stream singularity due 
to the angle of attack, will be found in the same manner 
as for the Legendre potential. We have 


4 
y= - [wv 5/2(m, 9 — 


9g’ 


i3/2)(7, 9) — (1+ K(a)-¥ 
a)] + ko(a)-V_i(n,8 — aw) (22) 


Here the coefficient c; has been taken as —4/9 in order 


to obtain in the present paper the same scale in the 


NOVEMBER, 1953 


SCIENCES- 
physical plane as in reference 1. The functions k,(q 
and ko(@) have been introduced, since the coefficient of 
the singularities at the point (yn = 0, 6 = a) are stil] 
undetermined. If in Eq. (22) we now carry out a de- 
velopment with respect to a and apply considerations of 
symmetry as was done for the Legendre potential, we 
obtain 

y = [(14/27)¥ 


(9, 0) + BW (9, 0) la (23 


where @ is a constant independent of a. 

One thus obtains a boundary value problem for y 
which contains an arbitrary parameter 8. For a fixed 
value of 8, the solution to this boundary value problem 
will exist and will be unique according to the proofs of 
Tricomi.* 
(21) that ¢, = O along the supersonic contour—i.e., 


In this formulation, the statement of Eq. 


that the superimposed solution does not change the 
value of x at the shoulder of the wedge—has not yet 
been expressed. The fulfillment of this requirement 
will determine the value of the parameter £. 

The future computations will be based upon the 
formulation of the boundary value problem in terms 
of the stream function. Here one has the firm founda 
tion of Tricomi’s theory. Furthermore, the particular 
solution ®:,,, which in terms of y appears as V 
in this form, less convergence difficulties, since the 


1, offers, 


singularity then is treated separately and enters into 
the solution of the boundary value problem only in the 
We shall, 
on the other hand, introduce the Legendre potential, 


form of the boundary values it induces. 


yg, during the computations because of its convenience 
in performing the transformation of the results to the 
physical plane. 


(B) Fulfillment of the Boundary Conditions 


The boundary conditions will be fulfilled in two steps. 
We consider first the boundary condition along the kori- 
zontal lines, 06 = 0 and @ = 6 (where & is the half-angle 
at the nose of the wedge), disregarding the condition 
along the characteristic KE’. This portion of the solu- 
tion will contain the required free-stream singularity. 
The second step will be to superimpose regular particu- 
lar solutions that fulfill the remaining boundary condi- 
tions along the characteristic KE’ without disturbing 
those that were considered in the first step. 

Let us consider the conditions along the horizontal 
lines. These conditions can be considered as symmetry 
requirements that may be fulfilled by a suitable system 
of images of the free-stream singularity placed along 


the @-axis. This portion of the solution is given by 


— 


14 , — 
p* ry = [3 a 2% (O's. ~ 0 + 


=é m 


Bini! SS (-1)"-fEm, — 1) |-< (24) 


m 
where 


Cw = 9(6 7 2m6)? tn (25) 
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Here the factor (—1)” is required, since the functions 


f(tm, — 4) and f({m, —1) are symmetric with respect to 
the line @ = 2,0. By examining the asymptotic be- 
havior of f({m, —4) and f({m, —1) for large values of 


the argument, one can show directly that the series in 
Eq. (24) converges uniformly in the region of the 
hodograph plane in which the boundary value problem 
is defined. 

We shall next superimpose regular solutions of the 
Tricomi Eq. (1) which leave the boundary conditions 
along AB and DE undisturbed but which give the 
correct boundary values along EE’. The regular par- 
ticular solutions to be used for this purpose correspond 
to Chaplygin’s particular solutions, and they are given 
by 

y’"” = g(n, A)-cos (A- 70/8) (26) 


where \ = m + (1/2) and where 


V |n| 0H), (ad) (n < 0) 


9 
72 VU nlJi/,(A8) + J -1/,(A8) ] (n > 0) 
gin, A) = (V3 
rs (7) - (n = 0) 
” (1/3)! 6 ¥ 
with z = (27/36)\n| ‘*. Here //,,,“” is the Hanker func- 


tion of the first kind of order one-third, and /+1/, is the 
Bessel function or order +'/;. This set of functions 
is chosen such that the nose of the wedge is not mapped 
to infinity in the physical plane and that the stream 
function is zero along the lines 6 = 6) and yz is zero for 
6 = 0. We shall now superimpose to the expression 
for the image system, Eq. (24), a linear combination of 
the regular functions given by Eq. (26), with suitable 
coefficients a,, and 0,, such that the remaining boundary 
conditions along F’G are satisfied. The stream function 


will therefore become 


97 


14 
p* Ly = ‘an ‘ >> (—1)"-f(Em — 4)+ 


7 anv | aA | 1 YS (-1)"f(Ems, — 1) + 


m=0 
(m) > 
bnew | (27) 
m=0 


Eq. (27) must be evaluated along the characteristic 
EE’ and the coefficients a,, and b,, must be determined 
in such a manner that y = 0 along this characteristic. 
Since up to the present 8 is unknown, this condition may 
be fulfilled regardless of the value of 8. Thus one ob- 
tains the two equations 

ni-* > (-1) "fm —4) + ¥ anv” = 0 (28) 
m 


m 


where y'”’ is defined in Eq. (26). 


2d (-)D "fm — D+ YE dnw’” = 0 (29) 


m 
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Fic. 6. Fulfillment of the boundary conditions along the Meyer 
expansion 


n = [(3/2)(% — 0)]” (30) 


In the subsequent discussion we shall restrict ourselves 
to the case of the coefficients ),,. 

Consider the boundary value problem for the deter- 
mination of that portion of the solution which is con- 
tained in the second series within the bracket expressions 
of Eq. (27). According to Tricomi,’® the solution of 
problem exists and is continuous. Therefore, it can be 
represented uniquely in the subsonic region by a series 
of Chaplygin solutions as given by Eq. (26). This 
series will converge throughout the subsonic region in- 
cluding the sonic line. In the supersonic region, an 
examination of the Bessel functions (26) in their asymp- 
totic form will show that the stream function will con- 
verge at least to the same degree as it did along the 
sonic line. Therefore, the representation of the stream 
function in terms of the above series of Chaplygin 
solutions will give the desired continuation of the sub- 
sonic solution into the supersonic region. In particu- 
lar, this argument shows that there will exist coefficients 
b,, in Eqs. (32a) and (32b) such that the given terms on 
the left side of the equations may be represented arbi- 
trarily close if one takes enough terms of the series on 
the right-hand side. 

It should be mentioned here that the continuation of 
the stream function into the supersonic region will rep- 
resent the desired physical solution only in the region 
AE’E in Fig. 5. It does not represent the analytic 
continuation of the desired solution for the flow down- 
stream of the limiting Mach wave—that is, in the re- 
gion between the limiting Mach wave AEF’ and the 
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positive y-axis. In this region, the representation gives 
a function that assumes, along the limiting Mach wave, 
the proper value of y and along the positive 7-axis, the 
value zero. Along the characteristics A’, discon- 
tinuities in the second and higher derivatives are per- 
mitted and will occur. This of course explains why the 
solution given in terms of the Chaplygin solutions does 
not give the multisheeted solution that must be ex- 
pected in the region between the characteristic A’ and 
the positive y-axis. 

The coefficients },, can be obtained in different ways; 
one can, for instance, minimize the integral over the 
square of the residues of Eq. (27), obtained with a 
finite number of coefficients 5,,.. This procedure has 
recently been used by Marschner® and proved satis- 
factory. We determined the coefficients 6, by a 
method similar to that of reference 1. Fig. 6 shows how 
well the condition along the characteristic is fulfilled by 
the coeflicients utilized in the present paper. Probably 
the method of reference 5 would give a somewhat better 
agreement, but basically the solutions can be con- 
sidered as satisfactory. 


The coefficients a,, were then determined in the same 
manner. These values are now to be inserted into Eq. 
(27). However, this does not give us the final answer. 
We must still determine the constant 8. For different 
values of 8, it is found that one obtains wedge profiles 
in the physical plane with different lengths of the sides 
when the solution due to the lift, Eq. (27), is super- 
imposed to that for the zero angle of attack. Since 
we are interested in the wedge with equal sides, we shall 
choose the constant 8 such that the solution due to the 
angle of attack does not change the location of the 
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Fic. 7. Change of the pressure distribution along the wedge 
profile due to the angle of attack. 
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corner of the wedge; the nose of the wedge will map to 
the origin of the physical plane in any case. Since the 
corner of the wedge maps to the characteristic EL’ of 
Fig. 5, this requirement is just that prescribed by the 
second condition, ¢, = 0, along this characteristic, 
The value of 8 was found to be 0.474. 

The continuation of the solution downstream of the 
limiting Mach wave is then computed by the method of 
characteristics. Here the result as given by Eq. (27 
is first transformed to the physical plane, and the solu- 
tion is then continued into the supersonic region. This 
result, when re-expressed in the hodograph variables 
will give the value of Y and ¢, downstream of the limit- 


ing Mach wave. 
(C) Calculation of the Pressure Distribution over the Wedge 


The solution of the boundary value problem given in 
the previous section will not give directly the change 
of the pressure distribution due to the angle of attack. 
One finds instead the amount by which a point with 
the conditions (n, 9) is shifted in the physical plane due 
to the angle of attack. Since the solution due to the 
angle of attack can be considered as a small perturba- 
tion of the basic flow over the wedge at zero angle of 
attack, the change of » may be expressed directly in 
terms of the hodograph solution due to the angle of 
For this purpose we shall introduce the 
The gradient of the Legendre po- 


attack alone. 
Legendre potential. 
tential in the hodograph plane gives the space vector 
in the physical plane, 


w(x +1)" "x=, 


wry = 


This equation holds true in the basic flow, as well as in 
the flow with the potential due to the angle of attack 
superimposed. We now want to find the changes 
6y and 60 of n and 6, which occur at a fixed point x, y 
by the superposition of the additional potential. In- 
serting the potential for zero angle of attack, ¢@, and 
the potential with angle of attack different from zero, 
@ + g, into the last equations, one obtains 


w*(k + 1)~ ‘x = &,(n, 0) 
wty = ®,(n, 0) 


and 


w*(x + 1)~''x = $,(n + bn, 0 + 60) + 
y(n + dn, 6 + 66) 
w*y = do(n + dn, 0 + 60) + go(n + Sn, 0 -+ 68) 


Linearizing the second pair of equations with respect to 
the superimposed potential g and with respect to 6, and 
66, one finds 


0 = ¢,,(n, 0)6n + $,6(n, 0)60 + ¢,(n, 8) 
Oo = $,0(N, 0)dn 7 Dea(n, 6, 660) + Ga(n, 6) 


Solving for 6n and 66 in these equations, we obtain 
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Fic. 8 
6, = —A (Deg &, i Pp ¥o) (31) 
do = —A~'(Pee¢e — P10) 
where 
A= D,,,° Pag D9 


These expressions give the change of @ and 7 to be ex- 
pected due to the angle of attack at a fixed point in the 
physical plane whose coordinates are given by the solu- 
tion to the flow for zero angle of attack. 

In order to find the change of the pressure distribu- 
tion due to the angle of attack, we must, according to 
Eq. (8), compute the value of 6 along the profile. 
Along the sides of the wedge, the expression for 6n given 
by Eq. (31) may first be simplified. We have the con- 
dition ¢, = 0 along the sides of the wedge which are 
lines of constant 6; thus, we have ¢,, = 0. Eq. (31) 
for én will, therefore, simplify to 


and the local pressure coefficient will then be 


(¢y/Pyn) (32) 





Pressure distribution due to thickness and to the angle of attack 


The total lift coefficient can be found from this equa- 
tion by means of an integration. 

Using the above formulas the pressure distribution 
and the total lift due to the angle of attack were com- 
In Fig. 7 the pressure distribution due to an 
this distribution is to be 


puted. 
angle of attack is shown; 
added to that for zero angle of attack on the upper side 
of the profile and subtracted from it on the lower side. 
Fig. S shows the change of the pressure distribution due 
to angle of attack of 1° for a wedge of 10 per cent thick- 
ness ratio as compared to the pressure distribution over 
the same wedge at zero angle of attack. The variation 
of the dC;/da at a = O (a in radians) with thickness 
ratio is illustrated in Fig. 9. The position of the center 
of pressure is independent of the thickness ratio and 
was found to be located at the 29.4 per cent chord 
position from the leading edge. In Fig. 10 the results 
of the present paper are shown in comparison with the 
values of the lift curve slope for supersonic Mach Num- 
bers for which the nose shock is attached. The latter 
results were computed by the shock-expansion method 
simplified by the transonic approximations. It ap- 
pears that the values at Mach Number | are consider- 
ably lower than what one would expect by comparison 
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Lift curve slope for various Mach Numbers 


to the supersonic results. Fig. 10 would then indicate 
that one must expect a considerable change of the lift 
curve slope through the range of Mach Numbers from 
1 to that for which the shock attaches to the nose. 
Moreover, according to reference 4, the slope of the lift 
curve at Mach Number | must even be horizontal; this 
is indicated in Fig. 10 by the short horizontal lines. 


This has been confirmed by recent investigations of 
Vincenti and Wagorier.® 


SECTION (IV)—FINAL REMARKS 


We have discussed in the present paper the flow pat- 
tern for slender profiles with a small angle of attack at 
Mach Number 1. It was found that the flow is uniquely 
determined by the boundary conditions and that no 
additional conditions such as the Kutta condition is re- 
quired. As a specific example, the flow over a wedge 
profile was computed. It was assumed for this case 
that the angle of attack was sufficiently small so that 
the essential unsymmetric effects are described by the 
linear terms of the angle of attack and that the effect 
of the local supersonic region arising at the nose could 
be neglected. 
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Appendix 
ESTIMATION OF THE INFLUENCE OF THE SUPERSONIC 
REGION AT THE NOSE OF THE WEDGE PROFILE 


The purpose of the present appendix is to show why 
for small angles of attack the effect on the lift of the 
supersonic region at the nose of the wedge can be neg- 
lected. There are several different methods by which 
this can be shown. 
examine the boundary value problem for the flow over 
Although 


A direct procedure would be to 


the nose with respect to the angle of attack. 
the boundary value problem cannot be formulated 
completely, one can still show in this way that the 
influence of the tip region on the lift force does not 
appear in the terms of first order of a development with 
respect to the angle of attack. 

We shall, however, proceed by another method that is 
simpler. It can be shown that the boundary value 
problem that was actually solved in the present paper 
corresponds to a flow over a wedge profile with a short 
flap that extends upstream from the nose (see Fig. 11). 
This flap becomes smaller and smaller as the angle of 
attack is decreased. The influence of the flow of the 
nose will then be estimated by examining the influence 
of this flap upon the overall lift. 

The boundary value problem that was computed is 
Here we have ¥ = 0 along the lines 
At the point A, which lies on the 


shown in Fig. 11. 
6 = 6 and —h. 
sonic line where 6 = a, a superposition of two prescribed 
singularities occurs which represents the free-stream 
singularity. The coefficients with which these singu- 
larities are multiplied are determined by the size and 
shape of the wedge. This boundary value problem 
has a unique solution in which no parameters are left 
arbitrary. When a = 0, a branch point of the stream 
lines will lie at the stagnation point. With an angle 
of attack, however, this can no longer be expected, and 
one will, in general, have a branch point of the stream 
lines that travels along the hodograph map of one of 
the sides of the wedge as the angle of attack changes. 
Which physical significance does such a hodograph 
The physical plane that belongs to this 
Here the zero stream line 


pattern have? 
flow is also shown in Fig. 11. 
forms a branch point not at the stagnation point but 
at a finite velocity. After the branch point, the two 
branches still proceed in the direction of —@, but along 
one of the stream lines the flow assumes higher velocities 
while along the other branch the velocity decreases. 
The flow near the branch point, therefore, resembles the 


flow over a flat plate inclined at an angle %. Along the 
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upper side of the flat plate the flow decelerates and 
finally reaches zero velocity at some point C. Here 
the zero stream line assumes the direction +%. We 
see, therefore, that the flow problem given by the hodo- 
graph formulation of Fig. 11 is the flow over a wedge 
profile at an angle of attack with the wedge having a 
short flap at the nose. The branch point of the stream 
lines at the tip of the flap. 
smaller as the angle of attack is decreased. 
fluence of this flap upon lift will be represented by the 
force that acts on the flap and furthermore by the 
changes of the pressure distribution over the wedge 
If the flap is small, 


The flap will become 
The in 


caused by the presence of the flap. 
then it is plausible that the force caused by the change 
of the pressure distribution over the wedge will be of 
the same order of magnitude as the force that acts 
directly upon the flap—that is, the ratio of the force 
on the flap and the “interference” force will be some 
finite number that does not tend to infinity as the flap 
becomes smaller. Therefore, one obtains an estimation 
for the entire influence of the flap if one determines 
only the force that acts on the flap. 

The flow velocities in the vicinity of the flap are ex- 
tremely small, and they become smaller as the angle 
of attack is decreased. Therefore, the force on the 
flap can be determined by an incompressible treatment. 
Since we have yy = O for 6 = +6 , the solution near the 
stagnation point can be represented by 

nO 


Y= da,-w-"” sin 


n=1 0 


[n — (1/2) ]( 3/6 1\ 70 
pS b,-w' Pee ‘cos [mn — 5 
n=1 A 


(A-1) 


where the constants a, and 6, are functions of the angle 
attack. If we examine the symmetry of the terms 
in the above equation, the form of the coefficients 
Eq. (A-1) can be rewritten 


a, and 6, can be given. 


as 


= 


y = > (Ano + A,1@" +...)4 


n 


‘ [10 — (1/2) ]( 9/be 
a > (bao + bya? +... )w' wre! eg: XK 
n 


1\ 70 
n — (A-2 
' 270 


, are constants independ- 


where d,0,@n1,---, andy, b,1,... 
For small values of a, only the region near the 


Under these conditions 


ent of a. 
stagnation point is of interest. 
the predominating terms in Eq. (A-2) are 


1/0 .-. 0 ./26 70 
y = Ajo°W’ -:sin + abyw’ COS 


The coefficients adj and by are numbers that are 


determined if one solves the boundary value prob- 


lem. 
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If one now introduces the complex velocity gq = 
w-e ", then Eq. (A-3) can be rewritten as y = JmQ, 
where 

. 26 
Q = —apg™’” + ta-dbyog" (A-4) 


The branch point of the stream lines is given by the 
conditions 


Oy/ow = 0; Oy /06 = 0 
which can also be written as 
dQ/dq = 0 


One finds that the branch point of the stream lines 
lies at 


260/ x| ‘ 27/00, +109 
g = a& "| dyo/ ayo) e 


Here the sign to be taken in the exponent is that of the 
sign of the expression },o/a,o. 

In the case of incompressible theory, the physical 
plane is related to the hodograph plane by the equation 


dz = dx + idy = dQ/q 


Therefore, we find from Eq. (A-4) that 
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z= vs (: = 4 oa all + 4 
Ao A 
l 1 jeer 
9 aby (= = ') q Tee | 


If one introduces the value of g for the branch point into 
this equation, then one obtains the length of the flap 
One finds 
if terms of higher 


as measured from the stagnation point. 
that it is proportional to a°®~°%” 
power of @ are neglected. 

If one now determines the velocity at the upper or 
lower surface of the flap for which z = 0, one finds that 
in the lowest order approximation the absolute value of 


206/m 


the velocity is proportional to @ According to 


Bernoulli's equation, the pressure difference between 
the stagnation pressure and the pressure corresponding 
to the velocity is proportional to w®. In the present 
case, the largest pressure difference occurring along the 


60/3 


. P 4 io ° 
flap is proportional to a We had previously seen 


(200/ x 


that the length of the flap is proportional to a” 
Therefore, the total normal force on the flap is found 


2+ (2600/7) 


to be proportional to a Consequently, the 
effect of the small flap at the nose of the wedge is of a 
smaller order of magnitude than the total lift of the 
wedge which is proportional to a. 
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Similarity Theory of Isotropic Turbulence 


J. C. ROTTA* 
Max-Planck-Institut fiir Strémungsforschung 


ABSTRACT 


The basic idea of this representation is that an arbitrary 
homogeneous and isotropic turbulent motion tends, in the 
course of time, toward a universal state, in which all statistical 
properties are determined by only three magnitudes. From 
dimensional arguments a general law of similarity for isotropic 


turbulence can be formed. The equation for the corresponding 


energy spectrum is derived from the well-known equation for 
the change of spectrum and has a one-parametric set of solutions 
But only one of these solutions is distinct, in that it satisfies the 
continuity condition in the range of the smallest wave numbers 
This distinct solution can be considered as corresponding to the 
universal state 

Comparisons of the theory with measurements of turbu 
lence produced by a uniform stream of fluid passing a mesh 
grid are discussed. According to this, it is at least probable that 
the universal state of homogeneous and isotropic turbulence is 
correctly described by the distinct solution of similarity 


(1) INTRODUCTION 


i pve TEMPORAL DEVELOPMENT OF a homogeneous 
isotropic field of turbulence can be followed by 
means of the equation for the propagation of the corre- 
lation first given by von Karman! or the equation for 
the change of spectrum, when suitable assumptions 
about the triple correlations or the transfer of energy 
between the wave numbers can be made. The deter- 
mination of the correlation function or the spectrum is 
a problem of continuity. That is to say, if the corre- 
lation function or the spectrum is known at any time, 
its further development can be determined theoretically 
up to infinity. Thus the difficulty lies, on one hand, in 
finding the correct expression for the energy transfer 
and, on the other, in fixing the initial distribution of 
energy. 

The approximate similarity of the correlation func- 
tions and spectra, already observed in the early days 
of turbulence research, led to the conclusion that the 
structure of isotropic turbulence approaches a shape in 
course of time which does not depend on the conditions 
under which the turbulence is produced. This obser- 
vation is again confirmed by recent measurements, 
reference 13. This state of turbulence could be con- 
sidered as a sort of statistical equilibrium. 

Under these circumstances it should be possible to 
determine the structure of the turbulence without 
having to consider its whole previous development, 
such that any vague assumption regarding the initial 
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distribution could be saved. Starting with such ideas 
von Karman first investigated the so-called self pre- 
serving correlation functions. Today it is realized that 
self preserving solutions do not exist except for the 
The assumption of self preserving struc- 
To overcome the contradic- 


final period. 
ture is too specialized. 
tions, a portion of wave numbers had to be excluded 
from consideration, and to explain the behavior several 
authors introduced further assumptions without any 
of them being satisfactory. 

The present paper deals with a more general treat- 
ment that does not start from the idea of self preserva- 
tion of the structure. It will be shown that this treat- 
ment leads to a unique solution that is valid for the 
whole range of wave numbers and Reynolds Numbers. 
2) GENERAL SIMILARITY OF ISOTROPIC TURBULENCE 


(2 


The basic idea is that an arbitrary homogeneous and 
isotropic turbulent motion tends, in course of time, 
toward a universal state, in which its statistical proper- 
ties do not depend on the special conditions governing 
the turbulence at the moment of its creation. This 
presumes a similarity of all isotropic fields of turbulence 
of the universe among themselves. It does not assert 
a similarity of the structure during the process of decay, 
as is assumed in the concept of self preservation of the 
i.e., the structure changes continuously dur- 
In this universal state, all statistical 


structure 
ing the decay. 
properties of isotropic turbulent motion in an incom- 
pressible Newtonian fluid are wholly determined by 
three independent magnitudes—-viz., a characteristic 
V(t), a characteristic length /(t), and the 
The characteristic 


velocity 
kinematic viscosity v of the fluid. 
velocity and length both are functions of the time /. 
From dimensional arguments we can formulate the 
general law of similarity for isotropic turbulence. The 
three-dimensional energy spectrum F(, ¢), which is a 
function of the wave number & and the time /, 1s given 


by 
F(k, t) = V>(U(Of(x, R) (2.1) 
where x is a nondimensional wave number 
x = ki(t) (2.2) 
R is a Reynolds Number 
R = V(t)l(t)/v (2.3) 
and f(x, R) is a universal function of x and Rk. This 


similarity includes as special cases all the similarity 
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hypotheses known up to the present time. In fact, 
only when such a universal state of turbulence exists 
does it seem to be reasonable to study the entire struc- 
ture without taking into account the initial conditions 
of the turbulence. The progress consists in that the 
universal function /(x, R) in Eq. (2.1) depends on two 
variables. The elimination of the second variable R 
reduces the consideration to the less general concep- 
tion of self preserving structure. 

Eq. (2.1) holds generally; however, V 
be chosen. It is convenient to put V proportional to the 
mean square root fluctuation u’ = 1/2 of any velocity 
The characteristic length / can, for ex- 


and / may 


component. 
ample, be put equal to the integral scale L or Taylor’s 
microscale A. It is also possible to introduce the energy 
dissipation ¢ per unit mass of fluid instead of one of 
these magnitudes. 

At this point we remember that locally isotropic 
turbulence is determined by two magnitudes only 
viz., the dissipation « and the viscosity v. From this 
the well-known law of similarity 


. 1 3 
F(k) = €’‘v “‘g(n) (2.4) 
follows, where ¢(7) is a universal function of its argu- 
ment 
n=kv'‘e “* (2.5) 
This particular case can also be deduced from Eqs. 
. r 1 1 3 
(2.1) and (2.2), putting V = e’‘vy’‘ and/ = e€ ‘‘y“* 


From Eq. (2.3) then results R = 1. Eq. (2.4), of course, 
holds only for a limited range of large wave numbers, 
The validity can be extended over the whole wave 
number range, when a third magnitude (say, wu’) is 
The function ¢ then depends again on 
yo? Goldstein" 


introduced. 
variables and ae “Sy 


two Viz., 7 
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Eq. (2.4), general- 
ized in this manner, expresses the same idea as Eq. 
(2.1). This consideration shows that there are different | 


also came to similar conclusions. 


ways to formulate the similarity, a nondimensiona] 
function of fwo independent variables is, however, 
necessary in each case to describe exactly the whol 
structure and to deduce the law of decay. 


(3) THe EQUATION OF GENERAL SIMILARITY 


In order to derive an equation for the corresponding 
similarity spectrum, it is necessary to transform the 
well-known equation for the change of the spectrum, 


Pa) *k *k 
y| F(R’, t) dk’ + I W(k’, t) dk’ + 
Ot/o of 0 


k 
q 2h *F(k’, t) ak’ = 0 (3,1) 
J0 





into an equation between the variables f(x, R), x, and 
R. Here the term II’(k,¢) of the energy transfer be 
tween the wave numbers can be written as 


Wk, t) = V%(t)w(x, R) (3.2) 


where w(x, R) is also a universal function of x and R. 
Since |’ and / are both functions of the time, the Reyn- 
olds Number R is also a function of time. However, 
we find, as a consequence of the supposed similarity, 


that the expression 


din l(t)/d1n V7(t) = —¢(R) (3.3) 


is a universal function of R only. This can be deduced 
from any two different integral moments of the equa- 


tion for the change of spectrum. 


Introducing the relations (2.1), (2.2), (2.3), (3.2), 


and (3.3) we obtain, from Eq. (3.1). 


ss fd | fray “(R +|) H(R) | R feax 

ash ae Oo) — ¢(R) fe 

RJo MI ORY 0? * ’ ™ 2 : J0 or * f | 
{ fd +2 ~ rR) | R ; f fd Uf soa / 4 . f° 2y/*fd ) =o Gs 
(Jo . 2 ' dRJo x J0 = RJo -_ x j 


Besides this equation, the nondimensional spectrum 
f must satisfy some additional conditions that depend 
on how V and / are defined. Although all the above- 
given statements and relationships are not influenced 
by the choice of V and /, it is convenient for purposes 
of the following analysis to define V and / so that 
simple forms result for these additional conditions. 


3 3 
V=5u? = f F(k, t) dk (3.5) 
— 0 


If we put 


7 : 1k : 
l=L=— F(k, t)‘ (3.6) 
4u'- Jo k 
the additional conditions to be satisfied are | 


In this case L is the transverse scale 


— mo F e af 
J fix, R) dx = 1; | f(x, Rr) “X = - (3.7) 
0 0 x oT 
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where g(7) is the transverse correlation function. The 
magnitude of ¢(R) of Eq. (3.3) then becomes 
a | 2fx dx + Rf w dx/x 
oT JO 0 ‘id 
((R) =1- = (3.9) 


8 ae 
2fx? dx 
0 


which can be derived from Eq. (3.1). 


(4) SIMILARITY SOLUTIONS FOR THE SPECTRUM 


From the great variety of solutions given by Eq. 
(3.4), those have to be selected which have no proper- 
ties depending on the arbitrary initial distribution. 

When R is extremely small, the energy transfer be- 
comes negligible, and Eq. (3.4) reduces to 


a *x 
I axitdy|_f fdx' - xs = 
0 J 0 
Bi * 
/ jax f 2x'*fdx’ (4.1) 
J 0 


ra x | 
| w al f fdx’ — &xf + (; 
0 0 - 


With the supposition that w represents an energy 
transfer from small to large wave numbers, these equa- 
tions have a set of special solutions that are also inde- 
pendent of the Reynolds Number and which are called 
by von Karman solutions of nonviscous sinilarity. 
For the portion of spectrum of x > x, the equilibrium 
spectrum of locally isotropic turbulence holds. The 
details of the solutions of nonviscous similarity depend, 
of course, on the relation between the energy transfer 
term |W(k,t) and the spectrum F(k,?¢). Batchelor® 
has shown that IW(k,/) commences with a 
power of k than F(R, ¢). 
nonviscous similarity have, for small values of x, the 


higher 
Accordingly, the solutions of 


form 
tix, ©) = const. x” (4.5) 
where 


m= (1/fF.) —1 (4.6) 


As already mentioned, from a spectral distribution, 
given at / = fy), the temporal development of the spec- 
trum can be deduced by means of Eq. (3.1). Thus we 
will find, in the case of decreasing Reynolds Number, 
all the appropriate solutions by a step-by-step integra- 
tion over the whole range of R if we start at a suffi- 
ciently high Reynolds Number with a solution of non- 


ae eet ‘andi’ 
a inf, fe -f w dx 
30 “ “xt 

w dx/x wd 
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In this equation R no longer appears. Hence, ¢(0) 
becomes a constant, and a solution is found which is 
independent of R, 
f(x, 0) = by ('/F — Fe — 4x? (4.2) 
The constants a and } depend on the definition of V 
and /. This special solution is identical with the self 
preserving solutions first suggested by von Karman 
and Howarth.' The magnitude of ¢(0) is determined 
by the preceding history of turbulence. 


At extremely high Reynolds Numbers, the influence 
of viscosity is restricted to large wave numbers. Con- 
sequently, a value x, can be defined such that the energy 
containing portions of the spectrum falls mainly within 
the range 0 < x < x, and the viscous dissipation is due 
essentially to that portion of the spectrum lying be- 
yond x. If V and / are defined by Eqs. (3.5) and 
(3.6), respectively, Eq. (3.4) for the range of 0 < x < x 
and Eq. (3.9) can now be transformed to 


(4.3) 


(4.4) 


viscous similarity, which is completed, for large wave 
numbers, by the equilibrium spectrum of locally iso- 
tropic turbulence. Since the first approximation to 
F(k,t) for small k 
throughout the decay, from a comparison of Eqs. (4.5) 


namely, F « k”—is constant 


and (4.2) it results that 


c(0) = ¢(@) (4.7) 


This simple relation associates the solution for ex- 
tremely small R with the corresponding solution of 
nonviscous similarity. 

The procedure described for obtaining the similarity 
A simplification of 
the equations is thus desirable. When V 
defined by Eqs. (3.5) and (3.6), respectively, the non- 
dimensional spectrum function f(x, R) changes fairly 
slowly with R. Further, the magnitude of ¢(R), for 
any value of R, does not greatly differ from the values 
corresponding to infinite and zero Reynolds Number. 
This follows from calculations of reference 2. Accord- 
ingly, the term 0f//OR may be neglected and ¢(R) may 


be made constant, 


solutions is somewhat laborious. 
and / are 


(4.8) 


Reintroducing the original functions F(k,¢), W(k, t), 
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and k in place of f(x, R) and x, we derive, from Eq. 





TABLE | 
(3.4), Variation of Characteristic Magnitudes with Time 
k ; 
2e(t) ee , | : Large Reynolds Small Reynolds 
er F(k , t) dk’ — kF(k, t) | — Number, Number 
3u'*(t)L Jo g -_ lor 
R— « R>0 
k *h 
f W(k', t) dk’ —v | 2k’F(k’, t) dk’ =0 (4.9) Energy per unit mass, 
0 , 0 (3/2)u’? CET STH . aig ts 
: : ‘ Integral scale, L x p2/(2+8) <7 
which seems to be a rather good approximation for the 9" “8"™' S*2" 
et ; ‘ fea ae F : ant oe 20) 
similarity solutions. The energy dissipation is Microscale, \ = (5 ) = ( ) 
‘eam Reynolds Number 
e(t) = | 2k? F(k, t) dk ae ee noes 
0 R = V3/2(u'L/v) c {2-O/ETO a (2-8/4 
The author first used Eq. (4.9) for an estimation of 
the spectrum of isotropic turbulence in 1948.4 Using ee i : ; 
S , the similarity condition of the final period when viscous 
the energy transfer relation, elas o oe , 
° dissipation prevails in the whole range of wave numbers 
rR - s “ke ° ° ° ‘ 
| W dk = f V F/p"”3 ax" f Db’2F dk’ (4.11) In Table | the variation of some, characteristic mag- Fr 
0 k 0 





nitudes with time in the large and small Reynolds 
suggested by Heisenberg,* solutions for different & have Number region are given. 
been investigated in reference 2. If desired, the solu- It may be explicitly pointed out that Eqs. (3.4) and ; 
tion of Eq. (4.9) can be iteratively improved by sub- (4.9) have not a unique similarity solution but a one- (5) 1 
stituting in the exact Eq. (3.4). Eq. (4.9) satisfies parametric set of solutions. The parameter & = 


§ T 
shart 
distt 
of d 


wav 


exactly, for large R, the condition of nonviscous simi- (1 + m) must be >1, otherwise the integral of the 


larity in the low and medium wave number range and — second condition Eq. (3.7) does not exist. In the next 
the equilibrium condition of locally isotropic turbu- section it will be seen which of these solutions is physi- | 
lence in the high wave number range and, for small R,  cally.significant. 
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(5) THE PHYSICALLY DISTINCT SOLUTION OF SIMILARITY 


The similarity solutions discussed represent that 
shape of spectrum toward which any initially given 
distribution F(k, fy) tends after a sufficiently long time 
of decay, when the initial distribution begins at small 
wave numbers with 


F(k, t)« R*~' (5.1) 


If we commence with any given distribution of F(R, to), 


for which the integrals 


f Pdb ond [era 
0 0 


exist, a step-by-step integration of Eq. (3.1), using the 
transfer relation Eq. (4.11), we can see how the simi- 
larity solution develops. The energy of large wave 
numbers decays much faster than that of small wave 
numbers and approaches rapidly to a local similarity 
according to the theory of locally isotropic turbulence. 
As a result, the wave number range, which contributes 
essentially to the kinetic energy, decreases continu- 
ously with time, and, finally, the remaining energy is 
distributed mainly over a wave number range that 
initially contributed only a small fraction to the total 
energy. When the range of wave numbers in which 
the energy is contained has diminished nearly to a 
range within which the spectrum did not initially de- 
viate essentially from Eq. (5.1), the spectrum has 
attained one of the above similarity solutions. The 
corresponding Reynolds Number depends on the value 
of the Reynolds Number in the initial stage. If the 
latter is not high, the procedure yields only similarity 
solutions having low R. If the initial value of R is, 
however, extremely high and the starting spectrum 
does not deviate too much from the similarity spectrum, 
the case of nonviscous similarity may be realized. A 
more detailed discussion on the physical background is 


presented in the Appendix. 


1 of turbulence behind a grid. Measurements are from refer- 


Theoretical slope, d In L/d In u’? = —(1/5). 


The preceding history of turbulence enters into the 
similarity solutions still in the form of spectrum at 
lowest wave numbers. Now we know from the con- 
tinuity condition at low wave numbers,’ as well as 
from Loitsianskij’s invariant relation, that the spec- 
trum of homogeneous turbulence begins with the uni- 
versal power law 


(5.2) 


F « k 


Thus only the similarity solution with the parameter 
£ = 5 is completely independent of any arbitrary con- 
dition of the initial state and can really be established 
for the whole wave number range. We conclude that 
this physically distinct solution corresponds to the 
universal state of isotropic turbulence. 

For extremely low Reynolds Numbers this solution 
is identical with that of the final period given by 
Batchelor and Townsend.’ The decay is u’? « ft” 
At extremely high Reynolds Numbers, the intensity de- 
cays according to the law u’? « ¢~'”’, which was first 
suggested by Kolmogoroff.* 

Fig. 1 represents the nondimensional spectrum for 
both cases, R = 0 and R = o, the latter being com- 
puted by the method given in reference 2. A more 
exact computation for the whole Reynolds Number 
range is being carried out and will be published later. 


(6) Decay oF TURBULENCE BEHIND A GRID 


To support theoretical deductions of isotropic turbu- 
lence, a great number of experimental investigations 
have been made in wind tunnels. Unfortunately, the 
case of turbulence produced by a flow through a mesh 
grid is not exactly homogeneous, though the deviations 
are, except for small distances x from the grating, not 
serious. Nevertheless, it is possible that some of the 
observations made are not causally concerned with the 
idealized case of a homogeneous isotropic field of tur- 


bulence. 
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This seems to come true for the decay of energy. 
It is seen that existing experimental results, in the 
initial period, show a law of decay, u’? « x~', which 
does not agree with the similarity solution with £ = 5. 
The experimental law of decay corresponds rather to a 
solution of Eq. (4.9) with & = 2. This solution has 
been suggested by Heisenberg® and is mathematically 
distinguished by a constancy of the Reynolds Number 
throughout the decay. At small wave numbers this 
theory predicts a linear relation between F(k, /) and k, 
which is, however, at variance with the continuity 
Any experimental evidence for the linear 
Since the 


condition. 
part of the spectrum has not been found. 
similarity solutions represent the remainders of any 
initial distributions that have, at small wave numbers, 
the form of Eq. (5.1), a plausible explanation for 
setting-in of the solution with £ = 2 can hardly be 
given. If this consideration is correct, it is not prob- 
able that the actual behavior of isotropic turbulence 
is exactly described by this solution. 


On the other hand, the measurements point to the 
fact that, in the early stage of turbulence behind a 
grid, the experimentally observed decay u’” « x7! 
generally holds. It is not too much to suppose that a 
characteristic phenomenon of turbulence behind a grid 
is here in question. This supposition is supported by 
the fact that the turbulent energy in the wake of a 
single rod decays also approximately by the same law. 
Behind a grid, each rod forms a wake in a similar way. 
Possibly the mechanism that governs the law u’? « x7! 
is already set in motion before the mixing of the wakes 
begins. Undoubtedly the phenomena in this stage of 
turbulence are not yet quite understood. But it is 
probable that they cannot be explained with the con- 


ceptions of homogenity and isotropy only. 


(7) COMPARISONS WITH EXPERIMENTS 


Comparing experimental observations, we must bear 
in mind that the turbulence needs a certain time until 
similarity according to Eq. (2.1) is achieved. Hence, 
only such observations that have been made at 
not too small distances from the grid can properly be 
considered. It is for this reason that the discrepancy 
with the decay in the initial period does not matter 
greatly. On the contrary, for extremely low Reynolds 
Numbers and extremely large distances, the theory is 
the same as that suggested by Batchelor and Town- 
send,’ which was found to be in full accordance with 
experimental observations. Therefore, what we need 
are comparisons at higher Reynolds Numbers and 
sufficiently large distances. 

From published experiments,’ it results that the 
decay law shows, at certain distances x (x/\/ ~ 100 
to 200), deviations from the initial decay u’? « x7! 
which seem to correspond to the law predicted by our 
This is only a qualitative comparison. An 


. 9 —10/7 : 
with the law u’ «x ‘, which 


theory. 


extended region 
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corresponds to the case with & = 5 at high Reynolds 
Numbers, cannot be observed because the Reynolds 
numbers of measurements are too small. But we are 
able to check quantitatively the relation (3.3) between 
L and u’*. It was pointed out that, for extremely 
high Reynolds Numbers, as well as for extremely low 


Further, it 


ones, the value of ¢(R) equals 1/& = 1/5. 
was mentioned that, for moderate Reynolds Numbers 
also, no great deviations of ¢(R) from 1/5 are expected. 
On Fig. 2, the longitudinal scale 2Z divided by the 
periodic spacing M/ of the turbulence-producing grid is 
plotted, to a logarithmic scale, against the intensity u” 
divided by the square of thesuniform stream velocity 
U relative to the grid. Only~such measurements of 
references 7 and 9 have been used, which were carried 
out at distances x > 100/17. 
It is seen that 


The straight lines drawn 


in have a slope of —1/5. 
d1n (L/d) Inu’? ~ —(1/5) 


when x/.\/ is not too small. This comparison com- 
pletely confirms the prediction of the theory and 


covers a range of Reynolds Numbers up to 
N/v ~ 380(7 3/2u'L/v ~ 45 
uA/v IY 3/2U L/t ») 


whereas the final period sets in, according to Batchelor 
and Townsend,’ at u’A/v < 5. 

A direct examination of the energy decay has been 
carried out over a wide range of distances and various 
grid shapes by Baines and Peterson.'? The compared 
measurements were made at high Reynolds Numbers of 
the order u’L/v = 1,000 to 1,600. These authors 
state that their data, for large values of x, fit the decay 
much better than the law vu” « x7'. 
10/7 to be in 


ro .—10/7 
law uu’? « x 


Also, Frenkiel'! found the decay u’? « x 
accordance with measurements of several authors. 
The most critical check of the decay is a comparison 
of measured values of the microscale \. Values of \° 
from reference 13, which have been determined by 
different methods—viz., Townsend’s differentiation 
method, Liepmann’s method of zero countings, and 
calculation from measured spectrums—are presented 
in Figs. 3 and 4. For distances x > 100 em. (x/1/ > 
60), the results of each method follow closely the slope 
d\?/dx = 7v/U. This is the slope predicted by the 
present theory with & = 5 at high Reynolds Numbers. 
These experiments agree with the initial decay u’” 
x! for a much less extended region than those of 
13 also found 


reference 9. The authors of reference 


that the quantity 


‘F (0) U 


C = (u’?/ug'?) 


does not show a definite trend. This is identical with 


din L/d1nu’? = —(1/5) 


since the frequency spectrum F,,(0) U at zero frequency 
is proportional to L. 

To check the theoretical spectrum functions, meas- 
urements made at sufficiently large distances from the 
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grid are scarcely available, except for those at extremely 
low Reynolds Numbers. However, the decay, espe- 
cially the slope of the curve (A*, xv/u), which corre- 
sponds to the second derivative of u’* with respect to 
x, will be perceptibly influenced already by such local 
deviations of the spectrum from the universal state 
which lie within the scatter of measurements. 
comparison with measurements of spectrum made at 


Thus, a 


smaller distances, where the decay law of the universal 
state has not yet been completely established, seems to 
be justified. 

The exact shapes of the theoretical spectrum and 
correlation functions are, of course, influenced by the 
particular assumptions regarding the energy transfer 
_term. If we decide to use Heisenberg’s relation (4.11), 
the similarity solutions are still affected by the choice 
of the constant x. For that reason Figs. 5 to 9 com- 
pare the solution of Eq. (4.9) for §& = 5 and R= © 
with measurements, the comparison then being inde- 
pendent of the arbitrary choice of a constant. As an 
alternative, the solution of Heisenberg’s similarity 
hypothesis ( = 2, R = @) is also plotted. 


The theoretical frequency spectrum* for — = 5 fits 
well the measurements of Simmons and Salter’? in the 
low-frequency range, Fig. 5. The theoretical curve 
plotted in a normalized form on Fig. 6 also agrees well 
with the measurements of reference 13 made at ex- 
tremely high Reynolds Numbers. The deviation of 
the theoretical transverse correlation function, Figs. 
7 to 9, from measurements by Batchelor and Townsend? 
are mainly caused by the final Reynolds Number. 
To illustrate this, the solution for R = 0 is also repre- 
sented on Fig. 7. The measurements join the theo- 
retical curve with € = 5, R = © more closely when the 
mesh grid Reynolds number UA//v is high. The 
agreement of the solution with = 2 and of the func- 
tion ¢ = (1 + ¢*)~' on Fig. 6, with the measurements 
is worse. 


From these comparisons of the theory with measure- 
ments of different origins, it can be concluded that the 
turbulence behind a grid obeys, with increasing dis- 
tance, more and more the laws that correspond to the 
discussed similarity solution with § = 5. Therefore, 
it is probable that this solution correctly represents 
the universal state of isotropic turbulence. 

* The frequency spectrum is obtained from the three-dimen- 
sional wave number spectrum by 


= n 2r es k? — (22n/U’)? 
UF, i —— F(k) dk 
l uJ 2xn/l ke 


where is the number of cycles per second. 
+ The transverse correlation function is related with the non 
dimensional spectrum function by the equation 


r SEF if fg. R) ( r xr 
y{—,R) = si 
g (: ) A | 9 x L in L + 


xr x) 
cos _ sin dx 
L xr £ 
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Appendix 


ON THE ESTABLISHMENT OF THE UNIVERSAL STATE 


Batchelor and Townsend’ have clearly shown that 
isotropic turbulence approaches to a universal state in 
course of time. This consideration neglects, however, 
the nonlinear terms in the Navier-Stokes equations 
and thus proves only that the universal state is estab- 
lished at extremely low Reynolds Numbers—i.e., in 
the final period of decay. We must examine the ques- 
tion whether there is a physical basis for the assump- 
tion that a universal state may be accomplished also 
in the case of any arbitrary Reynolds Number. 

Some study referring to the establishment of simi- 
larity has already been presented in reference 2. It 
could be shown that the large wave number part of 
spectrum decays considerably faster than the small 
wave number part. The direct consequence—viz., the 
continuous diminution of the wave number range con- 
taining the major part of kinetic energy and viscous 
dissipation—has been experimentally verified (com- 
pare Figs. 17, 18, 20, 21, 23, and 24 of reference 13). 
Further, using Heisenberg’s transfer relation, Eq. 
(4.11), it has been pointed out that local disturbances 
of the spectrum decrease more rapidly than the un- 
disturbed Disturbances at large wave 


numbers decay much faster than those at small ones. 


spectrum. 
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This means that there is a strong tendency of the 
spectrum to compensate every irregularity and dis- 
continuity. 

Now we are going to consider two spectra that coin- 
cide initially within the range 0<k < k, and differ 
from each other for k > k,. For simplicity we assume 
firstly, as has been suggested by Kovasznay," that the 
contribution to the energy transfer 


- 
J W(k’) dk’ 
0 


comes from a narrow wave number band in the imme- 
i.e., that W(k) depends only on 
In this case we see, 


diate vicinity of k 
F(k) and k in the vicinity of k. 
from the equation for the change of spectrum, 


(OF/dt) + W + v2k2F = 0 (1) 


that the differences of F for k > k, do not influence the 
change of F inthe range0 << k <k,. Consequently, the 
two curves of spectrum will still coincide in the range 
0<k <k, after an arbitrary period of time. 

The assumption that only an extremely small wave 
number band in the vicinity of k contributes to 


+) 
I W(k’) dk’ 
0 


is not exactly true. But one may say more generally 


that the essential contribution to 


nb 
/ W(k’) dk’ 
0 


comes, in fact, from a wave number band of restricted 
finite width in the neighborhood of k. This is in ac- 
cordance with Heisenberg’s transfer relation Eq. (4.11) 
and is the physical hypothesis underlying the theory 
of locally isotropic turbulence. Thus the change of 
F(k) within the range 0 <k < k, is indeed not com- 
pletely unaffected by the range k > k,, it will certainly 
not be affected, however, by the spectrum of wave 


x 


numbers much larger than ,. 

Next we add a consideration that applies to the 
initial stage. The energy contained in the region 
0<k<k, in which the spectra coincide may be a 
small fraction of the total energy only. In this case, 
according to Batchelor,’ it is possible to express the 
transfer term approximately in the form 


W(k) = n2k?F(k) (2) 


where 7 is the turbulent friction coefficient that is 
considered to be nearly independent of k within the 
range 0<k < k,. 
time and changes according to the decay of the part 
k>k,. Substitution of Eq. (2) into Eq. (1) and inte- 
gration gives 


n is, of course, a function of the 


7 


F(k, t) = F(k, i) exp| —? 2 | (n + v) dt’ (3) 
lo 
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which holds for the region 0 << k < k,. Consider again 
the two spectra Fi, fF, having different n (because of 
different shapes at higher wave numbers k > &.), 
From Eq. (3) it follows that the curve F(R) has the 
same shape for 0<k < k, in both cases, if different 
time intervals of decay are compared, such that 


"hh a 
| (m + v) dt’ = | (m2 + v) dt’ (4) 
to lo 


We can now discuss more thoroughly the behavior 
of the two spectra having properties defined above. 
During the first period, according to Eq. (3), no differ- 
ences occur between the shape of the two spectra in 
the range 0<k< k,, initially significant 
differences in the upper range k > k,. With proceeding 
time, the energy of the part k >k, decays and the 
range of k for which Eq. (2) is valid decreases. As 
an important fact, however, the initial difference 


despite 


between the two spectra decreases rapidly at the same 
time and a large part of spectrum with higher wave 
numbers quickly approaches a state of local similarity. 
Consequently, although the two spectra will no more 
exactly coincide for the entire range 0 < k < k, after 
a certain time of decay, the initially existing difference 
will not be able to influence the shape within a less 
extended finite region 0<k < kj <k,. When the 
energy initially contained in the range k > k; has died 
away, the remaining energy is entirely contained 
within 0 < k < k;, the shape of the spectra then being 
wholly determined by the initial conditions for 0 < 
Ss &,. 

From this consideration different 
spectra having, initially, the same form in the wave 
number range 0<k <k,; for instance, Fo k*' 
achieves with certainty equal forms over the entire 


result several 


energy and dissipation bearing wave number range 
after corresponding lengths of time. This is the same 
conclusion that was obtained above in Section (5). 
In case the initial shape, within 0 < k < k,, is of uni- 
versal type—i.e., F « k*—the spectra take a universal 
form. This is the physical picture underlying the idea 
that any homogeneous isotropic field of turbulence 
approaches a universal state with time. 

Any statement regarding the extension of the region 
O0<k <k; and the energy contained at the moment 
when the universal state is accomplished cannot be 
made. They depend on the initial conditions. How- 
ever, there is no limitation to the Reynolds Number. 
There is some indication that the given criterion is too 
strong. Because of the tendency of the spectrum to 
compensate irregularities and to approach to a quasi- 
equilibrium in the upper wave number range, the 
universal state actually is expected to be established at 
least when the maximum of the spectrum curve ap- 
proaches the limit within which the initial deviations 
from the form F « k* are small. 


(Concluded on page SOO 
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LUDWIG PRANDTL 
1875-1953 





and beloved person 


Gottingen Institute. 





The long-time director of the Kaiser Wilhelm and Max Planck Institute 
for Flow Research, Prof. Dr. phil., Dr.-Ing. e.h., Dr. rer. techn. h.c., D.Sc 
h.c. Ludwig Prandtl, died in Géttingen on August 15 in his 79th year. In 
deep sorrow we mourn the passing of this prominent scientist, this gracious 
As founder of flow research he has labored for this 


science all his life and has created the reputation and tradition of the 


Max Planck Society for Scientific Research 


Otto Hahn, President 








AS THE FIRST HALF-CENTURY of human flight nears 
its close, this official announcement informs the 
scientific world of the death of the architect and builder 
of modern fluid mechanics. His creative mind dis- 
cerned the essential key physical phenomena in complex 
fluid motions, established new concepts that are the 
daily tools of thought of aeronautical scientists and de- 
signers and formulated the analytical and quantitative 
procedures for mathematical study and optimization of 
proposed designs of aircraft. His influence on aeronau- 
tical progress within the last 30 years has been greater 
than that of any other individual except the Wright 


brothers themselves. 


Before Prandtl, aerodynamic theories were essentially 
mathematical exercises completely sterile so far as any 
application to practical problems was concerned. Ex- 
perimental aerodynamics was purely empirical. When 
Prandtl, after a short but successful service as me- 
chanical engineer, with Maschinenfabrik Augsburg- 
Niirnberg and 3 years as professor of mechanics at 
Hanover, accepted appointment to the chair of applied 
mechanics at Géttingen, he found a long-established 
tradition in process of further development by the great 
mathematician, Felix Klein. Klein believed that there 
should be a close connection between the general ab- 
stract mathematical theories and the experimental facts 
and practical applications. Prandtl proceeded to select 
suitable basic problems, to create new experimental ap- 
paratus, and to refine experimental methods, all looking 
toward a method of analysis embodying the essential 
physical aspects of the problem and in agreement with 
the data resulting from experiments. He drew around 
him a group of collaborators who were stimulated by 
his leadership to productive and creative extensions and 


applications of his own basic concepts. 


Prandtl’s most important contribution to aeronautics 
was his wing theory-—i.e., the discovery and computa- 
tion of induced drag by means of which rational design 
procedures were established to compute the influence of 
wing span, to compare monoplanes and biplanes, to es- 
timate influence of stagger and gap in biplanes, and to 
compute load distribution, down-wash, and so on. The 
development of the high aspect ratio monoplane wing 
of great efficiency was directly stimulated by Prandtl’s 
work. This theory was developed during the first World 
War and not known to others outside Germany until 
1921. 

Another concept regarded by many as equally im- 
portant is that of the boundary layer, historically ante- 
dating the wing theory. More correctly, Prandtl’s con- 
tribution was a theory of the motion of fluids with small 
friction, dealing with the whole flow field. Approximate 
methods of solution of the differential equations were 
given, and the large effect of a small friction in produc- 
ing separation of flow in adverse pressure gradients was 
explained. 

In these days of supersonic interest, we must rank 
highly Prandtl’s contributions to compressible flow 
first, his studies of the efflux of gases and vapors under 
high pressure. Some doubted that supersonic efflux 
speeds were possible, but Prandtl applied the Riemann 
theory of shock waves and the experimental schlieren 
method of Mach to give a full description and explana- 
These contributions were made in the 1904-1907 
period. Later, as propeller tips approached supersonic 
speed, he developed a linearized theory of compressible 
flow. Such terms as Prandtl-Meyer expansion and 
Prandtl-Glauert approximation remind us of his work in 


tion. 


this field. 


(Concluded on page 800) 
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On Laminar Boundary Layer over a Rotating 
Blade* 


H.S. Tan 


Graduate of Aeronautical Engineering, Cornell University, Ithaca, 


July 20, 1953 


" Begnanen A CYLINDRICAL BLADE extending infinitely in y- 
direction and rotating about z-axis with angular speed &. 
The boundary-layer flow about the blade is given by differential 
equations, 
uu, + vuy + wu, — 2 — x0? = —(1/p)pz + vu 
uv, + vy + we, + 2uQ — yQ? = —(1/p)py + vz (1) 
Uz + vy +w, = 0 
together with boundary conditions 
u=v=w=Q0, (2) 
u>U, v>V, worW, s+ \ e 
where U, V, IW are velocity components of outside potential flow, 
readily obtainable as 
U = yO 4,.(x, 2) 
7 = Q[d:(x, 2) — 2x] > (3) 
V = yOo,.(x, 2) 


_— 


in which ¢; denotes the potential for plane flow past the cylinder 
in uniform stream along x-direction.! From Eqs. (3), the pres- 
sure gradients pz, p, follow at once by using Bernoulli’s equation, 
(1/p)pz = x2? — UU. — VV: | 
ss > . (4) 

(1/p)py = yo? — U0*, \ 


Now, because of the nonlinearity of equations and coupling of 
dependent variables, a general solution of Eqs. (1) cannot be 
expected. However, if one restricts oneself to regions in the 
neighborhood of leading edge, where x/y < 1, it will be shown 
that a perturbation procedure can be set up so that all the suc- 
cessive approximations are linearized. Indeed, both u and v 
velocities are then expressible in terms of a single parameter 


n = (2/2) V yQ/xvas follows: 


- © 3 
yo > x\" f (9) 
v= 2 , y Ge (n) | 


To proceed, let us first observe that since Eqs. (3) give U,W~ 
y, V ~ x, we are justified in writing 

* This research was supported by the USAF under Contract No. AF 33 
(038)-21406, monitored by the Office of Scientific Research, Air Research 
and Development Command 


eu e+e +a +... 
vy =v) + yf) + 

3 , ) (6) 
w= w w)) + wi?) + 


peer e” + 2% + . 
where ( )" ~ O(x/y)". 
On substituting Eqs. (6) into Eqs. (1), collecting terms accord- 
ing to different orders in (x/y), and equating to zero respectively, 
it is seen that the mth approximation are given by differential 


equations, 
n n—1 \ 
do [u *)y,(n—8) wey r—)| + > vy, (n—-1-s) — | 
0 1 
2v"2—-D OQ — x75." = —(1 p)Pe\™ + viz," 
n n—1 ( 
7) 

lu yi (n—a) 4 ayla) ysin—a)] 4 Z yay, (ma) 4 
0 1 

2u@—D Q — yQ25," = —(1/p)py*—) + vse" 

u,\” W, n +. ty m1) = () 


together with boundary conditions 





am = = w™ = (), z=zQ@ | : 
- =: r, (8) 
ur —> UY), ve VO, wos W™, ss me | re 
SOLUTION FOR ROTATING FLAT PLATE 
In case of a rotating flat plate, Eqs. (3) and (4) give 
C=, Ve=-x, We=0 ‘ 
(9) 
pb: = py = 0 
so that Eq. (7) leads readily to the following results: 
(1) Zero-Order Solution: 
wuwx + wud = vu, | (10 
uz + w) = 0 \ 
with vw = w® = Oats = 0, and W > yQ, wW—>O0aszs— 


Now Eggs. (10) are exactly the well-known Blasius problem 
On introducing 


Z yQ 2 
= 3 0 = , 
’ 2 \ xv . 4 , 
(11) 
l yQv _ : 
w= (nF’ — F) 
2 * 
where F is a function of 7 alone, we arrive at equation 
F’'"' + FF’ =Q@ (12 


with F(O) = FO) = O and F’( 2) = 2. The solution of Eq 
(12) was obtained by numerical integration.? * 
(ii) First-Order Solution: 


wus) + uy, + wu, + wud = vu, } 
uv,) + wv) + 209) — yO? = w,,"! > (13 


uz™ + uw, = 0 : \ 


iW, 
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with wo) =v) = ww” =OQOats=0, and «a = ww) = 0, v = 
—«402. asz— 

By comparing Eqs. (13.1), (13.3) with Eqs. (10), it is easily 
verified that wu? = wi 0 are solutions for Eqs. (13). 

Now as Eq. (13.2) contains only v, it can be solved sepa- 
rately. Indeed, putting 


(14) 


v) = (x2/2)G'(n) 


Eq. (13.2) reduces directly to following differential equation for G 


G’"! 4 FG" — 2F'G’ = 8 F’ — 1) (15) 
with G(O) = G’(0O) = 0, and G’(“) = —2. Eg. 15 has been 
integrated numerically for G’ in references 4 and 5. 
(iit) Second-Order Solution: 
wus’) + uu, + vu? + wu?) + 
w'2)4.9 — QyVQ — x0? = pu," . 
: » P (16) 
uy, + wh = pvz,'? 
Uz*’ + We +v,) =0 
with vw? = vw? = w’’ = Oats = Oandass— 


By comparing Eq. (16.2) with Eq. (138.2), it is easily verified 


that vo‘? 0 is a solution of Eq. (16.2). Now, 
nu? = (x7Q/2y f’(n) 
xQ vx | 1 a (17) 
w?) = 2 (nf' —f)— 1 1(nG —G)-—f | 
y yQ 14 
satisfies Eqs. (16.3). Introducing Eqs. (11), (14), and (17) into 


Eqs. (16), we arrive at the differential equation for f as follows: 


+ Ff’ — 4F’f’ + 5F"f = 2(F’ — 4)G’ + F"G —8 (18) 
with f(0) = f’(0) = O and f’(o) =0 
(iv) Third-Order Solution: 
uu, + uu,9 + wn.) + wu” = vu,, 
u'r + By) + yDy DY) + gly (3) + (19) 
wo) + By] = ypu,, 
u®) + wv," = 0 J 
with uu‘? = 7 = uw) = (atz = Oandaszs— 


Comparing Eqs. (19.1) and (19.3) with Eqs. (13), it is easily 
verified that 1 = y'3 0 are solutions. 
Indeed, put- 


Now Eq. (19.2) can be solved for v‘*’ separately. 


ting 
(20) 


Eq. (19.2) reduces directly to the following differential equation 
for g: 


g’’’ + Fg" — 6F'g’ = 2G’ + 4)f' — 5G"f + GG” (21) 


with g(0) = g’(0) = Oand g’(~) = O. 


If we denote functions F, G, f, g, etc., by fo, fi, fo, fs, etc., respec- 
tively, then Eqs. (12), (15), (18), (21), etc., take on the following 


form: 

fr’ + fofo’’ = 0 

— = Sofi’ = 2fo'fr’ = fo’ — 1) 

Ja!" + fofe’’ — Afo'fe’ + 5fo' "fe = fo’ — 4)fi’ + fof: — 8) (22) 
fi'"” + fofs’’ — 6fo'fs’ = 21’ + 4)fe’ — of fe + Afi 


with f,(0) = f,’(0) = 0 for all m and with fo'( ©) = 2, f,'(o) = 


—2,andf,’(©) = Oform > 2. Weare thus able to write 


yQ x \2" 
nwt > (*) fan!(n) 
2-0 \y 
yvQ x\2ar! 
= - > Son 
> Oo \y 


It might be of interest to compare Eqs. (22) with corresponding 


(23) 


t 


+1'(7) 


set of differential equations for two-dimensional boundary-layer 
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flow expanded in powers of x, having a linear outside velocity dis- 
tribution U = U, + cx. Indeed this noted resemblance makes 
it possible to carry out our numerical calculation entirely after 


the manner of reference 7. 
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On the Evaluation of the Function /\(\/, o) 


G. Zartarian and H. M. Voss 
Aeroelastic and Structures Research Laboratory 
Institute of Technology, Cambridge, Mass 


July 27, 1953 


Massachusetts 


— THE FIELD of unsteady supersonic aero- 
dynamics, it becomes necessary to evaluate a function of the 
form 

+] ; 


° w \ 
exp. (—iau)Jo uj} u* du (1) 
0 M 


wherein M is the Mach Number, @ is a frequency parameter de 


fy (M, @) = 


x 


fined by 
@ = 2kM?/(M? 1) (2) 


and k is the familiar reduced frequency parameter in unsteady 
flow. For example, the forces and moments for the two-dimen- 
sional oscillating airfoil in supersonic flow are given by Garrick 
and Rubinow! in terms of functions that can be reduced to com- 
binations of f, with A = 0, 1, 2, and 8. The work by Miles? 
on the forces and moments on a rigid oscillating wide delta in 
supersonic flow requires, in addition, the function /4. 

Further work with the wide delta plan form shows that the 
elastic motions will require the values of f, for > of even 
For example, the moment due to parabolic cam- 


the camber 


higher order. 
ber will involve f;, and the generalized force in 
mode due to the camber mode itself requires evaluation of fe 
Similarly, spanwise modes (normal to the stream) lead to integrals 
of the form of Eq. (1) but with Bessel functions of higher order. 
Such integrals can be reduced by partial integration to the exact 
form of Eq. (1) but result in higher order \. The generalized 
force in a parabolic bending mode due to the mode includes 
fy up toand including fs. 

Previous methods of evaluation of the functions f/f, have 
relied on precise determination of fy and continued use of the 


recursion formula developed by Borbely,® 


M2 — | =f, (M. o) j 4 ( 1 ’ l ] j ; ( - ) 
) y @w) =< A) a Aw 
M2 lie : \' ; : | ; ‘ M 


I — iw ; l ‘ - l 
M e @1,(*) + 21 2d) fy, (AL, @) + 


1 : 
(1 — r)? - f _9(M, @) (3) 


wW 


A 


The function fo was evaluated by Borbely to five decimal places 


from the infinite Bessel function series, 
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— {M?—1_\* 
fo (M, &) = e~* > ( s) x 
n=0 M? 


— (ins inti 
Sal (ae 41) ONO) tarde) (4) 


Also included in reference 3 are fi, fz, and f; as calculated from 
the recursion Eq. (3) above, and these results are reproduced 
by Miles.?, Schwarz has determined fy to eight significant figures 
by expansion of the function as an infinite integral series that 
leads to a doubly infinite power series expansion. The results of 
Schwarz are continued in reference 1. 

A difficulty involved in use of the above is that the arguments 
M and @ for which the functions are tabulated are not con- 
venient to the problem at hand, particularly for the case of a 
control surface wherein modified f, functions are required. 

gx(M, @) = fy (MM, ox) 

hy(M, @) = fx(M, atl — x1}) 
where x; is the location of the aileron hinge line.!. Furthermore, 
it is found that, because of the inherent nature of the recursion 
formula, evaluation of f, for higher order \ by Eq. (2) results 
in continued loss in accuracy. 

For these reasons it is desirable to have available a method of 
evaluation of fy, which is accurate in itself to any degree and 
which permits simple variation of the frequency parameter, @. 
Such a form can be obtained by expansion of both the exponential 
and the zero-order Bessel function of the integrand as infinite 
power series, followed by term-by-term integration. There re- 


sults, 


~ (—1)*( —2)’ 1 \* 
A(Ma=-> >d \ ( ) x 


s=Or=0Or+t2s+At]1 


(a) + (5) 


This can be more simply expressed in terms of the real and 
imaginary parts of the function, 


— 1)” 
( a2,(M) (@)2" — 


fy, (M, @) = 
a ( “ n=O’ + 2n + 1 


7 (—1)" 
1 bon( M) (@)?" (6) 
free ree - a 


where the coefficients are themselves finite series, 


n 
I Ry 
on( M) = 
a2,(M) p> (j!)? [2(n — 7)]! (3) 
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Since both the exponential and the Bessel functions are analytic 
in the entire finite complex plane, their power series expansions, 
and hence the integrand of Eq. (1), is uniformly convergent 
throughout the finite complex @: plane, and in particular along 
the real axis. Then, by change of variable to au, the integral js 
everywhere regular, including the origin, and convergence of the 
term-by-term integration to the value of fy, (7, @) is assured 
for all values of 2 < © and M from the theory of functions (cf 
reference 5). 

The advantage of the series representation above is that the 
dependence of f, on A, JA/, and @ is completely separated, 
so that the coefficients a2, and bs, can be tabulated for a given 
Mach Number. 
\ and @ to an accuracy limited only by the tabulated coefficients 
It will also be noted that, for a given accuracy, the number of 


The function can then be determined for any 


terms required decreases with increasing \ and J/. 
Table 1 presents values of the coefficients a2, and b», for the 
range » = () to 6 and for most of the customarily used Mach 
Numbers. 
and the results were compared with the values obtained from Eq 


(4). In all cases the results agree to within one in the ninth 


This table has been used to calculate f, to nine places, 


place. Since the accuracy improves with increasing A, the ac- 
curacy for higher A is assured. 

In conclusion it is worth noting that one set of coefficients can 
be written in terms of the other—i.e., 


n 


nj Qin + (Qin + — 1) 
i * 2, el ~9) ars (0 
j=0 [2(n —j + 1)!! 


1/6, B; = 1/30, 


7/6, 


where the &,’s are Bernoulli’s Numbers, B; = 

B; = 1/42, B; = 1/30, By = 5/66, By = 691/2,730, By; = 
This form is often more convenient for calculation of the 

second set of coefficients when the first set has been evaluated 
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(7) 
n 1 1 1 2i 5 Knopf, K., Theory of Functions, Elements of the General Theory of Ana- 
bo,( 11) a yi ~ ( lytic Functions, Translated by Bagemihl, F.; Dover Publications, New York 
j=0 (j!)? (2m — 7) +1]! \2M 1945 
TABLE 1 
Table of Functions a2, be, 
M 5/4 10/7 3/2 5/3 2 5/2 
o=—> (50/9)k (200/51)k (18/5)k 25/8)k (8/3)k (50/21)k 

ao l l 1 l 1 1 

ay 0.660 ,000 , O00 0.622, 500,000 @:611., 311,111 0.590 ,000,,000 0.562 ,500 ,000 0.540 ,000 000 
as 0.128 ,066 ,667 0. 106,668 , 230 0. 100,308 ,642 0.088 ,691 ,667 0.073 ,893 , 229 0.062 ,066 , 667 
a6 0.011 ,369,333 0.008 ,419 , 900 0.007 ,599 , 832 0.006, 171,639 0.004 ,488,119 0.008 , 257 , 333 
as 0.000,571,718 0.000 ,377,178 0.000 ,327 ,040 0.000 , 244,416 0.000, 155,716 0.000 ,097 ,918 
a0 0.000 ,018 , 450 0.000 ,O10, 850 0.000 ,009 ,040 0.000 ,006 , 222 0.000 ,003 ,478 0.000 ,001 , 900 
ays 0.000 ,000 , 414 0.000 ,000 , 217 0.000 ,000, 174 0.000 ,000, 110 0.000 ,000 ,053 0.000 ,000 ,026 
bo l 1 ] ] ] ] 

be 0.326 ,666 , 667 0.289 , 166,667 0.277 ,777 ,778 0). 256 ,666 , 667 0.229, 166,667 0.206 ,6€6 , 667 
by 0.041 ,400 ,000 0.082, 501,563 0.029 ,938 , 272 0.025,358 ,333 0.019, 726, 5638 0.015,400 ,000 
be 0.002 ,712,190 0.001 ,895, 569 0.001 ,676,845 0.001 ,306, 163 0.000 ,888 , 788 0.000 ,600, 191 
bs 0.000, 107 ,937 0.000 , 067 , 226 0.000 ,057 , 138 0.000 ,040 ,977 0.000 ,024,451 0.000 ,014 ,326 
bio 0.000 ,002,881 0.000 ,001 ,600 0.000 ,001 ,307 0.000 , 000 , 863 0.000 ,O000 , 452 0.000 ,000 , 230 
bye 0.000, 000 ,055 0.000, 000 , 027 0.000 ,000 ,021 0.000 ,000 ,012 0.000 ,000 ,005 0.000 ,000 ,003 
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Spanwise and Chordwise Loadings for 
Rectangular Wings of Aspect Ratio near 


Unity 


William R. Laidlaw 

Department of Aeronautical Engineering, Aeroelastic and Structures 
Research Laboratory, Massachusetts Institute of Technology, 
Cambridge, Mass 

July 6, 1953 


, se LIFTING SURFACE INTEGRAL EQUATION has been formu- 
lated by Reissner! and for the rectangular wing may be stated 


l * dr dn 
4nb J_, dn (y — n) 
oy Vis — t+ —9F 


1 ' : 
4r Seed 1 On (x - 


In the present notation, the wing has a span of 2sb and a chord 
of 2b. The local vorticity y is related to the total pressure dif 


ference at a point by 


as follows: 
w(x, y) = 


didn (1) 


is — @) 


y = Ap/pU (2) 
The circulation I is the chordwise integral of the local vorticity 


atany spanwise station; hence, 


L(y) m 
— = Db v(é, y) dé (3) 
pl = 


In considering a wing of aspect ratio (s) near unity we may 
simplify the radical of Eq. (1) by the following considerations 


ry) = 


R= V(x — =)? + (y — 2)? = V ut + 0? 


where u = |x — ¢| andy = 





Therefore, 
R=(u +0) 


, uv I uv 2 | 
| (u+v)? 8L(u +0)? f 


Neglecting second-order and higher order terms 
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FORUM 


Now 
uv lj 
u+v 4 | 


Over the range of integration u and v are of the same order of 


j (5 
v)f ete 


magnitude, and, consequently, the second term of Eq. (5) may 
be dropped so that we have, finally, 
3 3 
R= i (s+) = (jx — & + ly — nl) = Ra 


4 (6) 


It may be noted by a comparison of Raq with R that for a fixed 
point (x, y) the maximum error (i.e., 25 per cent) occurs when 
The integral involving R is, however, singular 
y, and, consequently, when taking the prin- 


é—> x, Ory — ¥. 
at £ = x and 7 = 
cipal value, we find that near these singularities errors in R have 
a negligible effect upon the final result. 

Substituting Eq. (6) into Eq. (1) and simplifying, we have 


I dr dy 3 
- : - x 
l6rb > an(y — n) 8x 
> Wy dn 3 ; y(t, y) 
, dé 
ue oS {7 — 9) Sr - 4 (% &) 


z . . 
¢ = oX,n) = Pa y(t, ) dé 


w(x, y) = 


where 


We may next proceed to a solution of this equation by specifying 
a pressure distribution of the form, 
Ap 


= XY(%): V(¥) (8) 


pU? 
where the angular variables are defined as follows: 


x = —cos% y=scos¥ 


This assumed form of solution is restrictive in the sense that it 
makes no provision for the spanwise variation of section center of 
pressure, but from the standpoint of obtaining total lift and pitch- 
ing moment, spanwise and chordwise lift distributions, and 
approximate chordwise pressure distributions it should provide 




































































R=(u +i - [uv/(u + v)] (4) reasonable results. 
TABLE 1 
Fourier Coefficients for Flat Plate Wing in Symmetrical Flight 

Co Ci C2 C3 C4 Cs dCi C.P 

of of of of Af a dol c 
0.50 10000 | 0.0562 0.0117 0.7244 |-0.7848 |-0.8899 |-0.7643 |-0.6083 |-0.4969 0.8192 0.1220 
0.75 1.0000 | 0.0718 0.0156 0.8873 |-08225 |-0.8527 |-0.6732 |-0.5223 |-0.4148 LI747 0.1460 
1.00 1.0000 | 0.0864] 0.0195 1.0225 |-0.8312 |-0.8100 |—0.6024 |-—0.4768 |-0.3612 1.4974 | 0.16222 
23 1.0000 | 0.0998 0.0234] 1.1388 |-0826! |-0.7677 |—05469 !|—0.4164 |—0.3228 1.7910 01738 
1.50 1.0000; O1112 0.0269 1.2419 |-0.8141 |-0.7298 |—05022 |-03812 |-0.2938 2.0600] 0.1827 
1.75 10000] 0.1248] 00312 1.3334 |-—0.7987 |-0.6958 |—0.4706/-0.353! |-0.27/12 2.3046} 0.1897 
2.00 1.0000 | 0.1365 | 0.0351 | 1.4162 |-0.7812 |-0.6649 |—0.4348 |—0.3297 |—0. 2522 | 2.5306| 0.1953 
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Let 3.00 7 
| / 
_ = _ sin m¥ . ' ' 7 
V(¥) = }> Km (Ki = 1) (9) ones (ner. 21 | / i 
net m hey Z 
N ‘ 2.50 
ae cd . sin n® ‘ 
X(#) = c cot ~ oo z. fa (10) 
« n=1 n | 
Substituting Eqs. (8), (9), and (10) in Eq. (7) and performing the 2.00 ~ — 
necessary integrations, we have ‘i tena 7 
aaah 
2(V) ra C1 360 ; dA 
= 36 + oo + + 2(W) sin & + 
8s 2 2 8s 1.50 i el 
V in ( ie in (n +1) 4 YY LAWRENCE (REF. 3) 
3 Cn | sin (n — 1) sin (” D 
16 ini > | 1 7 +1 ] t Vf 
s : — n 
) n=1 n ‘i | VY | 
3 : cos nm ® w(P, V) 1.00 + — 
Y(v) a - : x ta - I— ( S : (11) VARIATION OF LIFT 
8 aa n U CURVE SLOPE WITH 
/ ASPECT RATIO 
where ‘sien / | | 
od _ sinmy ; Fi | 
av) = } Km ; | 
ina’ sin V / 
N ; i | 
: . sin (” — 1) ° 
dL, mm, etee~! 0 0.50 1.00 1.50 2.00 2.50 
= ASPECT RATIO (S) 
When the downwash function w(®, V) has been prescribed, the Fic. 1. Comparison of lift curve slope with the results of other 


unknown constants Ke, K3,..., Km, and @, 4, ..., €n may be 
determined by causing the above identity to be satisfied at a num- 
ber of points on the plan form equal to the number of unknowns 
(collocation scheme). 

Another and somewhat more desirable process is to use Eq. 
(11) in combination with a weighted integral procedure to gen- 
erate the required number of equations. 

For instance, if the integral operator 


LL 


(k = 0,1,2,... 


Ls k® sin V dbdv 


, N) 


is used on both sides of Eq. (il), a number of equations equal to 
the number of unknowns in the chordwise series may be de- 
These equations will involve the unknown c,,’s and the 
M 
quantity >» Km/m. We may designate this set of equations as 
m=1 
“‘Group I.” 
the operator 


veloped. 


Similarly, ‘‘“Group II’’ may be generated by using 


“rry 
4 > sin pW db dv 
LSI 
(p =2 


cod eee, 


In this set there are m equations equal to the number of unknown 
coefficients in the spanwise series, and once again the equations 
contain the unknowns c, and Kym. 

The coefficients are obtained by a method of successive ap- 
As a first approximation it is assumed that K,, 
= 0 for m 2 2. 
spanwise lift distribution. ) 
first approximation to the chordwise coefficients. 


proximations. 
(This is equivalent to assuming an elliptical 
Group I is then solved to obtain a 
These results 
are substituted into Group II to obtain the second approximation 
to the spanwise coefficients. The procedure is repeated until 
convergence is obtained 

The following aerodynamic properties of the wing may be de- 
rived readily by making use of Eqs. (8), (9), and (10). 

Spanwise lift distribution: 


Ll = pU*br[cyo + (a:/2)] Y(¥) (12) 
Total lift: 
L = pU?sb*( r?/2) [cy + (c:/2)] (13) 


theories. 


Section pitching moment about leading edge: 


—pU*b? = (« +aq- *) Y(¥) 












































m= (14) 
Total pitching moment about leading edge: 
eo nr? C2 
M = —pl 2sbh8 4 otan- 4 (15) 
Chordwise lift distribution: 
lL, = pU*sh(r/2)X(®) (16) 
0.30 
$ #00 
0.25 
—_. 
LAWRENCE (REF. 3) ——>- — 
aie a 
0.20 — 
Pa ae 
ea C__present 
,. THEORY 
£P 0.15 
Va 
“a , 
/ VARIATION OF CENTER 
0.10 L OF PRESSURE WITH 
. 7 ASPECT RATIO 
/ 
/ 
0.05 L 
/ 
/ 
/ $*0 
JONES (REF. 2) 
re) ——= © _— 
° 0.50 1.00 1.50 2.00 2.50 
ASPECT RATIO (S) 
Fic. 2. Comparison of the center of pressure location with the 


results of other theories. 
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Rolling moment: 


L = —pl?s*b3 Be («. a _) RK, (17) 
8 2 . 


Calculations have been performed for a family of flat plate 


rectangular wings at an angle of attack a, with aspect ratios 
ranging from 1/2 to 2. In this case, it is necessary to consider 
only the odd harmonics in the spanwise series; hence, 


Ko, K,, Kg, oree Kom i 0 


For a spanwise series of three terms and a chordwise series of 
six terms the results are presented in Table 1 Fig. 1 illustrates 
the variation of the lift curve slope dC,/da as a function of aspect 
ratio. Fig. 2 presents the effect of aspect ratio on the center of 
pressure location. In both figures comparisons with the results 


of Jones? and Lawrence’ are given. 
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A Note on Poisson’s Ratio* 


Morris Feigen 

Department of Engineering, University of California at Los Angeles; 
Now Hughes Research and Development Laboratories, Culver 
City, Calif. 

July 14, 1953 


OISSON’S RATIO DEFINES the ratio of lateral contraction to 

axial extension under an axial load and is constant as long as 
the action is wholly elastic. In the inelastic range Poisson’s 
ratio is no longer a constant as long as the axial load is increasing. 

In inelastic buckling of plates and columns it is assumed that 
the buckled state is a deflected state of equilibrium which is 
arbitrarily close to the uniform state of compression which pre- 
ceded it. In passing from the uniform state of compression 
to the buckled state, incremental axial extensions and lateral con- 
tractions (negative) occur. In this case the strains of interest 
must be the incremental changes that occur. In reference 1 it is 
these incremental changes that are used in the analysis of buckling 
of plates. 

This note proposes a new terminology in order that ratio of 
incremental contraction to extension be differentiated from the 
overall ratio of total contraction to extension. It also shows how 
the two above ratios may be expressed in terms of the variables 
of the uniaxial stress-strain diagram. 

If es is the total axial extension under an axial stress o, and if 
¢, is the corresponding total lateral contraction, then the usual 
definition of Poisson’s ratio is 


iy (1) 
This is defined as Poisson’s secant ratio, vs. 

If at any point on the stress-strain diagram an incremental 
Stress Ao, is added, there results an incremental axial extension 
Ae, and an incremental lateral contraction Ae,. The ratio of 
these incremental strains will be defined as Poisson’s tangent ratio, 


* This paper is an excerpt from a thesis submitted in partial fulfillment of 
the requirements for the degree of Doctor of Philosophy in Engineering, 


University of California at Los Angeles 
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Fic. 1. Poisson’s ratio for a Ramberg-Osgood stress-strain 
relation, Ee/oo.7 = (¢/00.7) + (3/7)(¢/ 00.7)". 
vr, where 
vr = Ae,/ Ae (2) 


Shanley? shows that Poisson’s secant ratio can be expressed 


as 
vs = (v'ex’ + vex”) /e (3) 


where e,’ is the elastic component of the total strain, ¢€,” is the 
plastic component of the total strain, v’ is Poisson’s ratio in the 
elastic range, and v” is Poisson’s ratio defining only the plastic 
extensions. The usual assumptions made are that e:’ + e.” = 
eand vy” = 1/2. 

In order to express Eqs. (1) and (2) in other terms, the following 


values are used: 


Or Or 1 l 
€r = ; €z' = ; é€z" = Or — 


Es Es E 


where £ is Young’s modulus and Es is the secant modulus. Sub- 
stituting the above into Eq. (3), Potsson’s secant ratio becomes 


Es Es 
psa Et ot _-—= (4) 


Similarly, if 
Aer = Ao,z/Er; Aez’ = Ao;/E 
Aex” Aoz[(1/Er) — (1/E)] 
Ae,’ = v’ Ae’; Ae,” = v” Ae” 


ll 


where Er is the tangent modulus, is substituted into Eq. (2), 
Poisson's tangent ratio becomes 


vr = v (Er/E) + v" [1 — (Er/E)) (5) 


vs and vr as expressed by Eqs. (4) and (5) are plotted on Fig. 1 
for a non-dimensional Ramberg-Osgood® stress-strain relation 
for v’ = 0.33 and v” = 0.50. The Ramberg-Osgood stress-strain 


relation is 
Ee/oo.7 = (a/00.7) + (3/7) (@/a.7)'" (6) 


where go.7 is the stress at which the secant modulus is equal to 
0.7 the elastic modulus. 
It should be observed in Fig. 1 that the tangent ratio ap- 


proaches !/2 more rapidly than does the secant ratio. 
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Two-Dimensional Correction of the Outlet 
Angle in Cascade Flow 


N. Scholz 

Institute of Fluid Mechanics, Technical University, Braunschweig, 
Germany 

July 14, 1953 


T IS WELL KNOWN that there are considerable difficulties in 
obtaining two-dimensional flow in cascade experiments, be 
cause of the boundary layers on the side walls of the cascade 
tunnel. These boundary layers generally have two different 
influences on the flow in the middle section of the cascade: (1) 
From the decay of lift at the side walls there arises a free vortex 
sheet that induces a secondary flow in the whole field. (2) Be- 
cause of the displacement of the boundary layers at the side walls, 
there is a contraction of the flow in the middle section of the cas 
cade, so that the two-dimensional equation of continuity is no 
longer fulfilled there. Both influences may be dealt with sepa- 
rately as long as they are small. 

In this note we want to deal only with the influence of the con 
traction in the middle section on the outlet angle of the cascade 
flow. The disturbances due to the vortex sheet have already 
been investigated by several authors.?~‘ 

In the normal cascade tunnel, the inlet flow comes through a 
closed channel, Fig. 1, so that the direction of this flow is given 
by the axis of the channel. Behind the cascade there are usually 
no tunnel walls so that there is a free jet. 
the side walls there is a contracting flow due to the displacement 
effect of the boundary layers on these walls (see Fig. lb). The 
rate of mass flow per unit of width and for one spacing of the cas 


In a plane normal to 


cade 


V= (1) 


f° w(y) sin B(y) dy 
eo V 
increases downstream for the middle section, Fig. lc, whereas 
the flow is deflected when going through the cascade. For a 
simple calculation we replace each of the blades by its load con- 
centrated at its aerodynamic center. The flow in plane 1 in 
front of the cascade and in plane 2 behind the cascade is assumed 
to be known from experiment. For reasons of simplification we 
assume that from plane 1 to plane c there is no deflection but only 
a contraction from Q, to Q,. Furthermore, we assume that in 
plane c the flow obtains an impulse in the y-direction, which is 
necessary for the deflection of the flow. Further downstream 
between planes c and 2 there is some more contraction from 
Q- to Q». But this second part of the contraction also gives an 
alteration of the outlet angle due to the momentum theorem 

We apply the momentum theorem for the y-direction for a con- 
trol surface in the middle section of the flow, Figs. la and 1b, 
which consists of two identical stream line planes having the dis- 
tance of the blades and two planes parallel to the side walls with 
the small distance Az. For the application of the momentum 
theorem it is not necessary to know the velocity distribution in 
plane c. It is sufficient that in plane c the incoming momentum 
in the y-direction is known. By definition, the incoming y- 
momentum is determined by the condition that the flow would 
attain the outlet angle 82 cor behind the cascade if there were no 


contraction downstream from plane c. Therefore, 


I, = p°Q, *W2 cos*COS B2 cor (2) 
Furthermore, the outgoing y-momentum in plane 2 is 
I, = p+ Q2-We-cos Bo (3) 


Because of the contraction of the flow between planes c and 2, 
there is an additional incoming momentum through the planes 
The incoming rate of flow is from the 
a — O.. 
in the y-direction is not known, b:.t it may be assumed to be be- 


parallel to the side walls. 


equation of continuity AQ = The velocity component 
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Fic. 1. Cascade flow for solid wall cascade tunnels. (a) Flow 
in the middle plane. (b) Flow between the tunnel side walls 
(c) Mass flow in the middle plane. 





tween those in planes cand 2. As an approximation we take the 
arithmetic mean value of these velocity components, and there- 
fore the incoming additional flow of momentum in the y-direction 
is 


AT = p(Qz — Q-) (1/2) (we cor COS Be cor + W2-COS B2) } 


Since there are no external forces on the control surface in the 
y-direction, the incoming flow of momentum is equal to the out 
going flow; therefore, 7, + AJ = J». Considering that 

Qc = t+ We cos’ SiN Be cor 

Q2 = 


one obtains for the corrected outlet angle, after some calculation, 


t-wo-sin Bo 





ctg Bz cor= (Q2/Q-) ctg Be 


Because ctg 8 = t-w-cos 8/Q, this means that the velocity com 
ponent in the y-direction is not influenced by the contraction but 
only by the x-component according to the increased rate of mass 
flow. For the evaluation of Eq. (5), not only a measurement of 
the wake in plane 2 would be necessary for obtaining Q» but also 
a detailed measurement for the flow in plane c between the blades 
at the position of the aerodynamic center, which is laborious 
As has been shown by unpublished measurements of the author | 


it is possible to relate Q, to Q; in plane 1. One obtains 
Oc - OV 


Q2 - O 


Kic); 0< Ke) <1 (6 
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ot fixed positions of planes 1 and 2, the.constant A(c) depends 
sly on the position of the plane c between the entrance and exit 
une. This means, that the contraction Q(x) — Qi, as indi 
ted in Fig. 1, is always approximately the same for a certain 
cade tunnel if it is related to the total contraction Q2 — Qi. 
the curve of the acceleration on the mean stream line is assumed 
,be independent of the local cross section as given by the dis- 
sacement due to the boundary layers on the side walls. This 
<evident if the distance of the side walls is not too small. 

The constant A, therefore, is of universal character and must 
e determined from experiments once and for all. The influence 
{the position of the aerodynamic center of the blade on K is also 
gnall for small total contraction. Therefore it is not necessary 
«know this position exactly; a rough estimate is sufficient. 

In some concluding remarks something may be said about how 
ther proposals for correcting the outlet angle are related to our 
results. Mortarino® assumes that the outlet angle is independent 
{the contraction. This would mean that there is no contrac 
tion between planes cand 2; therefore, Q; = Q2and K = 1. The 
other limit, A = 0, has been proposed by Katzoff, Bogdonoff, 
nd Boyet,® who take the x-component of the velocity of the exit 
cross section to be equal to that of the entrance cross section, 
(. = Q:. Another proposal' gives corrections for the entrance 
angle, as well as for the exit angle. This cannot be related to our 
result. Hawthorne’? from the free vortices obtains an x-com- 
ponent of the velocity at the position of the blade, which is the 
arithmetic mean value of the x-components in front of and be 
hind the cascade. This means 
Qe = (1/2) (Qi + Q2) 
andtherefore K = 1/2. 
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Mean Skin-Friction Coefficients for Biconvex 
Airfoils in Supersonic Flow 


J. Christopher Boison 
North American Aviation, Inc., Downey, Calif. 
July 13, 1953 


HE PURPOSE OF THE PRESENT NOTE is to show the effect of a 
chordwise pressure gradient (function of thickness ratio, 
t/c) on the mean skin-friction coefficient of biconvex airfoils in 
supersonic flight with a compressible laminar boundary layer. 
Weil’s procedure! was used in determining the variation of the 
local chordwise skin-friction coefficient on the airfoils, except that a 
linear chordwise velocity distribution was assumed. The present 
Tesults are shown in Fig. 1 as solid lines. Also indicated in the fig- 
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ure are results of an analysis that separately assumes the boundary 
layer to be incompressible, the local free-stream flow, however, 
still remaining supersonic. The increase in coefficient owing to 
compressibility in the boundary layer is quite evident. 


The ordinary nonlinear differential equation (Eq. (21), refer- 
ence 1) was solved numerically with the aid of a Magnetic Drum 
Digital Differential Analyzer (MADDIDA), where the ac- 
curacy of the solution was maintained at one part in a thou- 
sand. 


An approximation was used in the determination of mean 
skin-friction coefficients. The biconvex airfoil was assumed to 
be made up of two sections; the first 5 per cent of chord was a 
flat plate with incompressible boundary-layer flow, and the re- 
maining 95 per cent was a biconvex airfoil with a compressible 
boundary layer. Mean skin-friction coefficients were obtained 
separately for each section by integration of the local coefficients 
These mean coefficients were then weighted over a dimensionless 
chord, therefore giving an overall mean skin friction (Cry V R,,). 
The effect of R.., undisturbed free-stream Reynolds Number, was 
eliminated by ratioing the biconvex mean skin friction to the 
mean skin friction for a flat plate (Cpr V R,,) with an incom- 


pressible boundary layer. 


REFERENCE 
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On Theories of Plasticity and the Plastic 
Stability of Cruciform Sections 


P. P. Bijlaard and H. A. B. Wiseman 
Cornell University and The Pennsylvania State College, Respectively 
July 8, 1953 


A’ INTERESTING PAPER IN THIS JOURNAL,! which deals with 

the plastic torsional buckling of cruciform sections, is based 
on the contention that theoretical, as well as all direct experi- 
mental, evidence has shown that a simple flow theory would be 
valid rather than a deformation theory. Using simple flow theory 
and assuming initial distortions, the authors show that the re 
sulting buckling stress for a special case is in good agreement with 
that from deformation theory and, therefore, with experiments 
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In the first place we should like to point out that until now it 
has not been possible to predict, by pure mathematics or other- 
wise, how a highly complex polycrystalline material, with its 
flaws, dislocations, etc., will behave under the action of states of 
stress with varying ratios between the stress components. At 
most it may be concluded that for a quasi-isotropic material the 
actual behavior generally will be between that according to 
simple flow and deformation theories. The first writer calcu 
lated the plastic buckling stress from both these theories and 
showed that deformation theory leads to the smallest and thus 
to safe critical stresses.2 Later experiments appeared to be in 
good agreement with his results based on deformation theory, as 
was shown in reference 3. In certain cases deformation theory 
leads to impossible results, such as negative work done by posi- 
tive stress, as shown on page 734 of reference 2. However, in 
the special case of plastic buckling, where, in connection with the 
Shanley principle and small eccentricities, the shear stresses in 
all planes will increase during buckling, there is theoretically no 
reason why the final deformation should not be about the same 
as in a case where the ratio of the stress components remains con- 
stant—that is, in agreement with the deformation theory, as was 
stated on pages 532 and 533 of reference 3. Batdorf and Budi- 
ansky’s theory of slip,* ® which is theoretically sound, is based on 
the same principle and leads also to the conclusion that for plastic 
plate buckling the deformation theory leads to correct results. 


Although this theory was not confirmed by the special direct 
tests on compressed and then twisted tubes set up for that pur. 
pose,® these tests reveal that for finite superimposed shear stresses 
as occur in the case of buckling, due to unavoidable eccentricities 
the simple flow theory is also not correct. Moreover, recent 
tests conducted at The Pennsylvania State College’ with aluni- 
num alloy 14S-T4, show that over the entire range its behavior 
is in much better agreement with deformation than with simple 
flow theory. This is illustrated by Fig. 1, which is identical with 
part of Fig. 13, Set A, test 4, of reference 7, referring to shear 
stresses superimposed on pure tension. From more recent 
tests on 145-T6, the experimental curve (partly shown in Fig. 2 
usually falls to the right of both flow and deformation theories, 
so that in all cases the flow theory is much worse than the de- 


formation theory. 


It should be kept in mind that, in connection with unavoidable 
eccentricities, for the behavior in cases of buckling of cruciform 





sections it is not so important, as often implied, how the material 
behaves under superimposed infinitely small shear strains or 
stresses, as how it behaves under superimposed finite shear strains 
of the order of 0.001. For such shear strains, even from Fig. 12 of 
reference 6, the behavior becomes closer to deformation than t 
flow theory. This is in agreement with Dr. Budiansky’s feeling 
(from letter to first writer) that it would be unwise to pay too 
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Fics. 3 AND 4. 
constant. g;,* is the elastic, and a is the plastic, critical stress. 


o versus (0.°/E)/(o*/E) diagrams for a metal where, 


above the yield stress o*, the tangent modulus £; is a 
20; Fig. 4 referston = E/E; = ll 


Fig. 3 referston = E/E; 
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ych attention to the initial shear modulus as found from his 
sts. 

It would, of course, be nice if, as suggested in reference 1, flow 
eory would lead to the same results as deformation theory, but 
is is true only for the special case considered in reference 1 
20 and where 


om Fig. 9 of that reference, where n = E/E; 
10~%, the buckling load from flow theory 

This means that, 
ao versus o, /E diagram, with 2 = 20*/E, 
the critical stresses from both theories are indicated by the same 
However, from Eq. (20) of 
10~3, the buckling stress for 


o* = 2, with Yo 
jincides with that from deformation theory 
20 and for o,°/E = 


int P, as shown here in Fig. 3. 
ference 1, for 20 and Yo 
ther values of o-°/E is given by the dash line in Fig. 3, as calcu- 
ited conservatively for y = 0.05. Hence, for all other points 
than that considered in reference 1 (point P) there is no agree- 
The disagreement becomes more pronounced if n = E/E; 
Actually, the tests for 


nent 
ssmaller, as shown in Fig. 4 for n = 10. 
ruciform sections, as given in reference 8, refer to a stress-strain 
jiagram, which is practically identical to one for m = 10. More- 
wer, for equal initial angle rotations, yo is less for larger ratios 
‘’¢* than for smaller ratios, so that the disagreement with flow 
theory is still worse than is shown in Fig. 4. 
The obvious reason why (even by assuming initial distortions 
= 10~3, or 260) = 38’, which is ten times those that the authors 
f reference 1 assume as existent) the results from flow theory do 


not agree with the test results is, as explained above, that the 
material does not behave according to the simple flow theory but 
behaves in sufficiently accurate agreement with the deformation 
theory. Nevertheless, reference 1 is a welcome contribution to 
indicate that small deviations in the direction of the simple flow 


theory are canceled by unavoidable eccentricities. 
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Note on the Prandt! Number for Dissociated 
Air 


C. Frederick Hansen 

Aeronautical Research Scientist, Ames Aeronautical Laboratory, 
NACA, Moffett Field, Calif 

July 27, 1953 

I A RECENT ARTICLE, Moore! calculates a Prandtl Number 
for high-temperature air which varies from the usual 0.7 to 

values as high as 6, due to dissociation. Moore’s calculation 

depends on the assumption that the ratio of the coefficients of 

thermal conductivity and viscosity, A/u, is not affected by dis 

sociation. Crown? reports a different result—namely, that 

Prandtl Number for dissociated air is given by the same formula 


that Eucken deduced for normal gases,* 


FORUM 


Pr= 4y (97 — 5) 


In this expression the Prandtl Number, Pr, depends only on the 
ratio of specific heats, y, and cannot exceed unity. Crown's 
argument is based on the assumption that the difference between 
the heat capacities at constant pressure and at constant volume, 
Cp - 

According to kinetic theory of gases, the ratio K/y varies ap- 
proximately in proportion to C,,* § which is greatly increased by 
Thus, it is apparent from the definition of 


C,, is the gas constant R. 


dissociation. © 7 
Prandtl Number, 
Pr = C,p/K (2) 
that the values given by Moore are much too large. Crown’s re 
sult is also in doubt, since the difference, C, C,, for partly dis 
sociated air becomes many times the gas constant R.* 7 More- 
over, both Moore and Crown treat dissociated air, in effect, as a 
homogeneous gas rather than a mixture of atomic and molecular 
gases. Since mixing dissimilar gases can strongly influence trans 
port properties, such as viscosity and thermal conductivity,® it 
seems appropriate to investigate the atom-molecule mixing 
effect and to re-estimate the Prandtl] Number for dissociated air. 
A rather rigorous analysis of gas mixtures has been developed 
by Chapman,‘ Enskog, and others. However, this analysis re- 
quires knowledge of the properties of the mixture components, 
and for the atomic components of dissociated air such knowledge 
is uncertain. Therefore, the following discussion will be limited 
to a simple analysis, based on elementary kinetic theory, leading 
to a first approximation of the mixing effect. In view of the order 
of the approximation, air will be treated as a mixture of two pure 
components: molecular gas with mol fraction m, and atomic gas 
with mol fraction m». The coefficient of viscosity for the mixture 
may be expressed? as 
br 
oad (p: 12; + potiele) 


p= (3) 
where pi, Ui, and J; are, respectively, the density, mean random 
velocity, and mean free path of the ith component. If we invoke 
Eucken’s assumption that internal energy is distributed among 
the gas particles independently of their velocity distribution,’ the 
expression for coefficient of thermal conductivity becomes* °® 


ir 5 = ” 
K = se a Cit + Ci pit Ly + ra Ci! + C;* potigLe (4) 
32 2 “ 


where C;t and C2‘ are the translational heat capacities of the 
molecular and atomic components of the gas, Ci" is the internal 
heat capacity of the molecular component due to rotational and 
vibrational energy, and C2" is the heat capacity associated with 
dissociational energy of the atomic component. With the as- 
sumption that the mean free paths are the same for momentum 


and energy transport,° it follows that 


5 5 
, Cc + C,° rr Ct + C;' 
K 2 2 


——— (5) 
l+r 
where 
r = priiele/piiily 
In terms of the molecular weight ratio, 14,/ M2, 


= ( M, M) (ne /n;) 


P2/ Pi 
tio/ti; = (M,/M2)“ 
7 mSu V 2 + n2Si2 \ 1+ (M,/Ms2 


N2S22 V/2 + ny Si2 Vv ‘I ( M, 'M;) 


where 
Si; = x[(ox + o;)/2}? 
is the effective collision cross section for particles with diameters 


o; and a; 
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Figure |— Prondti number for oir os a function of temperature (pressure * 0! atmosphere) 
Combining Eqs. (5) and (2) defines the Prandtl Number as a 


function of the heat capacities, component mol fractions, molecu 
lar weight ratio, and the collision diameters, 
y¥(1 +7) C, 


¢ Cit + ci) +r (¢ Cot + c) 


The heat capacities and component mol fractions may be deduced 


Pr (6) 


from the equilibrium properties of dissociated air given in refer- 
ences 6 or 7, and /,/AM.is taken as 2, of course, for the dissociation 
of diatomic gases. The collision diameter ratio, o2/o;, is uncer- 
tain but probably lies in the interval from J/./\/,; to 1. This 
uncertainty does not influence the order of approximation, since 
the value of 7 is largely determined by the component mol fraction 
ratio, m2/n,, and is only weakly dependent on the collision diam- 
eter ratio or the molecular weight ratio as leng as these ratios are 
near unity. 

The Prandt! Numbers calculated from Eq. (6), using proper- 
ties of air taken from Krieger and White’ and a collision di 
ameter ratio of 0.5, shown in Fig. 1 of 
temperature of 10°! atmospheres. Similar 
results were obtained for pressures from 10° to 10 atmospheres, 
and only small differences occurred when collision diameter ratios 
of unity were used. 


are as a function 


for a_ pressure 


For comparison, the Prandtl Numbers 
given by Moore! and by Crown? are also shown in Fig. 1. Moore’s 
values are too large by a factor that is approximately the factor 
change in C, due to dissociation. On the other hand, the Eucken 
formula used by Crown |Eq. (1)] gives values of the same magni- 
tude as those calculated from Eq. (6), since neglecting the effect of 
dissociation on the difference C, — C, partly compensates for the 
neglect of the mixing effect. Admittedly, the atom-molecule 
mixing effect has been accounted for only to the first order of ap- 
proximation in the present analysis, and a more rigorous treat- 
ment ef the problem is in order. Nevertheless, it seems fair to 
conclude from this analysis that Prandtl Numbers for dissociated 
air remain in the region from about 0.5 to 1 
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A Double-Diaphragm Shock Tube to Produce 
Transient High Mach Number Flows 


Harry Bernstein 
Lewis Flight Propulsion Laboratory, NACA, Cleveland 
August 3, 1953 


— EQUATIONS GOVERNING THE OPERATION of the simple 
shock tube have been treated by several authors (e.g., refer 
ence 1). In the present note, the transient flow phenomena rf 
sulting from the combination of two shock tubes are considered 
This double-diaphragm shock tube consists of a long, constant 
area tube divided into three chambers (R, S and 7) by two thi 
diaphragms. Rupture of the first diaphragm is manual, while 
the secend is shattered when subjected to the pressure behind the 
strong shock caused by the disintegration of the first diaphragn 
The unsteady flow phenomena, for the required initial pressures 
Pr> Ps 

The gas originally in the middle chamber (region S of the x 
plane) is first compressed by the passage of a shock and then 


> pr, are depicted in the x-¢ diagram of Fig. 1. 


expanded isentropically by a rarefaction wave. Upon entering 
region s’ (the test region), the gas is at a high Mach Number an 
static temperature. The magnitudes of these quantities are, a 
in the case of the simple shock tube, wholly dependent upon the 


pressures, temperatures, and gases in the chambers prior to ru 
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ture 
Sample calculations have been performed for such a device 

employing a helium — air — helium gas combination, the results 

of which are presented in Figs. 2 through 4 

2 and 3 indicate that high Mach Number flows, free of 

gas condensation, may be produced. 


Figs 
For example, the use of 
diaphragm pressure ratios of 0.01 and 0.01 yields a transient flow 
of Mach Number 7.51 at a static temperature of 0.948 times the 
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temper: 
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The time duration of the} 


3 : { 
For the previously mentioned case, the duration ! 


of the flow is 1.01 X 10~‘ sec. per ft. of middle chamber length. | 
Hence, for a 50-ft. middle chamber length, a time duration 0 
5.05 millisec. will be obtained. 


gas temperature prior to rupture. 
transient flow is plotted in Fig. 4 as a function of the diaphragn 
pressure ratios. 


Time durations of this order ol! 
magnitude are sufficient for some types of model tests utilizing 
high-speed pressure or force-measuring instruments. 
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x-t diagram for double-diaphragm shock tube 
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The requirement of a middle section of such a large length, 
however, introduces a disadvantage in the form of large overall 
tube lengths. For the 0.01, 0.01 pressure ratio combination, the 
distance x, (see Fig. 1) was calculated to be 324 ft. for a 50-ft. 
middle chamber length. This is not the total tube length, for 
additional length must be provided to prevent wave reflections 
from the tube ends into the test region. Though no computations 
were performed, the indications are that the total tube length 
will be roughly twice x. 

The required tube lengths may be reduced, by a factor of ap- 
proximately 3, if argon, a gas of higher molecular weight, is used 
instead of helium in the R and 7 chambers. Such a substitution 
would, however, require slightly larger pressure ratios (Pr/ps and 
pPs/pr) across the diaphragms to produce the same Mach Num- 


ber. 
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A Note on the Contribution of Higher Mode 
Resonance to Bending in the Teetering Rotor 
Blade 


R. M. Carlson and D. R. Jacoby 
Structures and Aerodynamics Group Leaders 
Hiller Helicopters, Palo Alto, Calif. 


July 13, 1953 


Respectively, 


N A RECENT ARTICLE presented in the JOURNAL, Hirsch! dis 

cussed the resonance diagram for flapwise bending in a heli- 
copter blade. He pointed out that this diagram was invaluable 
as a means of predicting resonance in the preliminary stages of 
design, and so it is. It is the purpose of this note to discuss an 
interesting variation or expansion of the resonance diagram as it 
applies to the teetering rotor blade. 

The teetering rotor, as contrasted to the articulated rotor, is 
capable of performing the fixed-end modes of vibration, as well 
In general, this situation applies 
However, when 


as the pinned-end modes. 
when specifically considering a single blade 
considering the entire rotor, with a rigid attachment between 
blades, it is evident that there is a restriction that dictates how 
the rotor performs in the pinned and fixed-end modes. Of 
necessity, the rotor must behave in a symmetrical fashion when 
performing fixed-end modes and in an asymmetric fashion when 
performing pinned-end modes. Considering the fact that the 
teetering rotor behaves in this manner, it is convenient to utilize 
two resonance diagrams. One diagram plots the first, second, 
and third mode pinned-end rotating natural frequencies and a 
spectrum of the odd harmonic air loadings versus rotor angular 
velocity; another diagram plots the first, second, and third fixed- 
end rotating natural frequencies and a spectrum of the even har 
monic air loading versus rotor angular velocity 
For purposes of illustration, C.A.L. fiberglass blade No. 2 was 
considered as a teetering rotor, and the two resonance diagrams 
Fig. 1 is simply a portion 
For this illustration the 


were constructed (see Figs. 1 and 2) 
of the diagram presented by Hirsch.' 
effect of the eccentric flapping hinge on the rotational constant 
was not extracted. Fig. 2 utilizes the static natural frequencies 
employed in Fig. 1, modified to the fixed-end modes and the ro 
tational constants for the fixed-end modes.” 

It is of interest to note that, for this particular blade, by arbi- 
trarily considering it as a teetering rotor, its resonance character- 
The blade, as an articulated 
and fourth-order harmonic 


istics are considerably altered 
rotor, exhibited considerable second 
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content throughout the operating r.p.m. range. 
Fig. 2 indicates that as a teetering rotor this content would }, 


Inspection of 
appreciably reduced. It is also indicated that the eighth order 
third-mode resonance that was present on the Hirsch diagray 
is no longer present; however, it appears probable that th 
eighth-order resonance will manifest itself in the fourth fixed-ey, 
mode. It is evident from the above inspection of the resonances 
diagrams that it becomes increasingly more difficult to ayoj 
resonance in the operating r.p.m. range as the harmonic orde; 
increases. 

From the vibration point of view, the teetering rotor presents 
This 
note serves only to point out a situation that is not present j; 
the articulated blade but is of fundamental importance when a; 


no more or no less problems than the articulated rotor 


praising the teetering blade 
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Differentiation of Experimental Data by Means 
of Higher Order Finite-Difference Formulas 


Chi-Teh Wang and Daniel F. DeSanto 

Daniel Guggenheim School of Aeronautics, New York University 
New York 

June 13, 1953 


| A NOTE IN THE READERS’ Forum of the June, 1953, JOURNAL 

Kauffman and Shinbrot! presented a method for calculating 
the derivatives of obtained functions. The 
method consists of calculating two forward and two backward 
first-order finite differences at the given point with different in- 
tervals and then obtaining the ‘‘derivatives”’ by a graphical inter- 
The authors have demonstrated the accuracy of this 


experimentally 


polation. 
method by working out a test problem 

It appears to the present authors that good accuracy can also 
be obtained if the higher order finite-difference formulas are used 
in the ‘‘differentiation’’ of experimental data. This method has 
the additional advantage of being purely mechanical 
and 





Viz., no 
intermediate graphical interpolation is necessary there- 
fore, it is faster. 

To demonstrate this point, let us take the same example ana- 
lyzed by Kauffman and Shinbrot—namely, 


x = 6sin¢ + 0.5 sin 32 (0O<t< 7m) 


The values of x at ¢-intervals of 0.2 were computed as in the 
Kauffman and Shinbrot case. Using the values of x at two sta- 
tions ahead of and behind and at the station itself, the ‘‘deriva- 

T in terms of the higher order central differences 
are written as follows: 


dx ] om ee 9 
= Ax(T) — - A*x(T) 2 
dt J h 6 | 
d’x 1 = l ad 
= A*x(T) A‘x(T) 
dt? 1 h? 12 


tives” at? = 





where 
Ax(T) = (1/2) [x(T + kh) — x(7 h) 4 
Ax(T) = x(T +h) — 2x(T) + x(T — h) 5) 
A®x(T) = (1/2) [x(T 2h) — 2x(T + h) 4 | 
2x( T h) x(T —- 2h)) ( 
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t 
Fre. 1. 
Aix(T) = x(T + 2h) — 4x(T7 +h) + 
6x(T) 4x(T — h) + x(T — 2h) (7) 
and where hk is the ¢-interval between successive values of x; A, 


A?, A’, and A‘ denote the first, second, third, and fourth central 
differences, respectively. If the values of x at more stations are 
allowed to be used, Eqs. (2) and (3) can be extended still further. 


2) and (3), utilizing x-values at 


The derivatives given by Eqs. (2 
five neighboring points, are plotted with their exact values given 
The derivation of 


by (1) in Fig. 1. The agreement is excellent 
Eqs. (2) and (3) can be found, for example, in reference 2 

It is to be noted that it becomes impossible to use the preceding 
equations for points in the vicinity of the end points on the curve; 
this failing, however, also occurs in the case of the method pro- 
If it is necessary to evaluate derivatives at 








used 
d has 
.., 10 | 
here- 


ana- 


| the 


Sta- 





nces 





posed in reference 1. 
such points, then the derivatives must be expressed as series in- 
volving either backward or forward differences, depending on 


which end of the curve the evaluation is to be carried out. At 
the initial end of the curve, 
dx ] - A’x( T) A*x(T) A‘x(T) 
= AxT + -_— (8) 
dt }7 h 3 3 H 
d?x' l es 1] oe 
= = A*x(T) A‘x(T) +4 A‘x(T) (9) 
dt? J 7 h? 12 


where A, A’, A*, and A‘ denote first, second, third, and fourth 


forward differences, respectively, and are 


Ax(T) = x(T + h) — x(T) (10) 
A2x(T) = x(T + 2h) 2x(T + h) + x(T) (11) 

A(T) = x(T + 3h) 3x(T + 2h) + 
3x(T + h) xT) (12) 

A‘x(T) = x(T + 4h) tx(T + 3h) 4 
6x(T + 2h) in(T +h) +x(T) (13) 


In the case of the terminal end, the backward-difference formu- 


a 


las are 
(= 
dt 


] | ’ A2x(T) A3y( T) Atx(T) 
Ax(T) + { | > 
1 y 3 } 


(14) 
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FORUM 


d*x l 11 
( = A%x(T) + 43x(T) 4+ Aty( T) (15 
dt? Jz h? 12 


where 4, A’, A®, and A‘ denote first, second, third, and fourth 


backward differences, respectively, and are 


Axn(T) = x(T) — x(T — h) (16) 

A(T) = x(T) — 2x(T — h) + x(T — 2h) (17) 

Aiy(T) = x(T) — 3x(T — h) + 3x(T 2h) — x(T — 3h) (18) 
Atxy(T) = x(T) — 4x(T — h) + 6x(T — 2h) 

4x(T 3h) + xT th) (19) 


The derivation of the higher order forward- and backward- 
difference formulas can also be found in reference 2. The de- 
rivatives given by Eqs. (8), (9), (14), and (15) are plotted in 
0.0, 0.2, 2.8, and 


Fig. 1, together with their exact values at ?¢ 
3.0, again using five neighboring points as in the case of the cen- 
tral-difference formulas; the agreement appears satisfactory but 
not so good as that when the central difference formulas can be 
used. 

The authors found that in employing this method for the case 
in which x-values at the various values of ¢ were furnished correct 
to three decimal places, the root mean square errors in dx/dt and 
d*x /dt? were 0.058 and 0.593, respectively, for the condition where 
forward and backward differences were used to calculate deriva- 
tives at points near the ends of the curve. In instances where 
inner points only were investigated, the accuracy was consider- 
ably improved; the root mean square errors in dx/dt and d*x/dt? 
were found to be 0.055 and 0.037, respectively, for the condition 
where forward and backward differences were used to calculate 
derivatives at points near the ends of the curve. In instances 
where inner points only were investigated, the accuracy was con- 
siderably improved; the root mean square errors in dx/dt and 
d?x /dt? were found to be 0.055 and 0.037, respectively. 
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Note on Model Effective Span for Wind-Tunnel 
Wall Corrections 


Alan Pope 

Supervisor, Experimental Aerodynamics Division, Sandia 
Corporation, Albuquerque, N.M 

July 6, 1953 


N IMPORTANT PARAMETER in the determination of the wind- 
tunnel wall correction factor 6 is the ratio of the model span 
to tunnel width. Accordingly, agreement on the proper span to 
use can result in correction uniformity and better correlation in 
the results of tests on a given wing model in several tunnels 
Difficulty arises in selecting the model span. It is well known 
that the trailing vortex sheet shed by a lifting wing rolls up into 
two trailing vortices whose span is somewhat less than the geo- 
metric span of the wing. The question is, hence, whether to use 
the geometric span, the vortex span, or some value in between 
when setting up the image system, which in turn determines 6 
Since the wing is to be replaced by a simple horseshoe vortex 
that simulates uniform loading, there appears little reason to 
assume the vortex span to be equal to the geometric wing span 


Silverstein in reference 1 suggested reducing the geometric span 
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TABLE | 
Comparison of Wall Corrections, Eq. (1) and Other Sources 
Aspect Taper Aa, 
Source Ratio Ratio Plan Form Tunnel Aa Eq. (1) 
RM AdS2L04 $ 0 Delta 6x 6 0.55C, 0.554C, 
RM AS2L02 2 l Rect. 6x6 0.55C, 0.55C, 
, 0.3338 Swept 6 xX 6 0.55C) 0.55C, 
2 0 Delta 6X 6 0.98C, 0.92C, 
RM A5d2L 15a 3.1 0.388 Rect 6X 6 O.57C;1 0.555C1 
3.0 0 Delta 6x 6 0.55C, 0.55Cr 
3.0 0.40 Swept 6x 6 0.55C, 0.555C, 
L-603 §.87 0.66 Rect. «XX 1.02C, 1.04C; 
b, bg b, be |Eq. (1)] 
TR 770 1.0 xD 0.98 0.94 
0.5 ‘xD 0.88 0.89 
0.25 xD 0.83 0.86 


by a factor of 0.85 for rectangular wings and 0.8 for tapered wings 
This essentially suggests using the vortex span instead of the 
9 


geometric span. However, Swanson and Toll in reference 2 sug 


gest values somewhat larger than those of reference 1. In order 
to provide a simple but rational approach to the problem, the 
writer suggests assuming that, inasmuch as the complete roll-up 
occurs somewhat downstream of the model station, a reasonable 
compromise would be to select a value midway between the geo 
metric and vortex span-—that is, the span that should be used 
for computing the model span-tunnel width ratio should be calcu- 


lated from 


be = (bg + by) /2 (1) 
where 
be = effective span 
bo = geometric span 
by = vortex span 


Eq. (1) at least has the advantage of making an allowance for the 
widely different vortex spans experienced for a wide range of as 
pect and taper ratios. Still further, numerical checks of the wall 
corrections actually used in many tests reported by the NACA 
and others indicates that Eq. (1) is closely in line with constants 
(See Table 1.) 

Fig. 1 has been provided both to yield vortex spans for a wide 


used. 


range of plan forms and to illustrate the wide variation possible 
for the useful range of aspect and taper ratios. Use of Fig. 1 


(which was calculated by the method of reference 3 and is repro 
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duced from reference 4) and Eq. (1) for all wing tests will result 
in corrections that are both uniform and reasonable. Since this 
effective span is used for all the horseshoe vortices in the image 
system, it is equally applicable to corrections for wall effect aft 
of the wing. It also appears satisfactory for swept models 
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Hypersonic Viscous Flow over an_ Inclined 


Wedge? 


Lester Lees 
Associate Professor of Aeronautics and Applied Mechanics 
California Institute of Technology, Pasadena, Calif. 


July 29, 1953 


Dr. Shih-I 
vious work on hypersonic viscous flow over a flat plate? * to 
Unfortunately, his 


N A RECENT NOTE,! Pai attempts to extend pre- 
the more general case of an inclined wedge. 
results for the compression side of the wedge at large angle of 
attack are incorrect. The present note tries to clarify the situ 
ation and to define the regions of ‘‘strong’’ and ‘‘weak’’ inter 
action for this problem. 

Without going into detail, it turns out that the tangent-wedge 
approximation gives a good representation of the relation be 
tween pressure p* and flow deflection angle 6 along the surface 
of a body with an attached leading-edge shock in hypersonic 
rotational, planar flow.‘ Therefore, in the viscous flow problem 
p* and @ at the outer edge of the viscous layer are connected by 
the usual oblique shock relations. where @ = @2 + (dé/dx* 
Here, d5/dx* is the ‘‘induced”’ flow deflection and 6, is the slope 
of the body surface, equal to 6m + a@ for a wedge, where 4 is the 
wedge half angle and a is angle of attack. In considering the 
compression side of a wedge at large angle of attack Pai! relates 
P 
sum of 6, the shock angle without induced effects, and 6, the 


* to the corresponding shock angle 6s, and represents @s as the 


viscous effect. In this problem the quadratic term involving 


#2 cos? @ is small in comparison with the other two terms in the 


+ This work was supported partly by the USAF Air Research and De 
Contract AF-33(038)-250 (Princeton Un 
and partly by the Air Force and Army Ordnance under Contract 
GALCIT) 


velopment Command, under 
versity 
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xpression for p*/p.* = p, or 0,2 cos? 0 <K sin? O% |Eq. (19), 
ex] ] ] ] 1 


reference | But this condition rules out the possibility that the 
linear term involving 6, sin 26) could ever be large compared with 
the inviscid term. In other words, the pressure can never be 
proportional to 6, or to dé/dx*, and therefore the “‘similarity”’ 
condition for the velocity profiles and the associated relation 
§~x* ‘found by Pai do not, in fact, exist 

Actually, the hypersonic viscous flow on the compression side 
of a wedge at large angle of attack is a more general case of the 
weak interaction on a flat plate at a = 0° treated by Probstein 
and the present writer’ and by Maslen.* Here it is not the pres 
sure but the pressure perturbation, which is proportional to 
di/dx*. Although the flight Mach Number is high, the inviscid 
flow Mach Number is reduced considerably across the leading- 
edge shock, and the interaction effects are much smaller than for 
at the same local Reynolds Number and Mach 


Except possibly very near the lead 


the flat plate 
Number ahead of the body 
’ dé dx ° 


ing edge (see below 6n, and 


p* l dp dé 

: 2} 4 : 4 

pPr* p d0/ 6 =p adx* 
where the quantity [(1,/p)(dp/d6)|o 
oblique shock relations, and dé/dx* is given by the usual bound 


én is obtained from the 


ary-layer solution for zero pressure gradient, or 


dé 
dx* 


= dg V Cpl p?/V Rei 


The quantity dz is a function of the gas properties, surface tem 
perature, and Mach Number; Cz is the constant in the linear 
viscosity-temperature temperature relation; and ./,, Re, are 
Mach Number and Reynolds Number, respectively, based on 
inviscid flow quantites at the outer edge of the boundary layer 


“far” from the leading edge. In particular, if 


M.6n = M.(0w + a)>2 
(approx.) then 
p*/pe* = 1 + 2dp(Xe/Ae) +..., 
where 
xe = V Cp (M5°/V Resp 
Kp = Mp6p 


In this case solutions of the Prandtl boundary-layer equations are 
obtained in the form of asymptotic series in ascending powers of 
xe for any arbitrary thermal conditions at the wedge surface 
In order to compare this result with the flat-plate case, we recog 
nize that the viscosity ratio 


up/po* > JMO? 


if a linear viscosity-temperature relation is employed, and also 


M 6x — const., so that 


XB Ap* dx 
—_—— and — Fo — 2 
K; k? ps* K? 
where xy = VC M V Re, and K = M.6z. As indicated by 


the factor 2/A?2, the interaction effects are reduced considerably, 
and the solutions can be pushed to much lower local Reynolds 
Numbers than for the flat plate at zero angle of attack at the 
same flight Mach Number. T 

at not too high angles of attack, when the 
15) and the local Reyn 


For thin wedges 
Mach Number is sufficiently high (17 
olds Number is sufficiently low (Re, < 10° 
and \f..@ > 2 simultaneously, and we are dealing with the oppo- 


then dé/dx* > 6x 


weak interaction region is characterized by the condition 
small angles of attack 
In this case M9 = 1 


These question 


T In general, the 
XB OC and s« 
al Reynolds N 


and dé/dx* can be less than 


it exists also on thin wedges at 


when the lox imber is sufficiently high 


equal to, or greater than 4% 


and others are discussed in a comprehensive monograph on hypersonic 


viscous flow by R. F. Probstein and the present writer, now in preparation 
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FORUM 


site extreme case of ‘‘strong’’ interaction treated in references 


2,3,and 7. Except for the special case of the flat plate at a = 0° 
the geometric flow inclination 6, cannot be completely neglected 
in comparison with dé/dx*, as in the first approximation to the 
strong interaction region. For the more general case of the in 
clined wedge the generating parameter is suggested by considering 


the p*-@ relation for \/ 04> 2, 


15 \? On 2 37 + 1 
p= Jy + Me (4 YQ : ) } 
2 dx* di/dx* y+ 1 


87 l 1 
(y + 1)§ Ma%(dé a , Op ) 
db /dx* 
But 
Op M oO g A 
di/dx* = M(dé/dx*) ~~ -y/, 


where the flat-plate solution® 


6 : V Me 
x ( Rez ) 


is utilized to evaluate dé/dx* approximately. Therefore analyti- 
cal solutions of the Prandtl boundary-layer equations for the 
strong interaction region on an inclined wedge are obtained as 
asymptotic series of the following form:t 


pik po + psK? 
p = Pox 1 + } 
V x x 
u*(x, y) K 
= lola) + ula) + oo, Ge 
Ua* \ x 


where o ~ y*/(x*) ‘, and pi, po, ps are constants. The function 
uo(o) (and the corresponding 7) is the flat plate solution given 
in reference 3, while «:(0), 7\(o), etc., each satisfies a linear sec- 
ond-order ordinary differential equation with coefficients de 
pending only on the previous terms in the expansion. For 
Prandtl Number of unity the local heat transfer coefficient (for 


example) is expressed as 
ak ado + 
Vv x x 


»\0 a;k? 
= 


V: 


where Cy, is the classical value given by 0.332 VC 
and p'” = pox. 

So far as the expansion side of a wedge at large angle of attack 
obtained by 


Cu = Cu 


V Re, *, 


is concerned, preliminary experimental results, 
Nagamatsu and coworkers in the GALCIT 5-in. hypersonic tun- 
nel at .\f = 5.8, show that the pressure rise across the trailing 
edge shock is communicated upstream and causes separation of 
the laminar boundary layer almost from the leading edge. Thus 
it is doubtful that the usual boundary-layer theory can be ap- 
plied to this problem. 

All of the above results deal with the ‘‘direct problem" for a 
S. Tsien, 


“inverse 


body of simple geometry. At the suggestion of Dr. H 
Ting-Yi Li and the present writer recently solved the 
problem”’ of determining the airfoil shape required to produce a 
prescribed including hypersonic viscous 
effects. An the Karman-Pohl 


hausen method is employed and numerical computations are 


pressure distribution 


iterative scheme based on von 


now in progress. 
It must be emphasized that practically all of the current theo 
retical work on this problem is based on the idealization of a 


mathematically sharp leading edge. Since the leading-edge 


t Shen’s results for the insulated wedge obtained by the von Karman 


Pohlhausen method can also be expressed in the form of asymptotic series 


in ascending powers of 1/YU x [Eqs. (29) and (30) of reference 7, trans 


formed } 
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radius is always small but finite in real life, the above results are 
expected to be valid when the Reynolds Number per linear foot 
along the body surface is not too high. Otherwise, the hyper- 
sonic viscous effects are submerged in the inviscid flow effects as- 
sociated with a finite leading-edge radius of curvature. At high 
Mach Numbers, these effects may be propagated far downstream 
of the leading edge. Theoretical and experimental studies are 
now required to settle this problem. 
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On Strong Interaction for the Hypersonic 
Boundary Layer on an Inclined Wedge 


Shih-I Pai ; 

Associate Research Professor, Institute for Fluid Dynamics and 
Applied Mathematics, University of Maryland, College Park, 
Md. 


August 13, 1953 


iy THE PRECEDING NOTE,! Dr. Lester Lees criticizes my analysis 
for hypersonic viscous flow on the compression side of a wedge 
at large angle of attack.? Unfortunately, his criticism is based on 
his misunderstanding of my analysis, which was only briefly men- 
tioned in reference 2. In this note, I give a little more detail 
of my analysis in order to clarify such a misunderstanding and 
to show the relative region of application of my results and of 


those of Dr. Lees. 


For the compression side, the pressure is given by Eq. (19) of 
reference 2. Since the quadratic terms of Eq. (19) of reference 


2 are in general negligible and 6 is proportional to dé/dx, we 


~ 


have for the compression side the following expression of pres- 
sure: 
Pis« es db 
= yM_.,? sin 260 - (1) 
p p p dx 


For such a pressure distribution, I have shown in reference 2 that 
there is no similar solution. Dr. Lees’ analysis of reference 1, 
strictly speaking, holds only very far downstream. 
*(d6/dx) is not negligibly 


For a strong 
interaction region where the term J/ _ 
small, we cannot use the method of Lees. Some other suitable 
method should be used. One of the simplest methods, which 
would bring out the essential feature of the problem, is that 
given by Dr. Shen in his original analysis of the problem,’ in 


which for first approximation the fundamental equation is as 
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follows: 


n?p* = (2/m)é (2) 


where y = 6/L,§ = x/L,L = (M4, . *uw) (P4 M4), P*=p+ 
Riu and mis a constant. Eq. (2) is the same as that given by 
Shen in reference 3, but the characteristic length Z for the non- 
dimensional quantities » and é is different from that given by 
Shen. 
the compression side of a wedge at an angle of attack, the ex- 
pression of Eq. (1) may be used to calculate p*. 


Eq. (2) may be solved for any pressure distribution. For 


Substituting Eq. (1) into Eq. (2), we have 


n?[(dn/dé) + A] = Bé (3 


where A and B are constants for a given problem. 
Eq. (3) has a singularity at the origin. Near the origin we 


have 


ue (4 


errs 
A change of variable may now be made. Let 

= 1/; - 

g=& (5 


The solution of Eq. (3) may be written as follows: 


n= s? Do ang” (6 


The coefficients may be determined by substituting Eq. (6) into 
Eq. (3). Further improvements may be obtained by the pro- 
cedure outlined by Shen in reference 3. We will omit them 
here. 

It is easy to show that all the examples of similar solutions 
which I gave in reference 2 are the corresponding leading terms 
of the series similar to Eq. (6) for the corresponding expressions 
of p*. It shows the behavior of the singularity at the origin for 
the corresponding problems. That is why I mentioned briefly 
the similar solution for the compression side in reference 2 
Mathematically speaking, these solutions belong in the same cate 
gory. Unfortunately, Dr. Lees erroneously assumed that such 
a solution is applicable for the whole region of strong interaction 
All of his criticisms are based on an assumption that I did not 
make! In fact I do not think that any of the leading terms are 
good enough to be used for the whole region of strong interaction, 
including the case of a flat plate at zero angle of attack for which 
Dr. Lees made such an assumption in his analysis. According 
to my analysis [Eq. (11) of reference 2], the second term of the 
pressure is proportional to M _ . (dé/dx), which is not negligibly 
smallin that analysis. Further investigation of the effect of such 
a term should be made. 

In conclusion, my analysis is correct, and it discusses the case 
of strong interaction, while Dr. Lees’ results apply for the case of 
weak interaction. It is interesting to notice that Dr. Lees’ im- 
pression of my analysis is the same as his impression of Dr 
Shen’s original analysis‘ of this problem in the earlier stage of the 
investigation. Since the investigation of hypersonic viscous 
flow is still in a relatively early stage, I think that as time goes 
on all of the misunderstanding will be cleared up. 
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Effect of a Turbojet Engine on the Dynamic 
Stability of an Aircraft 


Hans C. Vetter 

Fl Segundo Division, Douglas Aircraft Company, Inc., El Segundo, 
Calif 

july 6, 1953 


SUMMARY 


It has been frequently noted that the calculated lateral damping of a 
turbojet-powered airplane is not always in complete agreement with flight 
This paper gives a summary of three major effects of a turbojet 
three effects are 


tests 


engine on the dynamic stability of an aircraft. These 


a) thrust component in the direction of wing lift; (b) ~Coriolis’’ reaction 


of air and exhaust gases in jet duct; and (c) coupling of lateral and longi- 


tudinal modes due to gyroscopic action of the turbine and compressor All 
these are essentially low-speed effects. 
(I) NOTATION 
x, 9,2 Cartesian coordinates 
¢, 0, disturbances and angle of roll, pitch, and yaw 
a,B disturbances in angle of attack and sideslip 
flight-path angle 
acl angle between jet duct centerline and rolling axis 
ar angle between gross thrust line and rolling axis (a7 acy if 
tailpipe is undeflected) 
q free-stream dynamic pressure 
Ss wing area 
b = wing span 
be mean aerodynamic chord of wing 
I moment of inertia tensor 
J polar moment of inertia of turbine and compressor 
W weight of airplane 
V true air speed 
Cr lift coefficient 
Cp drag coefficient 
C rolling moment coefficient 
Ce pitching moment coefficient 
Ge yawing moment coefficient 
Cr ram drag coefficient (ram drag Crq5) 
Cr gross thrust coefficient (gross thrust CrqS) 
v velocity of air in jet duct 
w engine r.p.m. 
a 
R xi + yj + 2k position vector 
—> 
Q di +07 + vk angular velocity vector 
—- 
G moment vector 
— 
H angular momentum vector 
i,j,k unit vectors 
p density of air 
A cross-sectional area of jet duct 
c distance measure (along jet duct axis) 
Ar increment due to thrust 
u() unit step function 
Subscripts 
p partial differentiation with respect to b¢/2V 
q partial differentiation with respect to ¢,6/2V 
partial differentiation with respect to by /2V 
1 static, trimmed condition 
The airplane axes used throughout are the conventional NACA stability 
axes 
(II) PowEer-On Lirt 
HIS EFFECT HAS, in the past, been generally passed off as in- 
significant. This is not so for aircraft with no flaps and rela- 


tively small lift-curve slopes. This effect amounts to a correction 


of the trimmed lift coefficient calculated from 


Cz = (W cos y¥1) ns (1) 


By considering the thrust in the force geometry of the airplane 


(see Fig. 1), we get 


; = Cy, + Cr, sin ar 
qi. ‘ 
, ( 


Co) 


W cos y1 | 


W sin y1 : 
Cr, 


= Cr, cos ar 
qs 
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194 


FORUM 





LIFT 





GROSS CENTER OF GRAVITY 


THRUST 








HORIZON 
— 
~ DRAG 
Gc 
~YRo 
SSs My 
Rus, 
“Zr 2 
‘Ne 
WEIGHT 


Z 
Fic. 1. Force diagram for trimmed flight. 
A solution of Eq. (2) for ar for a power-approach condition can 
as a 2° discrepancy with the value of a7 calcu- 
This means a change in the flight-path iner- 


show as much 
lated from Eq. (1). 
tias and in the stability derivatives. 


(III) ‘‘Corroits’”” DAMPING 


Those who have recognized the presence of this effect have 
probably ignored it and, in most cases, rightfully so. This effect 
can be significant only in airplanes with poor damping in yaw 
when flying at low indicated airspeeds 

“Coriolis’’ damping occurs because the air in the jet duct is 
traveling in a radial path with respect to the center of gravity of 
the airplane when the airplane is performing rotational oscilla- 


tions. By using the coordinate system of Fig. 2 and the formulas 
~t=d—> —<_ whe 
“Coriolis” moment = —2 RX(2X R)pA dt (3) 
é h 
eee ee 
“Coriolis’’ force = —2 (2 & R)pA dt (4) 
t h 
p(¢) A(¢) o(¢) = (rate of airflow) X [1 + ru(¢ +21] (5) 
we get 
ile as | 
ArG, = —2¢ Rie sin? ac, + 2Cr —, sin 2acr 
2 2 
ArC ae a 2 1C 2 | 
rli, = —Cr pe sin Zac, 4+ Cr he cos? acl 
. = a 2 a 
ArC,, = —Cr sin Zac, — 4CR sin? ac, 
b? b? 
A? 2 . 
ArCn, = —2Cr cos? ac, — 2Cr — sin 2ac; (6) 
b? b? 
Ary, = 4Cp(d? ‘ab) sin act 
ArCy, = 4Cp(d2 ab) cos acL 
ArCm, = —2Cr(A?/tw?) 
ArCrz, = 4C p(d2/ aly) COS acl 
ArCp, = 4C p(r? ‘at, ) sin acl J 
where 
A? = (h? d?) + r(h? — PP 
2 =al(h +d) +7(h—-))] 
rate of fuel flow (Ibs. per sec.) 
r = 


rate of airflow (Ibs. per sec. 
The most important one of these is A7C,,. ArCy, has been 
known to amount to as much as 12 per cent of the aerodynamic 


value of C,,, at high lift coefficients 
(IV) Gyro-CouPLING 


This is another effect, such as the first mentioned, which has in 
all probability been recognized but passed off as insignificant 
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However, with the trend toward larger and heavier jet engines, 
this effect is making itself felt more and more, especially in lighter 
type aircraft, 


In order to include-the gyroscopic terms in the equations of 


motion, we return to the basic equations of dynamics—namely, 
DH/Dt=H+QxH=G ( (7) 

H=I2+Joa { 

H, = Ir _ Tr + Jw cos acr, 
Hy = Iyé (8) 


H, = Tw _ I. — Jwsin act 
(Izy = Iyz = 0) 


By substituting Eqs. (8) into Eqs. (7) and neglecting all products 
of angular rates (perturbation method), we get 

Iz2@ — Izep — Jw sin acr6 = CiqSb 

Gy = Tyy6 + Jw cos acty + Jw sin acrd = CrgStw >} 
G, = Tw — Ta — Jw cos ac. = CrqgSb 


G; = 


(9) 


If we consider the quantity Jw cos acy, which appears in the yaw- 
ing equation of motion, as a yawing moment per unit pitch rate 

then the value of C,,, for one of the modern jet airplanes 
This is a defi- 


i.e., Ca, 
at low indicated air speeds is twice as large as Cp, 
nite indication of a sizeable coupling between the lateral and 
longitudinal modes of motion. This coupling, when included in 
the equations of motion, leads to the following general results: 
(a) an increase in the calculated lateral damping at low indicated 
air speeds; (b) a small decrease in the calculated longitudinal 
damping; (c) an insignificant change in the frequency of both 
modes; (d) small excitation of either mode when the other is 
disturbed—small enough, in fact, that a pilot seldom if ever 
notices it; (e) a small difference in the first few seconds of a calcu 
Fig. 3 shows the relative effects of 
coupling on the calculated damping for a modern jet airplane 


lated rolling time history 
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Note on Nonstationary Slender-Body Theory 
J.D. Cole 


Guggenheim Aeronautical Laboratory, 
Technology, Pasadena, Calif. 


July 14, 1953 


California Institute of 


AY” OF REVOLUTION performing arbitrary motion along 
its axis can be represented within linearized (acoustic 
theory by a distribution of sources that vary in time along thy 
axis of flight. It has been shown by Frankl! that for very slende 
bodies the instantaneous source strength is equal to the velocity 
times the rate of change of cross-sectional area along the body 
In this note an asymptotic expansion of the potential near th: 
axis is given, and the results for a special case are presented. 
The potential is represented by an integral of the retarded 
values of the source strength S(x, ¢), 
op Me eer 
ra ae, eo dé 


1 ri 
(1 
ae V(x — t)2? + 72 


in coordinates where the air is at rest at infinity (reference 2, p 
ee 
524). 


integral in Eq. (1) is carried out over only those values of £ which 


®x,7,t) = —- 


The source strength is zero outside the body so that th 


are common to the flight path and the surface of the upstream 
Mach cone 


t— [V(x — £)? + r2/c] 


T= 


as shown in Fig. 1. There are integrals for the recent past x; to 


x, and the distant past x; to x4, etc. 


The integral for the recent past can be analyzed for small r 
by splitting it into 


x—e(r) x-+e(r) "x2 
, , 
x1 x—e(r) x-+e€(r) 


since most of the contribution near the body comes from & = x 
For the middle integral, S can be approximated by its value at 
(¢ = x,r = 0); for the others, expansions of S and the square 
The integrals for the distant 


0, so that they 


root can be made in powers of r. 
past can easily be approximated by putting r = 
represent waves that are plane at the body. This procedure re 
sults in an asymptotic formula valid near the body for small r, 


oe r 1 , x— X% 
S(x, t) log = — S| xm, t — x 
Qn 2 Sor c 
Xe — x , 
log (x — x1) +S (x1 ) log (xe — »| 
P 
l * x2 ‘/ ; x — £| . 
S21 ét — + sign (x — £) X 
tor ss ri ti 


P(x, r,t) = 


x — € 
oe (« t— ° ) 
1 x4 C 
/ dé, etc., + O(e log €) (2 
x F 


tar | 


where x;(x, ¢) and x(x, t) are the values of £ in the recent past at 
which r = t# [(x — £)/c], respectively, intersects the boundaries 
of the flight path (Fig. 1); x; 
distant past; and e = e(r) such that /r+>e 

It follows from Eq. (2) and the boundary condition 4, = | 


and x, are similar values in the 
>rasr—0. 
y 


slope that 


t 
174, OA ™ 
S(x, t) = U(t) (x,t) = U(t) A x 4 U(r)dr (3 
Ox 0 
where 
U(t) = velocity of body in negative x-direction 
A(x) = cross-sectional area distribution of body at ¢ = 0 
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Just as Eq. (1) includes the source distribution for steady sub 
sonic or supersonic flow as special cases, given by a translation of 
axes, Eq. (2) includes the usual steady slender-body results (for 
example, Eqs. (10) and (21) in reference 3). 

As a simple application of Eq. (2), consider a cone decelerating 
Taking the equa- 
U(t) = 


where 


uniformly from supersonic to sonic velocity. 
tion of the leading edge as — = cr + br?, we have 
= 2ht; 
6 = semiangle of cone 

For V b X body length/c << 1, 


found 


the velocity is sonic at ¢ = 0. A(x) = 7wé*x? 


an approximation is easily 


2rd*c(x + ct) 


2) (x — ct) 


Sx, t) = 
x(x, t) = (1 
x(x, t) = eV (x+ct)/b+x+ ct 


At sonic velocity the nose is at x = 0), and the pressure coeffi- 


cient including square terms is 


Cp = 28% log [V2c/a(bx)'/*] — (82) (4) 


The logarithmic singularity as the deceleration » — O corre- 
sponds to the growth of pressure proportional to log time, a result 
of linearized theory at M = 1. Forany } ¥ 0 linearized theory 


gives a finite resistance at WM = 1. 
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Remarks on ‘‘Shock-Wave Effects on the 
Laminar Skin Friction of an Insulated Flat 
Plate at Hypersonic Speeds’’! 


M. Z. v. Krzywoblocki 

Professor of Gasdynamics and Theoretical Aerodynamics, University 
of Illinois, Urbana, III 

July 2,1953 


ACH PROBLEM IN ENGINEERING IN PARTICULAR, and in mathe- 
matical physics in general, may and can be reduced to the 
fundamental problems in these two fields: mathematics and 
theoretical physics. As the writer has shown,? the problem of 


the boundary layer reduces to the two fundamental problems: in 


J 
2 
o 


FORUM 


mathematics (existence and, possibly, uniqueness theorem) and 
in theoretical physics (the proper equations describing adequately 
the piiysical phenomenon in question). There is not much ad- 
vantage in going into details such as the interaction between a 
shock and the boundary layer without first solving the funda- 
mentals. The problem of interaction will technically introduce 
only some changes in the boundary conditions, since, from the 
mathematical standpoint, in a boundary value problem with a 
fixed differential system (partial or ordinary), the only changes 
that can occur are in the boundary conditions. If these boundary 
conditions are consistent with the boundary value problem in this 
or in that case, mathematically there will be no changes and the 
question reduces to some formal calculations. 

If they are inconsistent then a completely new boundary value 
problem must be attacked and solved. It is also known that, at 
the present status of knowledge in mathematics, one may be 
able to prove some existence theorems but not that of uniqueness 
with respect to the particular differential system and boundary 
value problem in the paper in question. Hence, the speculations 
of the authors on the uniqueness of the solutions (page 353)! 
Turning to the funda- 
that we 


simply cannot be seriously accepted. 
mental problem in theoretical physics; any statement 
can accept the Navier-Stokes equations as the universal system 
able to solve all of our problems of the boundary layer from zero 
altitude up to an arbitrary distance above the earth simply leads 
toanabsurdum. The solution of the problem of validity of vari 
ous forms of equations of motion and energy on various levels in 
various intensities of density of a gas is the most urgent problem 
of today. The kinetic theory of gases of a monatomic gas, even 
with all its deficiencies, applied to rarefied air may give errors 
of a smaller order than the Navier-Stokes equations of continuum 
mechanics applied to that range of rarefied gases. Of course, all 
such statements must be accepted cautiously, up to the moment 
when more data are available. Even now some results from the 
theory of shocks seem to indicate that such statements are true. 
The authors! presented some considerations on hypersonic flow 
going into the fundamentals. Hence some remarks 
The fundamental equations are the Navier 
Hence, logically, the results are valid only 
Some cases may 


without 
seem to be in order. 
Stokes equations. 
in the domain where these equations are valid. 
be discussed: (a) Certainly they are valid in the range of high 
density gases, but in this range there is no sense in talking about 
Mach Numbers of the order of 15, say. (b) Since the possibility 
of reaching such high Mach Numbers is in the range of rarefied 
gases, there may be a level at which the Navier-Stokes equations 
are approximately valid; but to show that, one has to show that 
any other system of equations, based on the kinetic theory of 
gases, say, containing higher order terms, gives an error greater 
than that of the Navier-Stokes equations. (c) There may be 
an altitude at which the equations based upon the kinetic theory 
of gases of a monatomic gas will give less error due to higher order 
terms than Navier-Stokes equations; in this case one has to use 
the kinetic theory of gases as the fundamental concept. (d) 
There may be the region in which the kinetic theory of polyatomic 
gases must be used. As is seen, to understand the results derived 
by the authors, a thorough discussion, perhaps beyond present 
possibilities due to the lack of sufficient data and information, 
must be presented. This is completely lacking. Of course, the 
higher order terms due to the kinetic theory of gases, if taken into 
account, may significantly change the numerical results of the 
authors. Another question is the thickness of the shock wave. 
In their discussion the authors use the calculations by A. E 
Puckett and H. J. Stewart based upon the Navier-Stokes equa- 
tions. But it is well known that the equations of continuum me- 
chanics are not adequate enough to supply information on the 
thickness of a shock wave. The only results comparable with 
tests were obtained by means of the kinetic theory of gases.* 

Since, moreover, the thickness of a shock wave increases con- 
siderably in the range of hypersonics, it seems that again this 
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problem was presented superficially in the paper in question. 
One has to understand deeply that the problem of the boundary 
layer in hypersonics cannot be solved merely by a simple exten- 
sion of Prandtl’s 50-year-old boundary-layer notion. A thorough 
construction of new fundamentals must be accomplished, and 
from this standpoint the address given by Professor Hirschfelder 
during the Third Midwestern Conference on Fluid Mechanics in 


Minneapolis was of great significance. 
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An Investigation of 





Ludwig Prandtl, 1875-1953 


(Concluded from page 779) 


Of his other major contributions, we have space only ~* leader of men. 


to mention his mixing lengta theory of turbulence, his 
work on heat transfer in moving fluids, his pioneer wind- 
tunnel designs, the Prandtl soap-film analogy for the 
study of torsional deflections and stresses in prismatic 
bodies, and his papers on creep and flow, on the theory 
of failure, and on the penetration of rigid bodies in a 
plastic medium. 

Those who had the privilege of friendship with 
Prandtl knew him as a modest, friendly, unassuming 
man, far from one’s preconceived notion of a great 


It is said that his lectures were not suit- 
able for beginners and that he was always surprised 
that others did not see the relationships so clearly as he. 
He taught by association and example. He gave much 
credit to his coworkers and deprecated his own contri- 
butions to aeronautical practice. His work will live 
long not only in the impress of his concepts on aeronau- 
tical practice but in the ever widening influence of the 
scientists and engineers who have felt the challenge of 
HuGu L. DRYDEN 
Director, NACA 


his personality and work. 











